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ABSTRACT 


A mathematical model is developed to analyze the seismic 
response of a drydocked vessel in three degrees of freedom; 
vessel rotation about its keel and vessel horizontal and 
vertical translations relative to the drydock cradle. Data 
from eleven actual vessel-drydock systems and the time 
acceleration history of an earthquake are implemented to 
predict vessel three degree of freedom response during an 
earthquake. Vessel seismic response and the resultant drydock 
forces are compared to vessel-drydock system failure criteria 
to determine stability of the system during earthquakes. The 
three degree of freedom vessel response model is compared to a 
one and two degree of freedom vessel response models and a 
model in which seismic loading is simulated by a single static 
force. The three degree of freedom vessel motion is shown to 
be the most accurate method for analyzing vessel-drydock 
system failure criteria. 
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1.0 INTRODUCTION 

1.1. Background 

In the design of a drydock for military vessels, an 
important factor is the environmental-generated forces that 
are found in the drydock cradle which supports the ship. 
Earthquakes and the wind are considered the two most 
significant environmental loading factors. This thesis 
explores the seismic loading aspect of the problem. 

Regardless of its location, any U.S. Navy ship that is 
nuclear powered and in drydock is considered to be in a high 
seismic risk area. Currently the Navy defines a high seismic 
risk area as one in which earthquake forces be approximated by 
a steady horizontal force of 0.2 g times the vessel mass 
acting at the center of gravity. The drydock can be analyzed 
to see if it can withstand this quasi-static loading. 

1.2 Previous Work 

In 1981 B.V. Viscomi studied the seismic response of a 
drydocked submarine using a single degree of freedom model 
[14]. This analysis assumed that the vessel was allowed to 
rotate about the keel. Viscomi used the time history 
acceleration record of the North-South component of the 1946 
El Centro, CA Earthquake, California Institute of Technology 
processing scheme. The vessel was analyzed for system failure 
(l1.e. the vessel lifting off of a row of side blocks). 

A thesis by C.F. Barker in 1985 [1] used a two degree of 
freedom model to study seismic response. Employing the 1946 


El Centro Earthquake, Massachusetts Institute of Technology 
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standard processing scheme, Barker analyzed rotation about the 
keel plus horizontal relative translation between the vessel 
and the drydock cradle perpendicular to the vessel longitudal 
axis. The two degree of freedom model permitted the analysis 
of drydock cradle failure: block sliding and tipping, which 
were not possible to analyze in Viscomi's one degree of 
freedom model, along with block liftoff. 

1.3 Contributions of This Thesis 

This thesis develops a three degree of freedom model of a 
vessel in drydock. Motion will be permitted in three 
directions: rotation about the keel and relative 
translations, horizontal and vertical, between the vessel and 
the drydock cradle. The three degree of freedom model will 
make possible the analysis of drydock cradle failure: 
additional block crushing forces due to vertical motion, which 
was not conceivable in the previous one and two degree of 
freedom models developed by Viscomi and Barker respectively, 
pong with block sliding, tipping and liftoff. 

The seismic input for the three degree of freedom model 
is the time acceleration history of the North-South component 
of the 1946 El Centro, CA Earthquake, Massachusetts Institute 
of Technology standard processing scheme. The MIT version of 
the El Centro Earthquake acceleration record is described in 
reference [9] , and is displayed graphically in Figure 1.1. 

1.4 Outline of This Thesis 

Section 2 of this thesis describes the vessel-drydock 
system for the reader and introduces terms which are used 


throughout the thesis. Section 3 describes the various 
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failure mechanisms found in the vessel-drydock system which 
can occur. 

Section 4 examines the modelling of the vessel-drydock 
system under seismic loading. Section 4.1 describes the 
quasi-static loading method, and Section 4.2 develops the one 
degree of freedom equation of vessel motion model. Section 
4.3 develops the two degree of freedom equations of vessel 
motion model. Finally, Section 4.4 develops the three degree 
of freedom equations of motion model. 

Section 5 discusses the system parameters that are 
required for the mathematical model found in Section 4. 
Included in this Section is the modelling of block stiffness. 
Eleven typical vessel-drydock system configurations are 
studied to be implemented into Section 6. 

Section 6 examines methods of evaluating the one, two and 
three degree of freedom models found in Section 4. Section 
6.1 evaluates the linear equations of motion models using 
modal analysis method. Section 6.2 evaluates both linear and 
non-linear equations of motion models using a fourth-order 
Runge-Kutta numerical analysis scheme. Section 6.3 explains 
the response spectrum method of determining the maximum value 
of vessel seismic response. Section 6.4 determines system 
response uSing quasi-static force method. 

Section 7 describes the development and testing of the 
one, two and three degree of freedom vessel motion computer 
programs, and contains the results for a seismic analysis of 
several actual vessel-drydock systems. 


Section 8 summarizes the response predictions for several 
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actual vessel-drydock systems obtained with quasi-static, one 
degree of freedom, two degree of freedom and three degree of 
freedom analysis methods, and conclusions are drawn. 

Appendix 1 is an example listing of the fourth-order 
Runge-Kutta computer program. Appendix 2 is an example 
listing of the modal analysis computer program. Finally, 
Appendix 3 describes the modal analysis of the two and three 
degree of freedom vessel-drydock systems and predicts maximum 
system response using the response spectrum method of analysis 


and modal participation factors. 
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2.0 THE VESSEL-DRYDOCK SYSTEM 

2.1 General Description 

In order to perform inspection, maintentance and repair 
on the outer hull of a ship, it is necessary to place the 
vessel ina drydock. A drydock is a concrete encasement built 
into the earth. At one end of the drydock is a flood gate to 
allow the entry and departure of the vessel. A cradle, formed 
by drydock blocks, is built upon the drydock floor in order to 
Support the ship once the water is pumped out of the 
encasement. The drydock must be able to hold the weight of 
the vessel and cradle. 

Due to the various hull configurations associated with 
the different classes of ships, the cradle component of the 
drydock-cradle-vessel system is varied to suit a particular 
docking situation. The cradle accomplishes the transfer of 
vessel weight to the drydock floor, provides stable support to 
the vessel and allows access to the vessel hull. This cradle 
is usually constructed of timber, concrete, or a timber- 
concrete composite. 

2-2 system Component 

The primary components of the drydock-cradle-vessel 
system are shown in Figure 2.1. This section defines terms, 
related to docking, used in this thesis to describe the system 


components. 
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Blocks - The units, consisting of timber, concrete, steel 
and other materials, which together make up the cradle. 


Cradle - A framework of blocks which supports the vessel 
when the drydock is dewatered. 


KG - The distance between the vessel baseline (ship's 
keel) to its vertical center of gravity. 


Keel Blocks - The center blocks, directly beneath the 
vessel's keel. 


Pier - A column built of blocks that extends from the 
Sittipenudseetoucdock Floor. 





Side Blocks - The blocks located to the right and left of 
the keel blocks. 


Ton - A long ton, 2240 lbs. 


g- Acceleration of gravity, 32.2 ft/sec?. 
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Vessel-Drydock System 
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3.0 VESSEL-DRYDOCK SYSTEM FAILURE MECHANISMS 

There are four failure modes that the vessel-drydock 
system can exhibit: block crushing, block sliding, block 
overturning and block liftoff (1.e. the liftoff of vessel from 


a row of blocks). 


3.1 Block Crushing® 


A block can support compressive stresses linearly until 
the blocks' proportional limit is exceeded. The proportional 
limit is a material property of the block and is the maximum 
compressive stress at which stress is still linearly 
proportional to strain. If the compressive proportional limit 


is exceeded, the block is considered to have failed. 


3.2 Block Sliding® 


Due to friction at the interfaces between the blocks and 
the vessel and between blocks with other blocks, drydock 
blocks will have the tendency to resist sliding when subjected 
to vertical and horizontal loads. This is true when in the 
absence of mechanical fasteners. 

First consider the blocks that form the keel pier. 
Figure 3.1 shows that the resistance of the keel pier to 
sliding is 

H =H jV (Biock—-bleck interface) 
and H =Ub5V (vessel-block interface) 
where 

H = horizontal resisting force 


= coefficient of friction, block-block interface 


y 
~~ 
| 


U5 = coefficient of friction, block-vessel interface 


vertical load on the blocks. 


= 
ll 
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Failure of the keel pier due to sliding will occur when 
the ratio of horizontal force to the vertical force is greater 
than the lowest coefficient of friction found in the block. In 
other words, sliding occurs when 

>(U} Or Uo). 

As with the keel piers, the resistance of the side piers 
to sliding will depend on the interface which offers the least 
resistance, the block-vessel interface or the block-block 
interface. However, this case is more complicated due to the 
complex geometry of the side pier as shown in Figure 3.2. 

The block-vessel interface of the side pier is examined 
and shown in Figure 3.3. An arbitrary force F at some angle 
from vertical is applied to the face for the side block cap. 
At this interface, the horizontal and vertical block force 


reactions are: 


Normal Force = H sind + V cos¢ 


Tangential Force H cos¢ - V sind 


where H = Horizontal Reaction 
V = Vertical Reaction 
@ = Block Inclination Angle. 


Using the formula for horizontal sliding resistance, H =U 
V, previously introduced in the keel pier case and the 
above equations yields 

H cos¢'- V sind = U5 (H sind + V cos). 


Rearranging the above equation gives 


H _ ¥4 cosd + sind 


V cos¢ - U5 sind 
(opm 

a) ay 

vy ~ 92 
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where U5 cosd + sind 


M2 ~ “cosé? - Hsin? 


For the block-block interface of the side pier, sliding 
resistance is found by taking the vertical and horizontal 
components of the applied force, F, and using them in the 
formula H =u; V. Like in the keel pier, failure of the side 
pier due to sliding will occur when the following is satisfied 

e > (, Orly!'). 

There is a critical angle o where no slippage occurs. 
Referring again to Figure 3.3, this angle can be calculated by 
comparing the maximum tangential force Fy may the interface 


can have without slipping and the corresponding normal force 


Fr- The following relationship holds 


Ft max = -EFn 


Finally, 

a = are tanu. 
Since the applied force F acts in a straight line through all 
blocks to the ground, the angle & for each interface must be 
calculated. If force F is applied outside of these angles, 


slippage will occur. 


3.3 Block Overturning!1 


The third failure mechanism is the overturning of a block 
due to an applied force. The line of action of this force 
must fall within the middle one-third of the base of the 
block, as shown in Figure 3.4, or the block will tip over. 

In order to show the region that the block will remain 


upright, the inclined force F must be broken into a transverse 
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component Fy, and a vertical component F,. Superposing the 
bending and axial stresses caused by F, and Fy, gives the 


resultant stress at any point in the block as 


Th gee 

Cac ieeeey A 

where L = block height 
A = base area, here bh 


I, = base moment of intertia, here 3/12, 
The minimum stress point on the bottom of the block 
occurs at point M (Figure 3.4), a distance -h/2 from the z 


axis. Hence, 
Ee 
el bh h bh 


When °,, at M is negative, compressive stress is present and 
the block is stable. When °o, at M is positive, tensile 
stress is present and the block will overturn since no 
fasteners hold block to the ground to develop tensile stress. 
The limit of block stability is when 
(Sx) m = 0. 
Thus, the following equation holds 
ae = oF. 

This condition exists when force F is applied at an angle such 
that the line of force lies within the one-third of the base 
of the block. If the line of force is outside this region, 
the block will tip. 

3.4 Vessel Liftoff! 

The fourth and final failure mode occurs when the vessel 
breaks contact with either the side or keel piers. This 


failure occurs when the dynamic deflection of a row of blocks 


is equal to, or exceeds the average static deflection of the 
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blocks. The static deflection of blocks is given by 


a 
Ss 2K a2 UX 
SV kv 
where W = submarine weight 


side pier vertical stiffness 


Ksv 


keel pier vertical stiffness, 


Kv 


The dynamic deflection can be due to vertical displacement y, 


rotational displacement 9%, or a combination of both. 
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FIGURE 3.1 
Keel Pier Forces 





FIGURE 3.2 
Side Pier Forces, General 
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FIGURE 3.3 
Side Pier Forces, Block-Vessel Interface. 
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FIGURE 3.4 
Block Tipping Forces 
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4.0 MATHEMATICAL MODELS OF VESSEL-DRYDOCK SYSTEM 

This portion of the thesis discusses the modelling of the 
drydocked submarine subjected to a seismic load. Models will 
include a quasi-static method, one, two and three degree of 
freedom mathematical models. The models developed in this 
section will be analyzed in Section 6 and 7 to predict the 


vessel seismic response. 


4.1 The Quasi-static Method of Modelling Seismic Response 


fala ee Approach? © 


Current U.S. Navy design method for submarine's seismic 
response will be examined in this section. The submarine is 
described as a rigid cylindrical body with its weight evenly 
dispersed longitudinally. 

The quasi-static force method replaces the earthquake 
motions by a force corresponding to the vessel mass time 
0.2 g. This force is horizontally applied to the submarine's 
center of gravity in the tranverse direction. The drydock 
blocking system is determined by the evaluation of the loads 
generated by the quasi-static force. 

4.1.2 Force Equations 

The application of the static force to the center of 

gravity of the submarine in drydock is represented in Figure 


4.1. The seismic overturning moment, Ms, is defined as 


Me = (A/g) (a) (KG) (2240) ft-lb 
where 
A = vessel displacement in long tons 
g = acceleration of gravity 
a = vessel's center of gravity acceleration. 
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Current U.S. Navy design practice states that the acceleration 
of the submarine's center of gravity, a, should be set equal 
to 0.2 g. Hence 
Ms = 448(A) (KG) ft-lb 
After determining the seismic overturning moment, the 


number of sideblocks required on one side of the vessel is 


Ne 





where N = # of side blocks required in the row 
A = contact area of a block in square inches 
Sp = proportional limit of the block cap in lb/in? 


L = distance between centerline of keel and side 
blocks, in feet, as shown in Figure 4.1. 


4.1.3 Vessel Response in the Quasi-static Load 
Method 


The application of a force applied at a vessel's center 
of gravity will generate reactionary side block forces in the 
vertical direction to oppose the seismic moment, as shown in 


Figure 4.2. Summation of moments about the keel yields the 


equation 
IM = Me - L(Fer + Fmr) + L(Fe] - Fmi) = 0 (Aran) 
where 
F. = static pier forces (left and right) 
Fy, = pier forces due to the applied moment (left and right), 
Due to symmetry in the vessel-drydock system, the 


applied moment resistive forces in the side blocks are equal 
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in magnitude but opposite in direction. Thus, Equation 4.1.1 


simplifies to the form 
M. = 2L(F,) 
where Fy = Fri = Fur- 


The force F, is equal to the side pier compression 


displacement, 6m, times its vertical stiffness, Kyg. 


compression displacement can be expressed as 
o, = Lb since 


where 9 = angle or vessel rotation about its keel. 


CO Shee 


THs 


(ae) 


Assuming the rotational angle is small, Equations 4.2.1 and 


4.1.3 are combined to yield 


2 a 
2L° K 
SV 


4.2 One Degree of Freedom: Rotational Response 


4.2.1 Approach 


(4.1.4) 


A one degree of freedom mathematical model of a submarine 


in drydock is developed in this section. The submarine is 


identified as a rigid cylindrical body with an even 


longitudinal weight distribution. The system damping and 


stiffness are provided by the drydock blocking arrangement. 


System excitation is due to seismic ground accelerations. 


Other assumptions are 


- The keel and side piers remain vertical during the 


ground motions, 


- no slippage between the block cradle and drydock, and 


- the keel and side piers bases accelerate at the same 


rate as the drydock (ground). 


The idealized model of the vessel-drydock system is shown 
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in Figure 4.3. In this system, vessel motion is restricted to 
rotation about the keel. The side piers are modelled as 
vertical due to typically small angles or inclination of the 
tops of the side piers. Only vertical displacements in the 
side piers are considered. 

The one degree of freedom mathematical model predicts and 
analyzes vessel rotational response about its keel caused by 
seismic ground accelerations normal to the submarine's 
longitudinal axis. The resultant angle or rotation, expressed 
in radians, and side pier forces, generated by the rotation, 
are compared to applicable pier failure criteria found in 
Section 3, to see if the blocking arrangment will remain 
intacted. 

4.2.2 One Degree of Freedom Equation of Motion Model 
For the one degree of freedom case, the keel is the 
origin with clockwise rotation, §&, positive. The simplified 
vessel-block system showing coordinate © at rest and when 
excited are shown in Figures 4.3 and 4.4 respectively. The 
moment equilibrium equation about the keel (origin) is: 


IMy = I,6 + MKGxg - MKGyg siné 
= B/2(Fisy - Fiagy) - B/2(Fosy + Foqv) + WKG sing (4.2.1) 


where I, = Mass moment of inertia of vessel about the keel 


Vessel Mass 


M 
Fisy,Faosy = static side block forces in the vertical direction 
(Figure 4.4) 


Figyv,Foav = dynamic side block forces in the vertical 
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direction (Figure 4.4) 


W 


vessel weight (Mg) 


and Xg,Yg = seismic ground accelerations in their respective 


adirections 


This equation can be simplified by the following: 


1.) The left and right side block rows have the same 
vertical spring constant, Ks,y, have the same static 
deflection and equal but opposite dynamic deflection due 
to symmetry, i.e. 
Fisv = Fasv = Fsv (4.2.2) 
and | Fygy!= | Foqv!= Fay: (ee as) 
2.) The dynamic side block force, Fgy, as seen from Figure 
4.5 equals to a force due to the modelled spring 
Gisplacement plus a dissipative force due to system 
damping. In other words, 
Fay(spring) = (B/2)Ksy siné 

and Fg,(dissipative) = C(B/2) 6 


or when combined, 


Fay = B/2 (Key sin® + C8) (4.2.4) 


where Key = vertical spring constant of the side blocks 


C = vertical damping coefficient of the side blocks 


Substituting Equations 4.2.2, 4.2.3, and 4.2.4 into 


Equation 4.2.1 and rearranging components yields 


I}, + (B22) C8 + [(B“/2)Ke, - W KG] siné 


= -MKG Xg + MKG Yg sin @ (4.2.5) 


An additional simplification is that since Sis very small, 


Sale) ses & (ele 


The final simplification is that damping is expressed as a 
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PoctlonsOor critical damping. To do this, allow 
(Ee72)Ge= & (Co, = 26 1,6, 


where — = viscous damping factor B-C/2C...) 


Cor critical damping coefficient = 21}, 


U) 


n system undamped natural frequency. 
Now, Equation 4.2.5 simplified to 

1,6 +2 Ip p68 + [(B2/2)Key - WKG]O 

= -MKG xg + MKG yg 6 (4.2.6) 

The non-linear tern, MKG yg¢, found in the equation of 
motion will be removed at this time to linearize Equation 
4.2.6 since 6 is assumed to be small. The effects of the 
non-linear term will be evaluated in Section 6 and 7. 

The final form of the linearized one degree of freedom 
equation of motion is 

m6 + C16 + kj8 = -mg Xg (4 20a) 
mere MM) = ly 
Cy = 26 Ip 


Kk) = (B“/2)Key - W KG 


M> = M KG. 


imematrix forn, 


- 
4. 
J 
cD 
T 


To simplify the explanation, the matrices will be redefined as 


[Al{y} + [B]{y) 
© m ltt) @) 

[A] = : IB} =| + 

_ my Cy oO ky 
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{E(t) } (Arana) 


where 





Equations 4.2.6 and 4.2.8 will be solved in Sections 6 and 7. 


4.3 Two Degree of Freedom: Rotational and Horizontal 
Translation Response 


4.3.1 Approach 


This section develops a two degree of freedom 
mathematical model of a drydocked submarine. As in the case 
of the one degree of freedom (Section 4.2.1), the model 
assumptions still hold, i.e. a rigid cylindrical body, keel 
and side piers remain vertical, no slippage between the block 
cradle and drydock, and the bases of the keel and side piers 
accelerate at the same rate as the drydock (ground). 

The idealized model of the vessel-drydock system is shown 
in Figure 4.6. In this system, vessel motion is restricted 
to rotation about the keel, and in horizontal transverse 
translation relative to the keel and side pier supports. The 
approach for determining vessel-drydock failue is identical to 
the one degree of freedom case except for the addition of 
sliding and tipping failure modes. 

The two degree of freedom mathematical model predicts and 
analyzes submarine rotational and transitional response caused 
by seismic ground accelerations normal to the vessel's 
longitudinal axis. The rotational response is in radians, the 
translational response in inches, and both responses are 
expressed as a function of time. The keel and side pier 


forces are compared to applicable pier failure criteria found 
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ime SeCeion 3. 
4.3.2 Two Degree of Freedom Equation of Motion Model 

As in the one degree of freedom case, the keel is the 
origin for the two degree of freedom system with clockwise 
rotation, 8, is positive. The linear relative translation 
coordinate, x, is defined with respect to vessel and drydock 
translation relative to the ground. Defining u as the 
position of the submarine keel relative to the ground, then 
the following relationships hold: 


u xX + Xq 
and u= xX + Xg. (Aen) 

The simplified vessel-block system showing coordinates (6 and 
xX) at rest and when excited are shown in Figures 4.6 and 4.7 
respectively. 

A balance of forces in the x direction, as shown in 
Figure 4.7, yields the first equation of motion 

DFy = Mu + MKG6= (Fish - Figh) - (Fash + F2an) 


- (F3sp + F3qnh) (4.332) 


where Fysh, Fosh, F3spn = Static block forces in the horizontal 

direction at their respective 
position 

Figh, Foaghr F3qh = dynamic block forces in the 
horizontal direction at their 
respective position. 

The above equation can be simplified by the following: 

1) The left and right side block rows have the same 


horizontal spring constant, K.g,, have equal but opposite 


static horizontal deflection and the same dynamic horizontal 


-34— 





deflection due to symmetry, i.e., 

Fish + F2sh = 0 

Fidh = F2dh - (4.3.3) 
Also, the keel block rows do not experience any net static 
horizontal deflection since the inclination angle of the block 
cap is zero. Hence 

Fach = 0. (Ace) 

2) The dynamic block forces, Fg, as seen from Figure 4./ 

equals to a force due to the modelled springs plus a 


dissipative force due to system damping. In other words 


Fah (spring) = Fidh (spring) + F2dah(spring) + F3dah (spring) 


KophX + Ken X + KypX 
Fah (spring) = (2Ksh + Kxn)x 
Likewise 
Fah (dissipative) = Cx* 
SO 
Fah = (2Kns + Knz)X + Cyx (Aran 5) 
where Kop, = Side pier stiffness 
Kyph = Keel pier stiffness 
Cy = system horizontal damping coefficient. 
Subscicucing Equations 4.3.1, 4.3.3, 4.3.4, and 4.3.5 
into Equation 4.3.2 and rearranging components yields 
Mj) X + mpo9 +c, X+ ky, X = -my] Xg (acre) 
Paere Mm), = M 
ees KG 
Cc) = Cy = system horizontal damping coefficient 
ki1 = (2Ksp + Kyp) = System horizontal stiffness 


Summing the moments about the origin in a similar fashion 
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as the previous degree of freedom case, as shown in Figure 
4.7, yields the second equation of motion, 


IMy = Ip + MKG u - MKG yg sino 
= B/2(Fisy - Figy) - B/2(Fosy + Foay) + WKG sin® (4.3.7) 


The second equation can be reduced as it was in the one degree 
of freedom case to 
1,8 + MRG X + C55 + (B2/2 Key - WKG)9 
= MKG Xg + MKG Yg°- (Ames: 8) 
Once again, the non-linear term, MKG Yg% found in 
Equation 4.3.8 will be removed leaving a linearized equation 
Since @ is assumed to be small. The effects of the nonlinear 


term will be evaluated in Sections 6 and 7. 


Now the system of equations for the two degree of freedom 
vessel-drydock system is as follows, 
my, X + myo + cz, X + kyy xX = -m} Xq (4.3.9a) 
and 
m5 © staal x a 6158 + ko99 = -m9] Xq (4.3.9b) 
Bere im) = M 
M21 = M12 = MKG 
m2 = I, 
Cc) = system horizontal damping coefficient 
Cy = system rotational damping coefficient 
Kki1 = 2Ksp + Kxh 
B2/2 Ke, - WKG. 


k22 


These two equations of motion are coupled in the mass times 
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acceleration term, which is known as inertial coupling. Also 
the equations are coupled in the damping term, because the 
system's two natural modes have translation and rotation 
involved in them. This model uses 5 percent of the critical 
damping coefficient as system damping. To evaluate damping, a 
modal analysis for the system is performed. This analysis is 
performed in Appendix 3. Substituting damping coefficients 
into Equations 4.3.9a and 4.3.9b, the equations of motion can 
be expressed in matrix notation as Equation 4.3.10: 

(M] (yt) + (C} (84) + CK] (vt) = (BN(t)) (4.3.10) 


where 
{y'} = Resnonse vector = {% 


m Meee 
[M] = Mass nacrix| 11 | 
a De 
Ss Cc 
fel = Damping Matrix = 11 =| 
ie 0 
Ky4 O 
[K] = Stiffness Matrix = 
7 O K55 


cls: 
late =—seismic Forcing Vector = S| 

-m 

2 Ss 


To simplify the explanation, the matrices will be redefined as 


[A] {vy} + [B] (vy) = {E(t)} (453-11) 


Oo 6[M] 
UES SS 
[M] [C] 


where 
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[B} 


{y} = 


bse 


tl 

(eee ee 
Foxe} 
e9 

- © 
ct 
ey 
ee ee” 

e 


Equations 4.3.6, 4.3.8, and 4.3.11 will be solved in Sections 


6 and 7. 


4.4 Three Degree of Freedom: Rotational, Horizontal, 
and Vertical Translation Response 


4.4.1 Approach 


This section develops a three degree of freedom 
mathematical model of a drydocked submarine. The idealized 
model of the vessel-drydock system is shown in Figure 4.8. [In 
this system, vessel motion is restricted to rotation about the 
keel, and in horizontal transverse and vertical translations 
relative to the keel and side pier supports. The approach for 
determining vessel-drydock failure is identical to the two 
degree of freedom case. 

The three degree of freedom mathematical model predicts 


and analyzes submarine rotational, horizontal and vertical 
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transitional responses caused by seismic ground accelerations 
normal to the vessel's longitudinal axis. The rotational 
response is in radians, the translational response in inches, 
and all three responses are expressed as a function of time. 
The keel and side pier forces generated by the responses are 
compared to applicable pier failure criteria found in section 
3 to check blocking stability. 

4.4.2 Three Degree of Freedom Equation of Motion Model 

As in the two degree of freedom case, the keel is the 

origin for the three degree of freedom system with absolute 
rotation, ©, and translation coordinate, x, is the horizontal 
vessel motion relative to the drydock. The linear relative 
translation coordinate, y, is defined with respect to vessel 
and drydock translation relative to the ground. This 
coordinate y describes the vertical motion of the vessel 
hull relative to the drydock and the drydock blocks. Defining 
v as the position of the submarine keel relative to the 


ground, then the following relationships hold; 


Yor v= Yo 

v=ytyYqg 
and veytYg (4.4.1) 
Likewise, u=x + Xg (aa) 


The simplified vessel-block system showing coordinates (x, y 
and 6) at rest and when excited are shown in Figure 4.8 and 
4.9 respectively. 

As developed in Section 4.3.2, the first and third 
equations of motion in the three degree of freedom model are 


Mx + MKGO + Cyx + (2Kgh + Kyn)X = -MXg (4.4.3) 
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and 
I, 9 + MKG x - MRGy® + c® 
+[ (B*/2)Kgy - WKG]o = -MKG xg + MKG ygé . (4.4.4) 
A balance of forces in the y direction, as shown in 
Figure 4.9, yields the second equation of motion 
-Fy = Mv = (Fisy ~ Fiav) + (Fasv ~ F2av) 


The above equations can be Simplified by the following: 

1) The vessel weight, W, must be equal to the summation 
of static block forces so that static equilibrium holds true, 
ae 

Fisy + Fosy + Facy = W (4.4.6) 

2) The dynamic block forces, Fq,, as seen from Figure 
4.9 equals to a force due to the modelled springs plus a 
dissipative force due to system damping. These dynamic block 
forces are only a function of vertical displacement, y, and 
velocity, y. The formulation of these forces is identical to 
the procedure used in Section 4.3.2 and will not be repeated. 


The final form of the dynamic block forces, Fyy, 1s 


Fay = (2Ksy + Kyy)y + Cy y (4.4.7) 


where Kgy = side pier vertical stiffness 


Kyy = keel pier vertical stiffness 
Cy = system vertical damping coefficient. 
Substituting Equations 4.4.1, 4.4.2, 4.4.6 and 4.4.7 into 
Equation 4.4.5 and rearranging components yields 


My + Cyy + (2Ksy + Kyy)y = -MYg. (4.4.8) 


The above equation along with Equations 4.4.3 and 4.4.4 


-40- 





formulate the three equations of motion for the three degree 
of freedom model. 

The non-linear terms, MKGy® and MKGy,°, found in 
Equation 4.4.4 will be removed at this time leaving a 
linearized equation since © is assumed to be small. The 
effects of the non-linear terms will be evaluated in Section 6 
and 7. 

Now the system of equations for the three degree of 


freedom is as follows; 


M1] x + M13 ° (24) x + kj] X¥ = -M)] Xg (4.4.9a) 
Moony Co Y + ko2 Gage ato Yg (4.4.9b) 
m336 + m3) X + cz 6 + k339 = -m3) Xg (4.4.9¢) 


mnere M)) = Moo = M 


M13 = M3) = MKG 


m33 = Ik 

C) = system horizontal damping coefficient 
Co = system vertical damping coefficient 
C3 = system rotational damping coefficient 


Ki1 = 2Ksn + Ken 

Kko2 = 2Ksy + Ky 

k33 = B“/2 Key - WKG. 
The first and third equations of motions, Equations 4.4.lla 
and 4.4.llc, respectively, are coupled in the mass times 
acceleration term, which is known as inertial coupling. Also, 
these equations are coupled in the damping term because two of 
the system's three natural modes have horizontal translation 
and rotation involved in them. The second equation of motion, 


Equation 4.4.11b, is uncoupled from the rest of equations. 
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Hence, the last natural mode of the system only depends on 
vertical translation. 

In order to introduce damping, this model uses 5 percent 
of the critical damping coefficient as system damping. To 
evaluate damping, a modal analysis for the system in question 
must be performed. This analysis is performed in Appendix 3 
System damping coefficients can be found using this modal 
anlaysis technique, but for now, the coefficients will be left 
in terms of C11, ©C13, C22, ©€31 and c33. Substituting these 
coefficients into Equations 4.4.lla, 4.4.11b, and 4.4.llc, the 
equations of motion can be expressed in matrix notation as 
Bouation 4.4.12: 

(M] (y'} + [C) (vt) + [K] (vt) = (Et (t)} (4.4.10) 


where 


x 
oe = RASS oMae Vetekeene = P| 
8 
1 He 
[M] =eMacss i flatrix = | 0 Mo 5 O 
a aed 33 
Kia O O 
[K] =——SEletness Maurix = Oo K55 O 
O O K33 
wee — Seismic Forcing Vector = “M11 Xo 
oe Yq 
-M34 see 
C O Cc 
iG] = Damping Matrix = tn 13 
O C55 0 
Sole” <a 
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To simplify the explanation the matrices will be redefined as 


[A] {vy} + (B] {vy} = (E(t)) GaAs) 
where 
[6 [M} 
[A] = 
[M] [Cl] 
-{M] 0 
ua = 
Coa K | 
ej 
{y} 
iA ES 
we 
O 
{E(t)} = 
Lee) | 


MeWations 4.4.3, 4.4.4, 4.4.8, and 4.4.11 will be solved in 


Section 6 and 7. 
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Vessel-Drydock System with Applied Quasi- static Force 
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Idealized One Degree of Freedom Vessel- Drydock System at Rest 
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FIGURE 4.4 
Idealized One Degree of Freedom Vessel- Drydock Systen, 
Excited 
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Dynamic Side Block Force 
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Idealized Two Degree of Freedom Vessel-Drydock System at Rest 
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Idealized Two Degree of Freedom Vessel-Drydock System, 
Excited 
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FIGURE 4.8 
Idealized Three Degree of Freedom Vessel-Drydock System at 
Rest 
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FIGURE 4.9 
Idealized Three Degree of Freedom Vessel-Drydock Systen, 
Excited 
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5.0 SYSTEM PARAMETERS 

In order to study the response of a vessel-drydock system 
to seismic loading, eleven typical configurations have been 
selected for implementation into the one, two and three 
degrees of freedom models developed in the previous section. 
This section explores these eleven systems along with the 
modeling of the pier blocks for stiffness. 

5.1 Vessel-Drydock System Parameters 

Eleven drydock system configurations have been chosen for 
analysis and are defined in Table 5.1. The vessel-drydock 
parameters, corresponding to each system, were taken from the 
appropriate NAVSEA drawings. 

ABE Oo 


Vessel-Drydock Configurations 


Block Longitudinal NAVSEA 

System Hull Type Block Spring Drawing 
il 616 Composite Seite 845-2006640 
2 616 Composite ie see 845-2006640 
3 616 Timber Sh Se 845-2006640 
4 616 Timber Toy fs 845-2006640 
5 616 Timber Side/ 16 ft 845-2006640 

Concrete Keel 

6 726 Composite 3) ie 845-4862749 
7 726 Composite U7 aw 845-4862749 
8 726 Composite Ie. aie 845-4862749 
9 688 Composite/Timber 12. £6 845-4403511 
10 637 Composite/Timber 2eet t 845-2140554 
ital 637 Composite/Timber tt 845-2140554 


=o 





Table 5.2 lists the vessel radii, weight, distance from center 
of gravity to keel, KG, moments of inertia about vessel center 
ee gravity, To.g. and keel, I,, and side block transverse 


spacing, B. 


TABLE 5.2 
Key Vessel Parameters 


ieGblil dl 
SYSTEM Radius Weight KG 


(inches) (Kips) (inches) (K- in- gec2 ees eac2) ras’ 


1-5 198 16,396 193 831,257 2,411,000 144 
6-8 252 37,656 223 3,097,535 7,949,000 180 

9 198 13,624 193 691,852 2,007,000 138 
10-11 190 9,529 174 445,592 1,193,000 138 


5.2 Block Parameters!1 

The models developed in Section 4 require the evaluation 
of the block's vertical and horizontal spring stiffness 
constants, Ky, and K,. The spring stiffness is the amount of 
force required to produce one unit of displacement in the 
respective direction. In order to facilitate the use of 
Hooke's law in evaluating spring constants, the assumption of 
homogeneous isotropic behavior will be used for the block 
materials: wood and concrete. Both the side and keel piers 
must be analyzed for K, and kK;,.- 

5.2.1 Vertical Stiffness K,11 
Using Hooke's Law, the linear relationship between 


vertical applied force and deformation can be written as 
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where F = vertical force 
A = area under applied force 
E = modulus of elasticity 
h = change in height 
ig 


= original height. 
Hence, the spring constant for a given material is 


(561) 

For the eleven drydock configurations, two types of 
standard blocks (composite and timber) are used. A standard 
composite block which is composed of a softwood cap, a 
hardwood middle portion, and a concrete bottom can be modelled 
as a series of three springs (Figure 5.1). The resultant 
block spring stiffness is 

Ky block = [(1/Keap) + (1/Koak) + (1/Keon) 17+ (5.2) 
where 

Kcap = softwood cap spring constant 

Koak = hardwood middle spring constant 

Keon = concrete bottom spring constant. 

All three constants listed above are calculated from 
Equation 5.1. A standard timber block is composed of a softwood 
cap and a hardwood body (Figure 5.2). Its block spring 
constant is 

Ky plock = [(1/Keap) + (1/Keax) 17+: (5.3) 

The total vertical stiffness for a drydock configuration 
can be computed by multiplying ky pioch. times the number of 


blocks. The keel pier has a total vetical stiffness of 
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Kky = Ky block * Dk 
and for the side pier 

Ksy = Ky block * Ms 
whre ny and ng are the number of blocks in one row for their 
respective pier. 

The values of total vertical stiffness for both keel and 
side blocks of the eleven drydock configurations are found in 
Table 5.4. The following parameters are used in the 
formulation of Table 5.4: 

A = cross-sectional area of softwood cap 


Ecap =m KS. 


eee ois675 Kei 

ean = 2000 ksi 

Hoap = 4 inches 

Hoak = 29 inches for composite block 
= 56 inches for timber block 

Hoon = 27 inches. 


5.2.2 Horizontal Stiffness K,11 
In order to determine the horizontal stiffness, two types 
of deformation must be looked at. They are the block-cap 
displacements due to bending and shear deformations. Modelling 
aS a continuous cantilever beam subjected to a concentrated 
lateral force at the cap surface (Figure 5.3), the bending 


displacement due to the applied force P of a composite block is 


P(H,+H,+H3) P(E, -E.) (H5+H3) 


ee 
B SE, 1 3E, E51 
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i (E,I-E,1,) 43 
+ jae (5.4) 

where H, = height of concrete (27 inches) 

Hy = height of hardwood (29 inches) 

H3 = height of softwood cap (4 inches) 

E,; = modulus of elasticity of concrete 

E> = modulus of elasticity of hardwood 

E3 = modulus of elasticity of softwood cap 

I = moment of inertia of block's cross-section 

I3 = moment of inertia of cap's cross-section. 
For the bending displacement of a timber block, Equation 5.4 
still holds true with the alterations: 


Hy height of hardwood (27 inches) 


E,; = modulus of elasticity of hardwood. 
In shear, deformation can be determined by modelling the 


composite block as an element subjected to shear stress at the 


top (Figure 5.4). The shear displacement becomes 
We con I wood) 2 
S AE AE 
Z 
i BY eel, ots 
AE 2 


(5.5) 
where A = area of cap's cross section 
Vwood = Poisson's ratio for wood (0.30 is used) 
Vaon = Poisson's ratio for concrete (0.15 is used). 
By redefining H, and Ej, Equation 5.5 yields the shear 


displacement of a timber block. 
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The horizontal spring constant of any one block can be 
determined through the Hooke's Law relation 

ky = P/(dp + dg) (5.6) 
The total horizontal stiffnesses for the keel and side piers 


are 


Keh = Ky * Np 
and Key, = Ky, *¥ Ne .- 

The values for the total area of the rows of blocks for 
each of the eleven systems are listed in Table 5.3. Also listed 
are the corresponding total area of the cap blocks for a row of 
side or keel blocks. The values of total horizontal and 


vertical stiffnesses of the keel and side piers for the eleven 


adrydock configurations are listed in Table 5.4. 
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Total Keel and Side Pier Stiffness 


System 
il 


2 


10 


ii 


Kyy 
50089159 


SVOCI So 
25495109 
28493109 
50089159 
26250695 
9S256695 
Sle) SUS ENe. 2) 
21758374 
Hoo 7509 


6577809 


TABLE 5.4 


Ksv 


7545899. 
3903052. 
4292468. 
2220242. 
2220242. 
S2 05736 
34542004 
26345596 
8229775. 
6765289. 


5464271. 
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in Pounds Per Inch 


Kkh 
14018741 
14018741 
5161684. 


5161684. 


14018741. 


24330635 


24330635 


24330635 


3941650. 


8072572. 


S10 257 2. 


Ksh 
2612721 


1351407. 
1282216. 
663215.3 
SGo2 ions 
16104611 
10676090 
Sil427 30. 
Z2os32 9. 
ZS 9500. 


P7258 057. 
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FIGURE 5.1 
Standard Composite Block Stiffness Model 
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FIGURE 5.2 
Standard Timber Block Stiffness Model 


-58- 








FIGURE 5.3 
Standard Block in Bending 





FIGURE 5.4 
Standard Block in Shear 





6.0 VESSEL RESPONSE TO SEISMIC EXCITATION 

This portion of the thesis discusses methods of 
determining the seismic response of the vessel-drydock system 
mathematical models developed in Section 4. The numerical 
analysis schemes to be implemented are a modal analysis method 
(Section 6.1) which looks at the linearized equations of 
motion, and a fourth order Runge-Kutta numerical method 
(Section 6.2) which analyzes the non-linear aspect and also 
verifies the linear solution. Section 6.3 incorporates a 
response spectrum analysis to solve the mathematical models. 
Finally, Section 6.4 generates system response for the eleven 
vessel-drydock configurations using the quasi-static force 
method. 

6.1 Vessel Response Using Modal Analysis 

In order to use the modal analysis method described 
herein, only the linearized equations of motion mathematical 
models can be evaluated. These are Equations 4.2.8, 4.3.11 
and 4.4.11 which correspond to the one, two, and three degree 
of freedom models respectively. Since all three linearized 
equations are in the form 

[A](¥} + (Bl{y} = (E(t)}, 
the three degree of freedom model will be analyzed. The other 
two cases are just a reduced version of the third. 

To solve the linearized three degree of freedom 
equations of motion (4.4.11) in matrix notation and obtain 
system responses as a function of time, the problem is 
solved as an eigenvalue problem with a procedure described in 


reference [13]. As a first step, the equations are decoupled. 
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To accomplish this, first consider the free vibration problem 
UAC se (EH Oehe = (6.1.1) 
Assuming Equation 6.1.1 has a solution in the form 
(yy = (y} et, 
Equation 6.1.1 can be expressed as 


ATA] {Wy} + [BJ] {v} = 0. (6.1.2) 


where Xi complex eigenvalues 


{ V} modal vectors. 
Solving the matrix determinant problem 
| es], ae Es (6 9153) 

yields the eigenvalues. In the three degree of freedom case, 
[A] and [B]) are 6 x 6 matrices hence the solution of Equation 
6.1.3 will yield 6 eigenvalues. These eigenvalues will be 
comprised of three sets of complex conjugates. Now the modal 
vectors can be found by substituting the eigenvalues back into 
Equation 6.1.2 and solving. The modal vectors are then 
combined to form the modal matrix [1] where 

CV} = C€v}_q (Vo C3 (Vd4 CV}5 (V6). 
The equations of motion can then be decoupled by using the 


modal matrix, i.e. 


Cw? (A)CVv) 


[—A—~] 


and 


[¥)? ([BI(v) = [-B~] 
where [A -_]} and [“B~-_] are diagonal matrices. 
Once the modal matrix is found, the forced vibration 


problem is analyzed. Assume the following 


[W) {2) 
[v]7+ (y), 


{y} 


ana {Zz} 


-~6]= 





where z is a complex coordinate system. Equation 4.4.1 is now 
expressed as 
[A] (¥] {2} + (B] [¥] (2} = (E(t)). 
Using Cw)t, the above equation can now be in the form 
[~A~] {2} + [—B~] (2} = (¥)T{E}) = (N(t)}. 
This will yield six uncoupled equation of motion in the form 
Qiq Z4(t) + bis 24(t) = Nj(t), 1= 1,6 
eee 21(t) — Ys 2:(t) = (1/aqy) Ny(t) (6.1.4) 
where V3 = -bjj/ajiji- 


The solution of Equation 6.1.4 is 


Zi(t) = a Se vi (t-T) echt tl 6. (6.45) 
Once Equation 6.1.5 is solved, the {y} matrix is found by the 
relationship 
a= UZ} 
Recalling from Section 4.4.2 the terms of the three 
degree of freedom equations of motion {y} matrix is 
x 


Y 


CDe 


va 


SEI ° 
NY 
6 
Hence, the desired vessel-drydock system response can be found 
in the fourth, fifth and sixth terms of the {y} matrix, i.e. 


fey and 0. 
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6.2 Vessel Response: Fourth Order Runge-Kutta Numerical 
Method 


A fourth order Runge-Kutta numerical scheme is chosen to 
analyze the non-linear equations of motion due to its 
simplicity and easy use of computer programming. It also will 
be used to verify solutions found by other methods. The 
non-linear Equations 4.2.6; 4.3.6 and 4.3.8; 4.4.3. 4.4.4 and 
4.4.8 correspond to the one, two and three degree of 
freedom models respectively. Since all three systems of 
equations are similar in form, only the three degree of 
freedom non-linear model will be developed. 

The three degree of freedom non-linear equations of 
motion are 

M)] x + m4 36 eeeoagi x + C13 6 4+ Ki] X = “M1 Xg 

Mao y + Cop ¥ + kao y = -Mg2 Yg 
and 

m336 + m3, X - m3, yO+ 0330 + cg, & + k32° 

= =-M3] Xg + ™3] Yg O , 
The coefficient of the above equations can be found in Section 


4.4.2. These equations can be rearranged into 


a m C ; E ; k 
e+ 735 --2t, - 425 -—x«-« (6.1. 2a) 
igh eet iW Tk J 
. C ee z 
y = wer V eee ye any, (Ce. 25) 
20 no? J 
em... 3) | 33 4 _ 32 33 1 Smee os Bie 22 ae 
ee abi as 33. 22 ee 
. aes x (G2 1c) 
: ke 
wee 


ea 





Equations 6.2.1a, 


equivalent 


xXx=R= 
yY=S se 
S=T = 
mM 
3 ie ee 
le 


e; (t,x,R,T) 


f,(t,y,S) 


6.2.lb and 6.2.1c can be made into an 


first order system by the following substitutions 


Gi ctgi5y Rs S,T) 


T = e, (t,x,R,T) 


S = f (ep y 7S) 


Z 
: m 
‘tag ~— 

irs 


R= go (t, diye oS ° 


Runge-Kutta formulas [8] for this system are 


= hR, 
= hs, 


= hT, 


1 
= h(R, +2¢ 

1 
= h(S, +2 


il 
= h(Tr ao 


= hA(S, + 


= h(T, + 


a9) 


{12 


£43 
11) “21 


£12) *22 


£13) 423 


1 
Doe sai 
1 

792) %32 


1 
2493) £33 


h{eo(t,Xn,Rny-Tn) } 

h{(fo(t,Yn/Spn) } 

ican cing, Yn/Rn/Sn/Tp) } 

h{eg(t,%, + 7ha1,Rn + st11/Tn + 5013)) 
H{f2(t,Yn + 2a208r + 2)? 

h{gg(t,0n + ee reas, 

Rn + fair: Sn + Lt, Tn + =£13)} 
h{e5 (t,x, + Fko1,Rpn + s¢21/Tn 43123) } 
h{fo(t,ypn + 2 221Sr ~ 199) ) 
h(gg(t,°n + Tha31Yn + sf22/Rn 5 sti 


1 24 
Sy + 3*22,Tn + 2%23)) 


-64= 





+ ata i €33)) 


an = h(Rp + £31) 41 = h{eg(t,xXp + Raa Rn 
Rao = h(Sp, + £32) ip = h{fo(t,Yn + Rao, Sy + Saad 
Ra = h(Tpy + £33) 43 = h{gg(t,0p, + R53, Yn + to Rn oF ont 


+ *33)) 


Sn + “32, Tn 

and lastly, 

Xn+1 = Xn t z( Py + 25, + 2h31 + ay) 

atl 7 *Yn * =( hy + 2&5 + 2732 + Rao) 

On+1 = Sn Ft =(hy3 + 2ko, + 2233 + 243) 

Rn+1 = Rn + =(£14 + 245, + 243, + £41) 

Sn+1 = Sn + =(£12 + 2459 + 2439 + £42) 

Tht = Tn + 2(413 + 2423 + 2233 + 443) 
where h is the time step. Hence, the desired vessel-drydock 
system response can be found and system failure analyzed by 
the iteration of the Runge-Kutta scheme. 

6.3 Vessel Response: Response Spectrum Analysis 

The degrees of freedom mathematical models developed in 
previous sections produce a time history of the seismic 
response of a drydocked vessel. These models use numerous 
computer iterations to find their time history in order to 
evaluate the system's maximum seismic response. Another way 
to determine the maximum seismic response of a drydocked 
vessel is the response spectrum analysis method. Simplicity 
and fast results are characteristic of this spectrum method. 
However, only maximum response is produced and no time history 
is generated. 


The response spectra are graphs of the maximum seismic 


response of single degree of freedom systems over a range of 
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system natural frequencies. These graphs are generated for 
particular earthquake and can contain acceleration, velocity, 
displacement, or combination of all three responses for the 
given system. Simplified response spectra are plotted for 
earthquake design analysis. The El Centro earthquake response 
spectra graph, shown in Figure 6.1, will be used to confirm 
the results of the mathematical models in this thesis. The 
simplified spectrum portion of this graph is described in 
reference [2] and updated in reference [7]. 

The simplified El Centro earthquake response spectra used 
described the maximum relative displacement of the system 
shown in the insert of Figure 6.1. The linear equation of 
motion for this simple system is 

mu + ci + ku = -m Tae f, (t) 


where u(t) = relative displacement of the mass with 
respect to the support 


and Ysa maximum support acceleration. 


The natural frequency of the system is given by 
we = k/m 
f= W/a2T , 
Using the response spectra graph and the particular frequency, 
f, of the system of interest, the maximum value of the 
relative displacement, u,, is determined. In equation forn, 
the relative displacements of a single degree of freedom 
system with 5% critical damping subjected to the El Centro 
earthquake ground motions are: 
for w < 2.24 rad/sec, 


Umax = 1-4 (Ysa)max = 11-62 inches, (Gece) 
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home 2.24 < | < 12.74, 


Umax = 1-9 (Ysa) max/ 26.03/w inches, iGo 2) 

POm > 12.74; 

Dee ean nG(yeayere io = 331.53/42 inches. (6.3.3) 
These formulas generated the simplified response spectrum 
shown in Figure 6.2. 

In order to use Figure 6.2 to verify the one, two and 
three degree of freedom vessel-drydock models, some 
manipulation must be performed. This is because the linear 
equations of motion of the models do not exactly match those 
of the system pictured in Figure 6.1. Examining the one 
degree of freedom linear equation of motion (Equation 4.2.7); 

T,6 + 2EI_W_o +[ (B2/2) Key - WKG]® = -M KG xg 
or 

6 + 2—w ,8 + ([(B2/2)Key - W KG]/I,)® = -(M KG/Ik) ere 
The natural frequency for the above equation is defined as 

Wn? = ((B2/2) Key - W KG)/I,-. 

Using Figure 6.2, this frequency gives a corresponding Uma,- 
Hence, the maximum rotational response of the system, , is 
then determined by the relation 

Ymax = M KG/Iy Umax: 
Table 6.1 lists the maximum rotational response for the eleven 


configurations found in Section 5. 
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TABLE 6.1 


One Degree of Freedom Vessel-Drydock System 
Maximum Response Using Response Spectrum Analysis 


System =n Smax 
(rad/sec) (rads.) 
ih Deo 0.01594 
2 3933 0.02619 
3 4.141 0.02149 
4 2.870 0.03100 
5 2.870 0.03100 
6 10.254 0.00698 
7 S.32/ 0.00860 
8 Wee loo 0.00987 
9 6.147 0.01445 
10 TeZo3 0.01299 
11 6.498 0.01450 


The response spectrum analysis to the two and three 
degree of freedom systems is applied by means of the modal 
method described in reference [2]. Concisely, the modal 
method requires that the system natural modes be determined. 
Then, the response of the system to a known forcing input can 
be developed by mode superposition. Treating each natural 
mode as a single degree of freedom system with its own natural 
frequency, this method is good for simplifying seismic 
response analysis. Since each natural mode acts as a separate 
system, the maximum response may be determined. Then, a 
conservative estimate for the maximum response of the original 
system can be made by adding up the maximum response of each 
mode. 

The method of adding maximum modal responses gives an 
upper bound for the maximum system response. Since there is a 
known response spectrum, the amount in which each natural mode 


contributes to the maximum response for a given input can be 
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determined. This is known as the participation factor of each 
given mode. A full description and derivation of partici- 
pation factors including implementation techniques can be 
found in reference [2], and will not be presented here. The 
participation factor for each natural mode in the two and 
three degree of freedom case is calculated to confirm computer 
generated results of maximum responses in Appendix 3. 

6.4 Vessel Response: Quasi-Static Force Analysis 

As described in Section 4.1, the quasi-static force 
method replaces the earthquake motions by a force 
corresponding to the vessel mass times 0.29. This analysis 
only allows rotation about the vessel's keel hence, the quasi- 
static force analysis is itself a one degree of freedom model 


with system response, in radians, determined by Equation 


iy .4, i.e. 
8 = M/L* Key 
where Ms = 448 * A * KG 
A = displacement in tons 
L = B/2. 


The results of this method are listed in Table 6.2. 
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TABLE 6.2 


System Response using Quasi-Static Force Analysis 


SYSTEM 0, IN RADIANS 
1 0.00809 
2 0.01564 
3 0.01422 
4 0.02749 
5 0.02749 
6 0.00199 
7 0.00300 
8 0.04419 
9 0.00670 

10 0.00051 
11 0.00637 
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FIGURE 6.1 
Actual El Centro Earthquake Response Spectrum 
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FIGURE 6.2 
Idealized El Centro Earthquake Response Spectrum, 5% Damping 
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7.0 COMPUTER SOLUTIONS 


In this section of the thesis, the computer programs 
needed to implement the mathematical models and the associated 
seismic response are developed. Due to the numerous 
iterations that is required for the modal analysis and 
Runge-Kutta methods to effectively determine vessel-drydock 
response, all computer work is performed at the MIT Joint 
Computer Facility (JCF). The required computer programs are 
written in Fortran 77. 

In order to check the computer generated results, the 
vessel-drydock mathematical models are subjected toa 
sinusoidal earthquake input. The advantage being that closed 
form solutions of responses to sinusoidal inputs can be 
calculated. Thus, closed form response can be compared to the 
computer generated responses. 

Upon obtaining the correct results for a sinusoidal 
earthquake, the El Centro earthquake's acceleration time 
history is applied to the vessel-drydock system. The maximum 
value of these computer generated results is compared to the 
maximum response for the system predicted by the response 


spectrum method found in Section 6.3. 


7.1 Modal Analysis Solution 
7.1.1 Computer Program Development 


As described in Section 6.3, the program approaches the 
problem by first assembling matrices A and B. The program 
then uses the International Mathematical and Statistical 


Library (IMSL) subroutine EIGZF resident on JCF to perform the 
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system free vibration eigenvalue problem on matrices A and B. 
The output of subroutine EIG2ZF is the modal matrix [WwW] of the 
vessel-drydock system. Reference [4] contains further 
information about subroutine EIG2F. 

The equations of motion can now be de-coupled by using 
the modal matrix. De-coupling the equations of motion leaves 
equations of the form 

ai; z3(t) + bss 2y(t) =NG(t) , i =i, m 
where m is equal to twice the number of the particular degree 


of freedom model, i.e. for the one degree of freedom ~ m = 2. 


The solution to the above equation is 


te 
l “bil/aii (t-T) | 
Z4 (TC) ge: ae € Nj (T) dcr 1 = 1,mn, 
AGL 
O 


where aj; and bj; represent the diagonal elements of the 


Giagonalized A and B matrices, respectively. The integral 


t 
ee) 
== N. (tT) at 


O 


can be evaluated using Simpson's rule, which is 


PAs(t) = [(Az(t - 2At) + Ny(t - 2At))e7 (Pii/4ii) 2dr 


+ 4N3(t - Atye7(Pii/aiiA+ Ns (t)] 


where At = data time increments 
and A(t - 2AtT) = LA(t) evaluated two time increments 
previously. 


The value z;(t) is evaluated as 


Ge eee / as) (0/3) (As (t)), 1 = 1,m. 
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To obtain the values of the system response for a 

particular time, the relationship 

ee ee 23} 
where {y} = system response vector is used, using the column 
vector [z] for the time at which results are desired. The 
values of block forces caused by the system response are 
calculated every 0.02 seconds for the first 30 seconds of the 
earthquake. Only the first 30 seconds are needed due to the 
nature of the El Centro earthquake the seismic acceleration 
after 30 seconds are small when compared to the first 30 
seconds. The resultant block forces and vessel motions are 
then checked against failure criteria established in Section 
3. The amplitude of the seismic acceleration is reduced until 
no block failure occurs. Hence, the maximum earthquake 
acceleration, in g's, that the vessel-drydock system can 
withstand without failure can be found. 

The three degree of freedom vessel response computer 
program is listed in Appendix 1. Since all three linearized 
models found in Equations 4.2.8, 4.3.11 and 4.4.11 are in the 
form 

PA PE) = {ECe) }; 
the three degree of freedom computer program is modified to 
give the one and two degree model results. 
7.1.2 Computer Program Input 

The computer program listed in Appendix 1 evaluates the 
horizontal and vertical translations and rotation of the three 
degree of freedom vessel-drydock system in increments of 0.02 


seconds. The program requires the following input: vessel 


Giz 





weight, WEIGHT, in kips, keel to center of gravity height 
(KG), H, in inches, vessel mass moment of inertia about the 
keel, Ip, in kip-in-sec*, time increments of ground 
acceleration data, DTAU, in seconds, side and keel piers 
vertical and horizontal stiffnesses, Key, Kpy, Kspy, and Ky, 
in kips/inch, gravitational constant, GRAVITY, in inch/sec?, 
side and keel pier dimensions of width and height, BASE SIDE, 
BASE KEEL, HTSIDE and HTKEEL, in inches, the coefficient of 
friction between the block-block and ship-block interfaces, Ul 
and U2, horizontal distance between the center of vessel 
contact with port and starboard side piers, BR, in inches, 
vertical/horizontal ground acceleration ratio, AMP, 
proportional limits of side and keel block materials, PLSIDE 
and PLKEEL, in kips/in?, total cross-sectional area of side 
and keel piers, SIDE AREA and KEEL AREA, in in *, percent 
critical damping, ZETA, vessel hull number, HULL, and finally 
the vessel-drydock configuration number, NSYS. The following 


is where these inputs can be found: 


Hable 5.1 NSYS and HULL 
Table 5.2 WEIGHT, H, I,, and BR 
Table 5.3 BASESIDE, BASEKEEL, HTSIDE, SIDEAREA 


and KEELAREA 
Table 5.4 Key, Ksy, Karp and Kep- 
The rest of the inputs are assumed to be 
GRAVITY = 384 in/sec? 


Ul = U2 = 0.5 


PLSIDE = PLKEEL = 0.7 kips/in@ 


and AMP = 0.5. 
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The value of 0.5 for the vertical/horizontal ground acceler- 
ation ratio (AMP) parameter is used because no vertical 
component of the El Centro earthquake acceleration was 
available at the time of this thesis. 

Damping coefficients equal to 5% of critical damping 
(i.e. ZETA = 0.05) are assumed for the entire system. This 
analysis for damping coefficients can be found in Appendix 3 
and is implemented directly into the computer program. 

Certain idealized system failure criteria from Section 3 
need to be modified to accurately represent the vessel-drydock 
system as it actually exists. As stated in Section 3.1, block 
crushing occurs when the stress on the block exceeds the 
block's proportional limit. The vessel rests upon soft wood 
caps placed on top of the keel and side piers when in the 
drydock cradle. Crushing of the cap is not considered a 
failure due to its small size. For the blocking system 
considered, the generated stresses are transferred through the 
cap to the top drydock block. 

Also, a side pier cannot tip inboard because the vessel 
hull is not physically attached to the pier and cannot pull 
the pier beyond it's upright position. The softwood cap on 
top of the side pier causes the vessel force to be applied 12 
inches from the inboard edge. As shown in Figure 7.3, this 
causes the static force vector to fall outside of the middle 
one-third of the pier base, as discussed in Section 3. Since 
the pier will only tip in the outboard direction, the 


resultant vessel force must fall in the outboard one-third of 
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the pier base. Thus, this failure criteria (from Figure 7.3) 


for the side pier is 


HORIZONTAL FORCE . 27 Otc tie 2,0 
VERTICAL FORCE Pier Height 


7.1.3 Computer Program Testing 


In order to validate the modal anlaysis method, the one 
degree of freedom mathematical model is checked with a set of 
sinusoidal ground accelerations for input. The input ground 
motion was 

E(t) = Xg = 323.95 sin(3.92t). (Gaikes 
This sample function is selected because its magnitude and 
frequency closely match that of the maximum acceleration 
portion of the El Centro earthquake. Also, vessel-drydock 
system configuration #1 is chosen for this initial analysis. 
The results of the one degree of freedom case will be compared 
to the closed form solution. 

The one degree of freedom linear equation of motion 
(Equation 4.2.7) is of the form 

e+ 2,6 + one a EC : 
The closed form solution of this equation is 


te 
MK ~Ew, (e-T) ; 
a(t) =- ECT) e sinw,(t-T) at 
Lay 
| eae 6 | A 


al 








= ISS eaecimu:to— Bit) —cosm. t] OTS Als ) 
Tw d d 
where £ E wit 
e 
_ d 
A(t) f E (tT) x: COSW4T AT 
O e 


252) — 





Ew T 


B(t) = E sin w a 
= ae 
= 





The results of a computer run with vessel-drydock system 
#1 and sinusoidal excitation (Equation 7.1.1) indicate that 
the maximum rotation angle 6 is 0.038243 radians at a time of 
2.02 seconds. To confirm this value, Equation 7.1.2 has to be 
evaluated 


—o 


9(2.02)=-U RS [a(2.02) sinwg(202) - B(202) coswg(2.02) ] 


1, Wg 
22.02 Ew tT 
= ae 
where A(2.02) - | B26..95 Sims. o2 Tt) 2. 02Eu COSWqg T AT 
0 e n 
202 
EwWnt 
° e ‘ 
and b(2.02) = Eye oS Sims, 927) “Zoran tat 8 
e 
0 


From Section 6.3, 


Wy = [((B2/2)Key - W KG)/Iy) 1/2 


5.580 rad/sec 


and Wg = Wy (l = p2yl/2 
= 5.573 rad/sec. 
With these values of natural and damped natural frequencies, 


the value of A(2.02) can be evaluated using the relationship 


te 





if Qat EGC osctudt = eat BORE eee ee eons ele 
0 2{[a°+(b-c) “J 
Te 
ft at {a sin(btc) t- (btc) cos (btc) t] 
: — A a 
2{[a +(btc) 7] 0 
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Substituting the values 


t = 2.02 
a= Wy = .279 
b = 3.92 

ama, C= UG = 5.573 


The value of A(2.02) becomes 
323.95(e°7536)-1 & - 18307 = -153.1629 


eamilarly, for B(2.02), 


[ Bees eG ot aene.aat ee eee) racos (boc) t] 
2[a +(b-c) 7] 
0 t 
pare [(b+c)sin(b+c)ttacos (b+c) t] 
2[a2+(b+c) 2] 


with the same values of a, b and c as before. 
Evaluating, B(2.02) becomes 
323.95 (e°7636)-1l * -.2191 = -40.3973 


Now, 


* 
1, 42-69797193-0 | _153.1629*sin(2.02*5.573) 


oo. 02) = 
2.54 2411000.*5.573 ~-(-40.3973%*cos (2.02*5.573) ] 





= .038247 radians. 

Comparing the computer generated solution (0.038243 
radians) with the closed form solution (.038247 radians), it 
can be seen that the computer program is calculating the one 
degree of freedom rotation formula correctly. 

One other simple method to check to see if the rotation 
displacement time history is correct is to look at the 
waveform. At steady state condition, 9®(t) will be a 
Sinusoidal with the same frequency as the sinusoidal ground 


acceleration. An amplitude different and a phase shift will 
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be present. The steady state solution [12] of a equation of 
motion in the form 


mx + cx + kx = F sinvt 


a —— ean = nit sin [.t-tan™ a 
kV (1-r2)24+(2¢r)2 


l-r 
W/W, 


(k/m) 1/2 = Natural Frequency 


where a 


On 


ll 


1/2 c/v km = percent critical damping. 


In this case, 


W, = 5.580 rad/sec 
r= 3.92/5.580 = .703 
0S an 

__ |x 
IF| = M RG —2—= 1051012.9 lb.in 
204 
k = 75071464 1b in. 
so now, 
8 steady state = -027402 sin(3.92t - .1381) radians (7.1.3) 


Figure 7.1 represents the sinusoidal ground acceleration 
(Equation 7.1.1) and Figure 7.2 which represents the rotation, 
0(t), of the vessel-drydock system as generated by the modal 
analysis computer program. Note the same frequency between 
Figures 7.1 and 7.2. The steady state portion of Figure 7.2 
is matched to Equation 7.1.3. In conclusion the model 
analysis method correctly predicts the one degree of freedom 


mathematical model. 
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7.1.4 Computer Program Results 


The modal analysis method can now be used with the El 
Centro earthquake data to predict the system response. Eleven 
vessel-drydock configurations (Section 5.0) are analyzed for 
maximum ground acceleration that the system can withstand 
without block failure. Also, the first block failure modes 
exhibited by the systems are found. 

For the eleven vessel-drydock systems, the linearized 
one, two and three degree of motion mathematical model 
responses are calculated and listed in Table 7.3, 7.4, and 
7.5, respectively. These tables can be found on pages 4] 
thru93. The program further calculates that vessel liftoff 
(Section 3.4) will be the first failure mode to occur when 
ground accelerations are greater than the maximum permissible 
ground acceleration listed in Table 7.3, 7.4, and 7.5. The 
maximum permissible ground acceleration is based on the 
reduction of the El Centro acceleration amplitude until no 
failure occur. However, the waveform remains unchanged. 

In order to get a visual picture of the various 
responses, vessel-drydock system #1 is selected to be plotted. 
Figures 7.4; 7.5a and b and 7.6a,b, and c are the appropriate 
responses for the one, two and three degree of freedom models. 
Only the first thirty seconds of the earthquake are required 
to give the maximum system response as seen by the plots. 

Note that the relative horizontal displacement (Figure 7.6a) 
and rotation (7.6b) time histories of the three degree of 
freedom model are indentical to the ones found in the two 


degree of freedom model. This is due to the uncoupled 
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vertical equation of motin found in the linearized three 
degree of freedom case. 

As a check for the computer generated results, the system 
maximum responses due to full magnitude El Centro earthquake 
accelerations are compared to the system maximum responses 
determined by response spectrum analysis in Table 7.1 and 7.2. 
These tables are generated by combining the computer results 
of various models and Tables 6.1 (Section 6.3) and A.3.3 
(Appendix 3). Since the response spectrum method is an 
approximate one (Section 6.3), the maximum 20% error between 
its results and the computer results in the one degree of 
freedom case is reasonable. Thus, the one degree program is 
generating accurate results. 

In the two degree of freedom case, the response spectrum 
method predicts the worst case if the system's two natural 
modes reach a maximum at the same time, and add together. 

This happened in the rotational case, but only twice in the 
relative displacement case. In the majority of the systems, 
the computer generated rotational response actually exceeded 
the rotational response predicted by modal analysis. The 
maximum 35% difference between the two results can be 
accounted for because of the approximate, simplified response 
spectrum mentioned earlier. The three degree of freedom case 
has the same order of error as found in the two degree case as 
shown in Table 7.2. Hence, all mathematical models perform 


satisfactory using the modal analysis computer program. 
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TABLE 7.1 


One Degree of Freedom Equation of Motion Response Comparison 
Modal Analysis Computer Results of Linear System 


1 o | 
System radians 
ak 0.01474 
2 0.01618 
3 0.01436 
4 0.02966 
5 0.02966 
6 0.00832 
7 020115 
8 0.01712 
9 ©. 01512 
10 0.01714 
ay) 


Response Spectrum Method Results of Linear System 


| 9 | 
System radians 
le 0.01594 
Z 0.02619 
3 0.02149 
4 0.032 
5 0.031 
6 0.00698 
7 0.0068 
8 0.00987 
9 0.01445 
10 0.01299 
les 0.0145 
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TABLE 7.2 


Three Degree of Freedom Equation of Motion Response Comparison 
Modal Analysis Computer Results of Linear System 


| ay | eC ee 
System (inches) (inches) (radians) 
al OZ 7971 07073174 0.01458 
2 0.22661 0.10948 0.01636 
3 0.53288 Ores 395 oO, OLS46 
4 0.58122 Oes352Z O03 1530 
5 0.24721 Oni] OS 0.02992 
6 0.36179 O20612 3 0.00814 
7 0.39118 0.06870 0.00879 
8 0.39889 0.07140 O701154 
9 0.61065 0.11749 OZ O01587 
LO 0.57808 OL Alea ule ee 8, ORO 547 
11 OO. Gis @ O. 11676 0.01734 


Response Spectrum Method Results of Linear System 


| x | | sae | been! 
System (inches) (inches) (radians) 
ag ORs 3275 0.10859 0.01646 
2 0.27141 O212225 0.02348 
c' 0.65142 0.19089 0.02350 
4 0.53788 0.21491 G02 5.510 
5 0.21468 0.12980 0.03202 
6 0245271 e..08025 0.00683 
7 0.46833 0.09708 0.00866 
8 0.46597 Gee 762 0.01009 
9 0.64366 O22 5562 0701552 
10 0.65864 0.13662 0.01342 
a1 0.67467 0.14955 201552 
NOTE: Two Degree of Freedom System Maximum Response 


correspond to x and 0. 
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7.2 Fourth Order Runge-Kutta Analaysis Solution 
7.2.1 Computer Program Development 


As described in Section 6.2, the Runge-Kutta scheme uses 
a set of first order differential equations to describe the 
second order vessel-drydock system. A simple example of this 
procedure is done on the following equation 

mx + cx + kx = f(t) + g(t)x 


let r = x and rearranging yields two first order equations, 


25 SS Ng 


r = -=r - = x semeds (ate yt mele x 
Now, the first order system can be implemented into a computer 
program. 

Briefly, the fourth order Runge-Kutta numerical scheme is 
one that makes use of predicted velocity and acceleration at a 
given time to compute the displacement and velocity for the 
next time increment. This can be done by the basic 
relationships 

Xn+]1 = Xn + Xn At 
and nt] = Xn + Xp At 
However, x and X are not found by differentiation but 

through the manipulations of the second order equation into a 
set of first order equations. This can be done by procedure 
shown in the previous paragraph. 

The first step in the computer program development is to 
place the set of first order differential equations into the 


Runge-Kutta formulas shown in Section 6.2. The evaluation of 


these formulas will predict the velocity and acceleration 
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that will be used in the next time increment. Now, the system 
response, both displacements and associated velocities, can be 
calculated using the previous information. This process is 
marched through time steps in our to produce a time history of 
the system response. 

Now, the values of block forces caused by the system 
response are calculated every 0.01 seconds for the first 30 
seconds of the earthquake. The first 30 seconds are only 
required to analyze for maximum seismic response. The 
resultant block forces and vessel motions are then checked 
against failure criteria found in Section 3. The amplitude of 
the seismic acceleration is reduced until no block failure 
occurs. Like in the modal analysis method (Section 7.1), the 
maximum earthquake acceleration, in g's, that the 
vessel-drydock system can withstand without failure can be 
found. 

The Runge-Kutta analysis gives the ability to analyze 
both the linear and non-linear equations of motion for the 
one, two and three degree of motion mathematical models. The 
non-linear three degree of freedom vessel response computer 
program is listed in Appendix 2. The program is developed 
around Equations 4.4.5, 4.4.6 and 4.4.8 and those found in 
Section 6.2. Since all linear and non-linear models developed 
in this thesis can be arranged into the Runge-Kutta formulas, 
the non-linear three degree of freedom comptuer program is 


modified to give desired model results. 
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7.2.2 Computer Program Input 


The computer program listed in Appendix 2 evaluates the 
horizontal and vertical translations and rotation of the 
non-linear three degree of freedom vessel-drydock system in 
increments of 0.01 seconds. The program requires the 
identical input found in Section 7.1.2 for the modal analysis 
method. Hence, the program inputs needs not be listed again. 
Like in the previous case, Appendix 3 provides the damping 
coefficients to be implemented into the computer program. 

7.2.3 Computer Program Testing 

In order to validate the fourth order Runge-Kutta 
numerical scheme, the linear one degree of freedom 
mathematical model is checked with a set of sinusoidal ground 
accelerations for input. The same input ground motion that 
is used to validate the modal analysis is implemented in this 
case, l.e. 

E(1t) = Xg = 329.95 sin(3.92t). (Tae) 
This sample function is selected because its magnitude and 
frequency closely match that of the maximum acceleration 
portion of the El Centro earthquake. Also, vessel-drydock 
system configuration #1 is chosen for this initial analysis. 
The results of the linearized one degree of freedom case will 
be compared to the modal analysis and closed form solutions. 

The results of a computer run with vessel-drydock system 
#1 and sinusoidal excitation (Equation 7.2.1) indicate that 
the maximum rotation angle 6 is 0.038236 radians at a time of 
2.02 seconds. Comparing this the computer generated solution 


with the modal analysis solution (0.038243 radians) and the 
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closed form solution (0.038247 radians) found in Section 
7.1.3, it can be seen that the computer program is calculating 
the linearized one degree of freedom rotation formula 
correctly. Also, Figure 7.7 which represents the rotation, ® 
(t), of the vessel-drydock system as generated by the 
Runge-Kutta computer program is essentially identical to 
Figure 7.2 generated by the modal analysis computer progran. 

In conclusion, the fourth order Runge-Kutta method 
correctly predicts the linearized one degree of freedom 
mathematical model. Now, the computer program will be used 
for the linear and non-linear one, two and three degree of 
freedom systems. 

7.2.4 Computer Program Results 

The fourth order Runge-Kutta analysis can be implemented 
with the El Centro earthquake data to predict the system 
response. Eleven vessel-drydock configurations (Section 5.0) 
are analyzed for maximum ground acceleration that the ystem 
can withstand without block failure. Also, the first block 
failure modes exhibited by the systems are found. 

For the eleven vessel-drydock systems, the linear and 
non-linear one, two and three degree of motion mathematical 
model responses are calculated and listed in Table 7.3, 7.4 
and 7.5, respectively. These tables can be found on pages (I 
thru 93. The program further calculates that vessel liftoff 
(Section 3.4) will be the first failure mode to occur when 
ground accelerations are greater than the maximum permissible 


ground acceleration listed in Table 7.3, 7.4 and 7.5. 
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In order to get a visual picture of the various linear 
and non-linear responses, vessel-drydock system #1 is selected 
to be plotted. Figures 7.8, 7.9 a and b, 7.10 a, b and c are 
the appropriate response for the linear one, two and three 
degree of freedom models and Figures 7.11, 7.12a and b and 
7.13 a, b and c are for the non-linear models. Note that 
these figures are almost identical to Figures 7.3, 7.4a and b 
and 7.5 a, b and c generated by the modal analysis method 
(Section 7.1.5). This reinforces the equivalency of the 
results found by the two numerical shemes. This can also be 
seen by comparing the system maximum responses of the 
Runge-Kutta method to that of the modal analysis method. Both 
are listed in Table 7.3, 7.4 and 7.5. The maximum error for 
any given case is 5%. Hence, all mathematical models perform 
satisfactory using the fourth order Runge-Kutta computer 


program. 
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TABLE 7.3 


Modal Analysis Method Results of Linear System 


© max Time Maximum Acceleration 

System (radians) (sec) (g's) 
1 0.003390 6.64 0.0759 
2 0.003883 9.28 0.0792 
3 0.006033 9222 0.1368 
4 0.006822 8.68 0.0759 
5 0.004513 8.68 0.0462 
6 0.001999 536 0.0792 
7 0.002422 5.84 0.0924 
8 0.002760 8.98 0.0792 
9 6.005138 7.54 0.0990 
10 0.004537 8.98 0.0990 
11 0.004972 7.92 0.0957 


4th Order Runge-Kutta Numerical Method Results of Linear System 


Omax Time Maximum Acceleration 
System (radians) (sec) (g's) 
aL. 0.003390 6.64 0.0759 
2 0.003882 O27 0.0792 
3 0.006032 9.22 0.1386 
4 0.006822 8.68 0.0759 
5 0.004153 8.68 0.0462 
6 0.002004 5255 0.0792 
7 0.002426 SoS 0.0924 
8 0.002763 8.97 0.0792 
9 0.005141 Tee 0.0990 
10 0.004543 8.98 0.0990 
11 0.004972 7.92 0.0957 
4th Order Runge-Kutta Numerical Method Results of 
Non-Linear System 
Omax Time Maximum Acceleration 

System (radians) (sec) (q's) 
1 0.003391 6.64 0.0759 
2 0.003883 9.27 0.0792 
3 0.006032 9.23 Os S6 
4 0.006825 8.68 0.0759 
5, 0.004154 8.68 0.0462 
6 0.002004 53 0.0792 
7 0.002427 Seco 0.0924 
8 0.002764 8.97 0.0792 
9 0.005139 7.54 0.0990 
10 0.004546 8.98 0.0990 
11 0.004970 7.91 0.0957 
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TABLE 7.4 
Two Degree of Freedom Equation of Motion Response 
Modal Analysis Method Results of Linear System 
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TABLE 7.5 


Three Degree of Freedom Equation of Motion Response 
Modal Analysis Method Results of Linear System 


Maximum 
*max Time Y max Time max Time Acceleration 
System inches sec. inches sec. radians sec. 's 
1 0-067 132 8.38 0.017698 5.68 0.003501 9.02 0.0792 
2 0.061211 Dip wo OF02518a 8.06 0.003765 oO 26700), 0759 
3 0.207826 eZ 0.060041 DO 0.006033 O22 Ge Or 267 
4 0.122056 Deez 0.038540 5.64 0.006574 8.66 0.0693 
S O-032158 D3 70 O-015215 Sy oie) 0.003891 8.72 0.0429 
6 0.090448 5.78 0.019865 DOL 0.002037 5. SOMO Os > 
7 0.105619 5.46 0.018550 6.68 0.002375 59.86 0.0891 
8 0.091745 5.98 0.016424 D536 0.002654 7.62 020759 
9 0.201516 6.06 0.038773 Deol’, 0.005238 7.60 0.1089 
10 0.161863 6.04 OF 02356 51h Sirs 0.004332 7.84 0.0924 
11 0.179020 6.08 0.033861 5.40 0.005030 7 540.0957 
4th Order Runge-Kutta Numerical Method Results of Linear System 
c y 8 Maximum 
max Time max Time max Time Acceleration 

System inches sec. inches sec. radians sec. 'S 
1 0067373 8.38 0.017858 Soy 0.003504 07 O72 me 0 7.92 

Z 0.062228 5.39 0.025102 8.06 0.003761 Sf CBO SYS) 

5 0.207015 43 0.059736 Do 0.006304 9.25 0.1287 

4 0.121529 Dass O7038651 D205 OF006577 8.67 0.0693 

5 0.032442 DOL 0.014985 5.39 0.003891 8.79 0.0429 

6 0.094623 5278 0.019463 Syeei 0.002130 oye ie Os Ohopsile: 

7 0.106292 D4 / 0.018152 6.69 0.002378 D.60/ 0.0891 

8 0.091996 5.99 0.016558 rae: 0.002656 dace 0.0759 

9 0.200428 6.06 0.039205 5.40 0.005233 7.60 Q.1089 
10 0.161009 6.04 0.033887 5.39 0.004338 7.84 0.0924 
11 0.180863 6.08 0.034608 Dy 40 0.005033 55 OF 0957 


4th Order Runge-Kutta Numerical Method Results of Non-linear System 


Maximum 
* max Time ymax Time © max Time Acceleration 
System (inches) (sec.) (inches) (sec.) (radians) (sec.) (g's) 
1 OF 067253 8.38 0.017858 7, 0.003507 9.02 0.0792 
2 0.062050 5.39 0.025102 8.06 0.003763 97.27 © 020759 
3 0.206406 5. 7 0.059736 5.64 0.006037 93250 YOU 1267 
4 0.121484 aay ee 0.038651 26D 0.006580 8.67 0.0693 
5 0.032443 Syoneydl 0.014985 5.39 0.003891 8.79 0.0429 
6 0.094523 Sao 0.019463 Di 7 0.002130 Dao/ UnOSS6 
ii 0.106348 Say) 0.018152 6.69 0.002379 5.87 0.0891 
8 0.091968 5.99 0.016558 Dio 0.002656 7262. 0207 59 
9 0.200131 6.06 0.039205 5.40 0.005234 7.61 0.1089 
10 0.160892 6.04 0.033887 5.39 07004337 Fai UR Oi 
11 0.180674 6.08 0.034608 2.40 0.005032 7255. O70957 
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7.3 Discussion of the Models and Associated Results 

Both the modal analysis and the Runge-Kutta numerical 
methods give satisfactory results in determining the system 
Iaximum response in the various models and configurations. 
The maximum error in the comparison of the two methods for 
any given run is 5%. This is excellent when considering that 
each method has completely different basic concepts behind 
them. The modal aniysis method runs faster to obtain the 
result due to less computer operations. But, the Runge-Kutta 
method is much easier to program into the computer and is more 
flexible. In all, both methods can be used to predict the 
seismic response of the linearized one, two and three degree 
of freedom models. 

In the comparison of the one, two and three degree of 
freedom responses, certain conclusions can be drawn. The one 
degree of freedom model need only be used if system maximum 
rotational response is the main concern to drydock system 
design. This model gives results sufficiently close to the 
rotation generated by the other two models. Since vessel 
liftoff of the side block is the first block failure to be 
exhibited by the vessel-drydock system, the one degree of 
freedom model adequately predicts the maximum earthquake 
amplitude the system can withstand. Most of the vertical 
displacement at the side pier is due to rotation and not of 
relative vertical displacement as found in the three degree 
model. Hence, the one degree of freedom system response 
should be used to determine vessel roation and maximum 


earthquake amplitude. 
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The three degree of freedom system results should be used 
for a parameter study on block failure modes and for block 
design. This is because only the three degree model produces 
all required displacements for estimation of generated forces 
in both keel and side blocks. However, when the relative 
vertical translations generated are small, the two degree 
model can be used in block design,. The one degree model is 
truely inadequate for block design since crushing and vessel 
liftoff of the side blocks are the only failure modes. 

The non-linear equations of motion mathematical models 
needs to be addressed. As shown in Table 7.3, 7.4 and 7.5, the 
non-linearity of the system caused less than 1% difference 
between the linear and non-linear responses. Under the eleven 
respresentative vessel-drydock configurations, the non-linear 
terms (ie. y6 and Yg®) should be removed from the equations of 
motion. This linearization of the equations of motion will be 
verified in the proceeding discussion. 

The next step in the studying of the non-linearity of the 
system is to determine, when the non-linear terms mentioned 
earlier should be included in the equations of motion. The 
non-linear equation of motion for the one degree of freedom 


model is 


Q 


k = — 
G6 1 -M KG .. MaKG =. 
ee et a = ne es 
T I, I, g I, g 





O + 


mg 


Let ie =a Sinwe 
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M KG a/I, = 8 
A° = kj/Iy 
and 2£\ = C,/I,. 
Substituting and rearranging the Equation of motion gives 
Q + 2EN6 + )26 = -g sinwt(1 -yé). (7s) 
where A = system natural frequency 


B 


normalized maximum amplitude 

Y = ratio of vertical/horizontal amplitude. 
Equation 7.3.1 not only represents the vessel-drydock system 
but any other system than can be reduced to this form. A 
parameter study is performed using 4 8 and Y to determine 
whether or not to use non-linear terms in the equations of 
motion. The selection criterion will be when the non-linear 
response (ie. Y # 0) is different to the linear response (ie Y 
= 0) by 10%. If for a given \’1,8 and y this criterion is 
exceeded, then the non-linear term,Y°, should be included. If 
not, exclude them. For this study, letw = 3.92 rad/sec and & 
= 0.05 to emulate the sinusoidal excitated system with 5% 
damping discussed in Section 7.1.3. Using a Runge-Kutta 
numerical scheme, Figures 7.14 a, b, c, ad and e are generated 
by varying the three parameters. Note that a negative ¥ 
indicates a phase shift of 180 between the vertical and 
horizontal sinusoidal excitation. The linearity selection 
criterion of 10% difference holds true for the region marked 
"N-L" for non-linear for a given ®. In the case of the one 
degree of freedom vessel-drydock systems, the linearizing of 
the equations of motion can be shown on Figure 7.14b since 6 = 


0.5,Y= 0.5 and i ranges from 2.870 rad/sec to 10.254 rad/sec 
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(Table 6.1). 

In order to verify the linearity selection criterion, a 
linear point selected off of Figure 7.14b of 8 = 0.5, Y = 5.58 
andY = 20.0 is implemented into the Runge-Kutta computer 
program in Section 7.2. The nonlinear response equals 
0.015839 radians compared to the linear response of 0.014741 
radians when using the El Centro earthquake input. This is a 
difference of 7.5% which is less than 10%. The linearity 
selection criterion holds true. Now, a non-linear point of & 
= 0.5, \ = 8.327 and y = 15.0 is implemented into the 
Runge-Kutta scheme. The non-linear and linear responses are 
.009652 and .008664, respectively. This is a difference of 
11.4% which is greater than 10%. The linearity selection 
criterion is validated. Thus, vessel-drydock system should be 
checked before selection of linear or non-linear models. 

As mentioned earlier, Equation 7.3.1 can be used to 
validate systems other than the vessel-drydock systems 
subjected to sinusoidal ground acceleration whose frequency 
matches that of the El Centro earthquake, ie. & = 3.92 
rad/sec. this can be done by modifying the appropriate 
linearity selection criteria figure with 

= 8 92 
and replotting the figure. This essentially changes Equation 


7.3.1 into a non-dimensional of 


6 + 2EA, 6, + X72 6, = -8 sinty (1 - Y7 93) 
where 0) = Wu? 
Oy = O/ uw) 
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8) = b/w 2 
Yy = Y 





tC, = %t 


and in this case, ¥ = 3.92 rad/sec. 
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FIGURE 7.1 
Sinusoidal Ground Acceleration 
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FIGURE 7.2 
Vessel Seismic Response in Rotation using Modal Analysis with 
Sinusoidal Ground Acceleration 


= I 








H/3 : H/3 H73 | 


FIGURE 7.3 
Side Pier Forces, with Softwood Cap Block 
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FIGURE 7.4 
Vessel Seismic Response in Rotation for One Degree of Freedom 
Model using Modal Analysis 
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FIGURE 7.5a 
Vessel Seismic Response in Relative Horizontal Translation for 
Two Degree of Freedom Model using Modal Analysis 
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FIGURE 7.5b 
Vessel Seismic Response in Rotation for Two Degree of Freedom 
Model using Modal Analysis 
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FIGURE 7.6a 
Vessel Seismic Response in Relative Horizontal Translation for 
Three Degree of Freedom Model using Modal Analysis 
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FIGURE 7.6b 
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FIGURE 7.6¢c 
Vessel Seismic Response in Rotation for Three Degree of 
Freedom Model using Modal Analysis 
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FIGURE 7.7 
Vessel Seismic Response in Rotation using Runge-Kutta 
Numerical Analysis with Sinusoidal Ground Acceleration 
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FIGURE 7.8 
Vessel Seismic Response in Rotation for One Degree of Freedom 
Model using Runge-Kutta Analysis 
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FIGURE 7.9a 
Vessel Seismic Response in Relative Horizontal Translation for 
Two Degree of Freedom using Runge-Kutta Analysis 
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FIGURE 7.9b 
Vessel Seismic Response in Rotation for Two-Degree of Freedom 
Model using Runge-Kutta Analysis 
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FIGURE 7.10a 
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Three Degree of Freedom Model using Runge-Kutta Analysis 
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FIGURE 7.10c 
Vessel Seismic Response in Rotation for Three Degree of 
Freedom Model using Runge-Kutta Analysis 
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FAGURE 7.L1 
Vessel Seismic Response in Rotation for Non-Linear One Degree 
of Freedom Model using Runge-Kutta Analysis 
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FIGURE 7.12a 


Vessel Seismic Response in Relative Horizontal Translation 
for Non-linear Two Degree of Freedom Model using Runge-Kutta 


Analysis 


= dh ikler= 





| 


Re ON TN AD IANS 


sO 


TAU ie UN) Se CIC N IDES 


FIGURE 7.12b 
Vessel Seismic Response in Rotation for Non-linear Two Degree 
of Freedom Model using Runge-Kutta Analysis 
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FIGURE 7.13a 
Vessel Seismic Response in Relative Horizontal Translation 


for Non-linear Three Degree of Freedom Model using Runge- 
Kutta Analysis 
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FIGURE 7.13c 
Vessel Seismic Response in Rotation for Non-linear Three 
Degree of Freedom Model using Runge-Kutta Analysis 
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Linearity Selection Criteria, 5% Damping, Beta = 0.5 
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Linearity Selection Criteria, 5% Damping, Beta = 1.0 
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FIGURE 7.14d 
Linearity Selection Criteria, 5% Damping, Beta = 5.0 
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FIGURE 7.14e 
Linearity Selection Criteria, 5% Damping, Beta = 10.0 
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8.0 SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS FOR FURTHER 

STUDY 

8.1 Summary 

This thesis investigates four models of the vessel- 
drydock system for seismic analysis. These models are 
implemented to determine vessel response to seismic loading 
for eleven representative submarine-drydock systems. Using 
Equation 4.1.4, the quasi-static force model generates various 
vessel responses listed on Tabe 6.2. For example, the 
submarine in system #1 will rotate about its keel 0.00809 when 
subjected to a strong earthquake. The one degree of freedom 
model predicts a maximum vessel rotation for system #1 of 
0.01471 radians in response to the El Centro earthquake ground 
acceleration history. This rotation of 0.01471 radians is 
verified by the two and three degree of freedom models. See 
Tables 7.2 and 7.3. This maximum seismic response of system 
#1 can also be confirmed by the response spectrum method as 
described in Section 6.3 and reference [1]. The response 
spectrum method gives the rotation as 0.01594 radians. All 
eleven systems have similar differences between the quasi- 
static force model and the various degree of freedom models. 

More important than the maximum vessel response due to 
the El Centro earthquake of .33g magnitude is the maximum 
permissible earthquake acceleration that a vessel-drydock can 
withstand without failure. The maximum permissible 
accelerations and associated vessel responses for the one, two 
and three degree of freedom are listed in Table 7.3, 7.4 and 


7.5, respectively. For example, the one degree of freedom 
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model for system #1 predicts a maximum permissible 
acceleration of 0.0759 g's with maximum vessel rotation of 
0.00339 radians. The two degree of freedom determines a 
maximum permissible acceleration equalling .0759 g's with 
maximum vessel rotation and relative horizontal displacement 
of 0.003358 radians and 0.064566 inches, respectively. 
Finally, the maximum permissible acceleration of 0.0792 g's 
with associated vessel responses (rotation, relative 
horizontal and vertical displacements) of 0.003504 radians, 
0.067373 inches, and 0.017858 inches. The slightly higher 
maximum permissible acceleration found in three degree 
model is due to the beneficial effect of the vertical 
displacement delaying failure of the system. However, as 
shown by the other systems, this is not always the case. In 
summary, the range of magnitudes known as maximum permissible 
accelerations that the eleven systems could withstand without 
liftoff is 13% to 42% of the magnitude of the El Centro 
earthquake depending on system and model used. 

8.2 Conclusions 

The quasi-static method currently used by the U.S. Navy 
for seismic response analysis underestimates the block forces 
caused by an earthquake the magnitude of the El Centro 
earthquake. The one degree of freedom model adequately 
predicts the vessel rotational response to seismic loading but 
is ineffective in the analysis of possible failures to the 
vessel-drydock system. The two degree of freedom model does a 
good job of predicting vessel rotational and horizontal 


translation response to seismic loading provided the response 
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due to vertical acceleration is small, which is the case for 
the eleven systems analyzed. Finally, the three degree of 
freedom model predicts the vessel rotational, horizontal, and 
vertical translation responses and can be used to completely 
investigate various types of failure to the vessel-drydock 
system in all eleven vessel-drydock systems examined. The 
vessel would experience side pier liftoff failure during an 
earthquake with the magnitude of the El Centro earthquake (ie. 
0.33 g's). The various vessel-drydock systems would not 
remain intact during this magnitude of earthquake. Also noted 
that the fourth order Runge-Kutta numerical scheme should be 
used to evaluate system response because of its simplicity, 
results accurancy, and ease of introducing terms into the 
Equation of motion. 

8.3 Recommendations for Further Study 

The seismic response of drydocked vessels needs further 
investigation in three areas. First, a three degree of 
freedom model which will allow vessel liftoff should be 
studied so that other failure mode will occur. Second, a 
study of the vessel-drydock system parameters needs to be 
carried out to explore ways in which drydocked vessel seismic 
response can be decreased. Third, an in depth study of the 
drydock block itself in order to provide failure modes along 
with stiffness and damping characteristics for implementation 


into mathematical models. 
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APPENDIX 1 


NEAR IZED THREE DEGREE ©F FREEDOM SYSTEM RESPONSE USING 
MODAL ANALYSIS METHOD 


imeederetd,ib,1¢,1,mm,n,1job,1z2,ier,hull,.nsys,flagl10 
integer flagl,flag2,flag3,flag4, flag5,flag6,flag7, flag8 
real a(6,6),b(6,6) ,beta,wk(72),weight,h,Ik,gravity 
real m(3,3),cx(3,3) .k(3,3),sidearea,keelarea,plside 
real baseside,basekeel,htside,htkeel,critl,crit2,crit3 
real an(6,6),ac(6000) ,dtau,maxx,maxt,maxy,timex,timet 
Healeyii fiz nro nt) ni2,nti3,ampacc,mass,ampacmax 

real kvs,kvk,khs,khk,base,ht,counter 

real timel,time2,time3,time4,time5,time6,time7,time8 
Peale o00)» t (6000) ey (6000) 1 1(7),s51(7) .pi(7) ,XSCL(6) 
real bbb,ccc,w12,w1,w22,w2,w32,w3,model,mode3,capwidth 
real mmxl,mmangl,mmx3,mmang3,br,amp,plkeel,crit4,ul,u2 
real timey,mmmmm1,mmmmm2,mmmmm3,mmmmm4 

Gamplexeailta(6) 2 (0,6) ,aa(6,6),bd(6,6) ,aa (6,6) ,bb(6,6) 
complex g(6),v(6).,.yy (6.6000) ,ABC (6) ,zt (6,6) ,betal 


CHARACTER*40 XLABEL, YLABEL, YYLABEL, YYYLABEL, DEC 


READ IN VESSEL AND DRYDOCK DATA; VESSEL WEIGHT,KG,I(ABOUT KEEL), 
ee eb eve Nie@resDATA POINTS, VERTICAL STIFFNESS OF SIDE AND 
Pome Plenog HORIZONTAL STIFENESS OF SIDE AND KEEL PIERS, 
GAVITATIONAL CONSTANT,SIDE BLOCK BASE AND HEIGHT, 

KEEL BLOCK BASE AND HEICHT, 

BeCehaereOck Avy pLOeCK-HULE PRICTION COEFFICIENIS, 

SIDE AND KEEL BLOCK'S PROPORTIONAL LIMIT, 

SUDEwE TER VESSEL, CONTACT AREA, KEEL PIER-VESSEL CONTACT AREA, 

CAP BLOCK INCLINATION ANGLE. 


read (44,*) weight,n,Iik,dtau,kvs,kvk,khs,khk, gravity 
read (44,*) baseside, basekeel,htside,hntkeel,ul,u2 
read (44,*) br,amp,plside,plkeel,sidearea,keelarea,zeta 
read (44,*) hull. nsys 

write (6,*) ‘do you want response plots? (y or n)' 
read(5,15) dec 

format (a) 

CO210,71=1,,6 

de-117j=1,6 

a(i,j)=0.0 

b(i,j)=0.0 

continue 

econcinue 

dep 12,1=1,3 

dOnl3,j)-L,3 

m(i,4)=0.0 

Ki(a7 j)=0.0 

exi(1773)) =0...0 

continue 

continue 
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CALCULATE SYSTEM PARAMETERS 
mass=weight/gravity 
beta=asin(sqrt (br**2/(4*h**2) ) ) 


m(1,1)=mass 
m(1,3)=h*mass 
m(2,2)=mass 
m(3,1)=mass*h 
m(3,3)=Ik 


k (1,1) =(2*knhs+khk) 
k (2,2) =(2*kvstkvk) 
Kero) =—((C.o*kvs* (br **2.0))—(welght*h) ) 


DETERMINE NATURAL FREQUENCIES OF SYSTEM 


Pee omni hts) nos) "kK (l, 2) —-m(1,3) *k(3,1)-m(3,1) *k (1,3) ) 
Pmt mm (Ss 23) —m (1,3) *m(3,1)) 
Cee os) —K (1S) *k(3,1)) /(m(1, 1) *m(3, 3) —m(1, 3) *m(3.1) ) 


NATURAL FREQ. MODE #1 


Dele (beech (bob **2.0-4*ccc) ) /2 
wl=sqrt (w12) 


NATURAL FREQ. MODE #2 


w22=k (2,2) /m(2,2) 
w2=sqrt (w22) 
NATURAL FREQ. MODE #3 


Wo7eet pope tsor e(bbb**2.0—4*ccc) ) /2 
w3=sqrt (w32) 


MODE SHAPE #1 & #3 


mogel—(m(1>3) *wlz2-k(1,3)) /(-m(1,1) *w12+k(1,1) ) 
mode3=(m(1,3) *w32-k(1,3))/(-m(1,1) *w32+k (1,1) ) 
Dene Ne. Gil 7 Ci3s,>C31,C3s3 
mmxil=m(1,1)+m(1,3) /model 
mnangl=model*m(3,1)+m(3, 3) 

mmx 3=m (1,1) +m(1,3) /mode3 

mmang 3=mode3*m(3,1)+m(3, 3) 
mmnmmnl=2*zeta*mmxl*wl 

mmmmmn2=2*zeta*mmx3*w3 

mmmmnn 3=2*zeta*mmangl*wl 
mmnmmm4=2*zeta*mmang3*w3 


ex (1,3) =(mmmmm1—-mmmmm2) / (1/model1-1/mode3) 
cx (1,1) =mmmmm1- (cx (1, 3) /model) 
cx (2,2) =2*zeta*m(2,2) *w2 
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cx (3,1) =(mmmmm3-mmmmm4) / (model-mode3) 
cx (3,3) =mmmmm3- (cx (3,1) *model) 


C SE. Urea AND 2 MATRICES FERFORM EICENVALUE PROBLEM, AND 
C DECOUPLE EQUATIONS OF MOTION 
a(1,4)= m(1,1) 
ate) — mis) 
a(2,5)= m(2,2) 
a(3,4)= m(3,1) 
a(s.6)— m(3, 3) 
a(4,1)= m(1,1) 
a(4,3)= m(1,3) 
a(4,4)= cx(1,1) 
a(4,6)= cx(1,3) 
a(sn2)= m(2, 2) 
ats 7>)—7ex (2, 2) 
a(6,1)= m(3,1) 
a(6,3)= m(3, 3) 
a(6,4)= cx(3,1) 
a(6,6)= cx(3,3) 
DG j= -m (1,1) 
Dai) — —-n (1, 3) 
BiZr2)= —m(2Z,2) 
Bie.) = =—m(3,1) 
bie, 3) =. —m(3,3) 
b(4,4)= kK(1.,1) 
b (4,6)= k(1,3) 
Eis, 5)— K(2,2) 
b(6,4)= k(3.1) 
b(6,6)= k(3,3) 


do 100, i=1,6 
domll07s j=1-6 
any j) = Oe (i, }) 
ae) continue 
OC continue 


ia=6 

ib=6 

ic=6 

1z=6 

n=6 

1 job=2 

Cal lmetge © (bia, an,ib,n,1job,alfia,betal,z,iz,wk,ier) 
WRITE (6,*) wk(1),ier 


DO 200,1=1,6 
aoeZo), j=1,6 
zt(j,i)= z(1-3) 

710) ak continue 

200 CONTINUE 
Gall mult (2t,a,aa) 
call multc (aa,z,ad) 
Gall multe (26, byb>) 
Ga emul tGe (bb, 27 ba) 
do 204,1=1,6 
Veer * ba (i, 1) /ad(i-.i) 
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204 


300 


10000 


50000 


(2 i @ et @ | 


S03 


continue 
READ IN ACCELERATION DATA 


do, S00-n—1,, 500155 

reacew(45,-) ac(hyeec(n-.) -ac(nt2) ,ac(nt+3) ,ac(n+4) 
continue 

ESTABLISH FAILURE CRITERIA AND FLAGS 


critl= min(ul, (u2*cos (beta) +sin (beta) )/ 
(cos (beta) -u2*sin (beta) ) ) 

Crit2—Mim (ule) 

crit3= (0.66*baseside-12.0) /htside 
crit4=basekeel/ (6*htkeel) 

ampacc=1.0 

counter=0.0 

ampacmax=0 .0 

continue 

flag1=0 

flag2=0 

flag3=0 

flag4=0 

flag5=0 

flagé=0 

flag7=0 

flag8=0 

flagl0=0 

do 50000 i=1,6 

ABC (i) =0.0 

continue 

maxx=0 .0 

maxt=0.0 

maxy=0 .0 

mm=0 

SOLVE FOR Y,THE COLUMN MATRIX WHICH IS THE SOLUTION IN 
THE IMAGINARY COORDINATE SYSTEM WHERE THE EQUATIONS OF MOTION 
ARE, DECOUPLED 


Gemcol, 1=3,3501, 2 

do 302,i=1,6 

Gia) —(ABC(i)+ac(1—2) /2.54) *exp(v (i) *2*dtau) 

+4.0* (ac (1-1) /2.54) *exp(v (1) *dtau) tac(1) /2.54 

ABC (1) =g (4) 

Nee) (ze ( 4) somp=2t(i,5)+2t (1,6) *h) *mass/ad(i,i)) *q(1) 
*ampacc*—dtau/3 

continue 

mm=mm+ 1 


USING THE MODAL MATRIX ( 4th,5th, AND 6th ROWS ), OBTAIN 
VALUES FOR TRANSLATIONS AND ROTATION 


G6es0s, 1=1,6 

Si (i) =z (4,1) *tyy (1,1) +si (1-1) 

continue 

x (mm) =si (6 

if (abs (x(mm)).gt.abs(maxx)) then 
timex=dtau* (1-1) 
maxx=x (mm) 

endif 

Cor2047 1-16 
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304 


300 


400 


410 


fe) —=2 (On * yy (1s. 1) ar (i= 1) 

continue 

t (mm) =ri (6) 

if (abs(t(mm)).gt.abs(maxt)) then 
timet=dtau* (1-1) 
maxt=t (mm) 

endif 

Gense5, 1=1,,6 

Bigij=z (5,1) “vy (12 1)+pi{i-1) 

continue 

y (mm) =pi (6) 

if (abs(y(mm)).gt.abs(maxy)) then 
timey=dtau* (1-1) 
maxy=y (mm) 

endif 


CALCULATE VERTICAL AND HORIZONTAL FORCES CAUSED BY VESSEL, 
TEST FOR FAILURE 


CALCULATE FORCES ON SIDE/KEEL BLOCKS 
rfl=kvs* ((weight/k (2, 2) )-y (mm) - (br /2) *t (mm) ) 
rf2=kvs* ((weight/k (2,2) )-y (mm) + (br /2) *t (mm) ) 
rf3=kvk* ( (weight/k (2,2))-y (mm) ) 
hf1l=khs* (x (mm) ) 

hf2=khs* (x (mm) ) 

hf3=khk*x (mm) 

TEST FOR SIDE BLOCK SLIDING 


if (flagl.eq.1) then 
go to 400 
Siccmlnmthnil.1t.0.0.and.rfl.gt.0.0-and.abs(hfl/rfl).gt.critl) then 
timel— dtau* (1-1) 
flagi=1 
Slccmimemnicsoc.0.0.and.rf2.gt.0.0.and.abs (hf2/rf2).gt.critl) then 
timel=dtau* (1-1) 
flagl=1 
endif 
x1=x (mm) 
yi=y (mm) 
t1l=t (mm) 
continue 


BES sECR KEES BLOCK SLIDING 


if (flag2.eq.1) then 
go to 410 

elsemlcm i nis.gt.0-0.and.abs (hf3/rf3).gt.critzZ) then 
eimez—dtau* (1-1) 
flag2=1 

endif 

x2=x (mm) 

y 2=y (mm) 

t2=t (mm) 

continue 

TES! FOR SIBE BLOCK OVERTURNING 


if (flag3.eq.1) then 
go to 420 

Sisco) ieee ts o.0-.and.rfl.qt.0.0.and.abs(hfl/rf1).gt.crits3) then 
time3= dtau* (1-1) 
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420 


430 


440 


450 


flags=) 


Glsenit (hfi2.gt.0.0.and.rf2.gt.0.0.and.abs(hf2/rf2) 
time3=dtau* (1-1) 
flag3=1 

endif 

x3=x (mm) 

y 3=y (mm) 

t3=t (mm) 

continue 


TEST FOR KEEL BLOCK OVERTURNING 


if (flag4.eq.1) then 
go to 430 
Glicemiem(rfs.6c.0.0.and.abs(his/rf3) .gt.crit4) then 
time4=dtau* (1-1) 
flag4=1 
endif 
x4=x (mm) 
y4=y (mm) 
t4=t (mm) 
continue 


Ueeperen olDE BLOCK LIFTOFF 


bie(tiage.eq.1) then 
go to 440 

else if (rf1l.1t.0.0 .or. rf2.1t.0.0) then 
timeS5=dtau* (1-1) -- 
flagS=1 

endif 

x5=x (mm) 

y5=y (mm) 

t5=t (mm) 

COne nue 


HeSotegOnencel, BLOCK LIPTOFE 


if (flag6.eq.1) then 
go to 450 

elseelimirrs.lt.0.0) then 
time6=dtau* (1-1) 
flag6=1 

endif 

x6=x (mm) 

y6=y (mm) 

t6=t (mm) 

continue 


TEST FOR SIDE BLOCK CRUSHING 


if (flag7.eq.1) then 
go to 460 

else if (rfl.gt.0.0 
flag7=1 
time7=dtau* (1-1) 


-and. (rfl1l/sidearea) .gt.plside) 


else 1f (rf2.gt.0.0 
flag7=1 
time7=dtau* (1-1) 


eand. (rf2/sidearea) .gt.plside) 
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then 


then 


then 





460 


470 


sO 


e 


60000 


+ 


+ 


JIE} 
i) 


4000 
4050 
4100 
4150 


4200 


endif 

x7=x (mm) 
y 7=y (mm) 
t7=t (mm) 
continue 


DESD er Ok OREER BLOCK CRUSHING 


if (flag8.eq.1) then 


go to 470 

else if (rf3.gt.0.0 .and. (rf3/keelarea).gt.plkeel) then 
flag8=1 
time8=dtau* (1-1) 

endif 

x8=x (mm) 

y8=y (mm) 

t8=t (mm) 

Continue 

Genes nue 


PEGiORE SULTS 

fofe\) else 22.) 

continue 

We(ohe.£O. N’) THEN 
Wimeltetc,. je 1 am finishing.’ 


GO 10 996 
endif 
Weitetoe = 0 am plotting. - 
XSCL(1)=0.0 


XSCL (2) =30.0 — -, 

XLABEL='TIME IN SECONDS' 

YLABEL='ROTATION IN RADIANS' 

YYLABEL="RELATIVE HORIZONTIAL DISPLACEMENT IN INCHES' 
YYYLABEL="RELATIVE VERTICAL DISPLACEMENT IN INCHES' 

Cb eeOet eine (7, 1s00,OXSCL(XSCL) ~OISCL (—2) ,OXLAB (XLABEL). 
QYLAB (YYLABEL) , QLABEL (4) ) 

Gree tere (1717 1500, OXSCL(XSCL) ,OISCL (—2) ,OXLAB(XLABEL) , 
QYLAB (YLABEL) , QLABEL (4) ) 

ereeer Perth (1,1, 1500,OXSCL (XSCL) ,OISCL(=2) ,QXLAB(XLABEL) , 
QOYLAB (YYYLABEL) , QLABEL (4) ) 

Gouto,z0000 

CONTINUE 


if(ampacc.eq.1.0) then 


write (46,4000) nsys 

hemmawidleaw-cox, **** Svstem ‘,12,1x,'*****) 

write (46,4050) hull 

hemmed / oO, **) Hull ~,13,1x,'*** ) 

write (46,4100) 

foumatiix,//,2e6x, * Ship Parameters *°) 

write (46,4150) 

format (1x,/,5x, 'Weight',8x,'Moment of Inertia'’,9x,'K.G.") 
write (46,4200) weight,!Ik,h 

fommat(ix, 19. 151%, kips >ix,fl1i.1,1x, kips-in-sec2', 


Foxtel LA, ins ) 


4250 


write (46,4250) 
founae('s,//,20x, * Drydock Parameter *°) 
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4300 


4350 


4400 
4450 


4470 


4475 


4500 


4550 


4600 


4650 


4700 


a7 50 


eto 


4780 
4800 


4850 


4875 


e995 


4990 


4900 


200 


S000 


S8fs}8) 
|100 


write (46,4300) 

format(1x,/,1x,'Side Block Height',3x,'Side Block Width’, 
+3x,'Keel Block Height',3x, ‘Keel Block Width’) 

write (46,4350) htside,baseside,htkeel,basekeel 

fogmatiex to. oc ins 711x,to.), 1x, ins’ ,11x,f6.1,1x,‘ins', 
+9x £6... ,1x, ins} 

write (46,4400) 

format (1x,/,1x,'Side-to-Side Pier Distance’) 

write (46,4450) br 

feommacilx, 716.151. ins ) 

write (46,4470) 

format (1x,/,' Total Side Pier Contact Area’ 
+,3x,'Total Keel Pier Contact Area’) 

write (46,4475) sidearea,keelarea 

temiac(i< ox cll ix dnd » 14x, flil.dix,'in2*) 

write (46,4500) 

fogarty ios lock—-on-block Eriction Coeff ’,3x, ‘Hull-—on-Block 
+ iP hleciom COoett ) 

write (46,4550) ul,u2 

format (1x, 10x, £7.3, 23x, f7.3) 

write(46,4600) ~ -- 

format (1lx,/,1x,'Side Pier Proportional Limit',3x,'Keel Pier' 
+7 Proportional Limit) 

write (46,4650) plside,plside 

fegmac (ix, 10x717.35,1x, Kips/in2'15x,f7.3,1x, ‘kips/in2") 

write (46,4700) 

femmaw ix. /-ix, sice Pier Vertical Stiffness’ ,3x,'Side Pier‘, 
+" Horizontal Stiffness’) 

write (46,4750) kvs,khs 

formar tlie ix, Kips/in ,11x,f1l.1,1x,'kips/in') 

write (46,4775) 

format (ioc,/,1x, Keel Pier Vertical Stiffness’ , 3x, 
+'Keel Pier Horizontal Stiffness’) 

write (46,4780) kvk,khk 

Pomiatielc oii lex, Kips /in’, ix, f11.1,1x, ‘kips/in') 

write (46,4800) 

format(1lx,//,20x,'* System Parameters and Inputs *') 

write (46,4850) 

feymac(ix7, ix Inout Eorcing Function is Horizontal Component’, 
+' of the 1946 El Centro’) 

write (46,4875) 

format (1x,20x,' Earthquake Acceleration Time History.') 


write (46,4995) 

format (1x,/,1x,'Vertical/Horizontal Ground Acceleration Ratio' 
+,3x,'Data Time Increment’) 

write (46,4990) amp,dtau 

fommatiix,1Ox.t6.39,to5,16.3,1X, sec’) 

write (46,4900) 

format (1lx,/,1x,'Gravitational Constant',3x, ‘Percent System 
+ Damping’) 

write (46,4950) gravity,zeta*100. 

femiatCh-w*. to.2,1x, in/sec2',10x,f6.2,1x,‘'%') 

write (46,5000) 

fommat (ix,/7,2>5x, ‘Mass Matrix’, /) 

do 5100 i=1,3 

Wiemce (46,5050) m(i,1).m(i,2),m(i,3) 

Poumaciioe 115.4,/5x,f15.4,Sx,f15.4) 

continue 

write (46,5200) 


oa Oe 





5200 
5)/25)8 
23500 
5400 
5450 
212) 0/8, 
5700 


6000 


Lommel ./,25%, Damping Matrix -7/) 
dowosoe 1-153 

Wimnce (46,5250) “come ) , ex to2)7 ex (1,3) 
Hommat (ix, f15.4,5x,f15.4,5x,f15.4) 
concinue 

write (46,5400) 
femmaciix,/ 25x, Sceitiness Matrix ,/) 
do S500 i=1,3 

weate (46,5450) k(i71),k (1,2) ,k(i,3) 
ftonmat (1x,f15.4,5x,115.4,5x,f15.4) 
continue 

write (46,5700) 

format (1x,//) 

WRITE (46,6000) 

FORMAT (1X, 'Undamped Natural Frequencies',t35,'Mode #1',t50, 


+'Mode #2',t65,'Mode #3') 


<(8)0)a5 


+! 


6002 


write(46,6001) wl,w2,w3 

Pomiactix tol tieox, Tad/sec ,t4o,f7.3,1x, rad/sec',t62,f7.3, 
rad/sec') 

WRITE (46, 6002) 

FORMAT (1X, ‘Damped Natural Frequencies',t35,'Mode #1',t50, 


+*Mode #2',t65,'Mode #3') 


6500 


ens 


O500 


a 


25000 


WRITE (46,6500) wl*sqrt (1-zeta**2) ,w2*sqrt (1-zeta**2) ,w3*sqrt (1-zeta**2) 
fom tix7t al 1) eo7Ix, Tad/sec ,t46,f£/7.3,1x, 'rad/sec',t62,f£7.3, 
rad/sec') 

endif 


write (46,10500) ampacc*100 
fopmactix<.///71x, Ber Earthquake Acceleration of °,£6.2,° % ' 


uOtetne Els Centro ,/) 


write (46,25000) 
foumatiix, Maxiamums/Failures ,t26,°X (ins) *,t36,"Y (ins)',t5S1. 


+'Theta (rads)',t65,'Time (sec) ') 


B10 
olalpl 


au Z 


aS 


wr ite (46, 25001) 

Bema tt ————————=—= cco, ===> >= TR a a as Picis)3 
------------ ',t64, '-----------') 

write (46,310) maxx,timex 

femmat (ix, “Maxamum x ,t25,f£9.6,t65,f5.2) 

write (46,311) maxy,timey 

fommatenix, Maximum Y ,t3S,f£9.6,t65, fS.2) 

write (46,312) maxt,timet 

PeoGuatlioce aMaximum Rotation ,tS0O,f£9.6,t6S,f5.2) 


Peevmagi.eq..) then 
flaglO=flagl10+1 
write (46,313) xl,yl.tl,timel 
Pomiatm(ixa aide block sliding ,t25,f9.6,t35,£9.6,t50,£9.6, 


+t65,£5.2) 


314 


endif 


wemcrlag2.eq.1) then 
flagl0=flagl10+1 
write (46,314) x2,y2,t2,time2 
hommes. Meel olock Ssiiding ~t25,.f9.6,t35,f9.6,t50,f9.6, 


+t65,£5.2) 


endif 


if (flag3.eq.1) then 


=o 6> 





315 


316 


Si7 


318 


Suh] 


320 


11000 


flagl0=flagi0+1 
write (46,315) 
format (1x, 'Side block 
HEOS, fo. e) 
endif 


if (flag4.eq.1) then 
flagl0=flagl0+1 
write (46,316) 
format (1x, 'Keel block 
+EGO 71.2) 
endif 


iieCflags.eq.1) then 
flagl0O=flagiO+1 


write (46,317) 
format (1x, 'Side block 
EO Selo. 2) 
endif 


if (flag6.eq.1) then 
flagl0O=flagl10+1 
write (46,318) 
format (1x, 'Keel block 
+t65,{5.2) 
endif 


Him flags.eq.1) then 
flagli0=flagi10+1 
write (46,319) 
format (1x, 'Side block 
TEGO, 15.2) 
endif 


if (flag8.eq.1) then 
flagl0O=flagliO+1 
write (46,320) 
format (1x, ‘Keel block 
mtos tows) 
endif 


if(flagli0O.le.0) then 
write (46,11000) 


Xo, Vou.rs,cimes 
Ovemreunming. 4 tzo, f9.6,t35,19.6,t50,£9.6, 


x4,y4,t4, time4 
SvenrtuUnning —-tzo,19.6,t35,19.6,t50,f9.6, 


x5,y5,t5,timeS 
liftoff’ ,t25,f9.6,t35,f9.6,t50,f£9.6, 


x6,y6,t6, time6 
liftoff’ ,t25,£9.6,t35, £9.6,t50, £9.6, 


Kx7>yV?,t7 time? 
erushingae,tz5, £920, t35,19.6,t50,£9.6, 


x8,y8,t8,times 
Emuciing meet elo-to-6,Co5.19.6,t50,19.6, 


HOnmatiion/,1%, No failures occurred. ‘) 


if (counter .egq.1.0 
go to 60000 
endif 


if(counter.eq.0.0) then 


ampacmax=ampacc 
ampacc=ampacc+.1 
counter=1 .0 
write (6,*) 
endif 
endif 

if (ampacc.le.ampacmax) 


if(counter.eq.1.0) then 


ampacc=ampacc-.01 


sand. 


flagl0.le.0O) then 


‘I am in the secondary looping stage.' 


go to 20000 


else if (counter.eq.0.0) then 


ampacc=ampacc-—.1 
endif 


tll 8) Ie 


20000 


300 


200 
100 


300 


200 
100 


Geace 10000 
continue 
stop 

end 


subroutine mult (a,b,cyy) 
real b(6,6) 

complex a(6,6),cyy (6,6) ,a(6) 
go 100, i-176 

do 200, j=1,6 

do 300,k=1,6 

@ (k) =a (i,k) *b(k, j) 

continue 

cyy (1, J) =d (1) +4 (2) +4 (3) +d (4) +4 (5) +4 (6) 
continue 

continue 

return 

end 


subroutine multc (a,b,cyx) 

Gomplex a(6,6).,b(6,6) -cyx (6,6) ,a(6) 
do. 109, i=1,6 

do, zZ00, j=1.6 

do 300, k=1,6 

Qk) =a (17k) *b(k, 3) 

continue 

cyx (1,5) =d (1) +d (2) +d (3) +d (4) +d (5) +0 (6) 
continue 

continue 

return 

end 


es 





CCN GO GY Cra GO) 


IS 


Wg) 
abZ 


APPENDIX 2 


NONSDINEAR THREE DECREE OCF EXREEDOM SYSTEM RESPONSE 
USING FOURTH ORDER RUNGE-KUTTA METHOD 


integer NN,1,mm,n,hull,nsys, flagl10 

integer flagl,flag2, flag3, flag4, flag5, flag6, flag7, flag8& 
real beta,weight,h,IkK,gravity 

real m(3,3).,cx(3,3).k(3,3),sidearea,keelarea,plside 
real baseside,basekeel,htside,htkeel,critl,crit2,crit3 
real ac(6000) ,dtau,maxx,maxt,maxy,timex,timet 

mealenli niZ2 nis nil hi2 nis, ampacc.mass, ampacmax 

real kvs,kvk,khs,khk,base,ht,counter 

real timel,time2,time3,time4,time5,time6,time7,time8 
real x(6000) ,t (6000) ,y (6000) ,XSCL (6) 

real bbb,ccc,wl2,wl,w22,w2,w32,w3,model,mode3,capwidth 
real mmxl1,mmangl,mmx3,mmang3,br,amp,plkeel,crit4,ul,u2 
real timey, mmmmml1,mmmmm2,mmmnm3,mmmmmn4 

Peano dao (5), (5) >C(5),D(5),5 (5) .F (5) ,G(5) ,HH(S) 
CHARACTER*40 XLABEL, YLABEL,YYLABEL, YYYLABEL,DEC 


READ IN VESSEL AND DRYDOCK DATA; VESSEL WEICHT,KG,I (ABOUT KEEL), 
PMEeNeKEMENT OF DATA POINTS, VERTICAL STIEENESS OF SIDE AND 
Ropero HORIZONTAL STIFENESS OF SIDE AND KEEL PIERS, 
GAVITATIONAL CONSTANT,SIDE BLOCK BASE AND HEIGHT, 

KEEL BLOCK BASE AND HEICHT, 

BECCK OI pEeGh AND BECCK-HULL ERICTION COEFFICIENTS, 

SIPESAND KEEL BLOCK S PROPORTIONAL LIMIT. 

SUObeETER=VESoEE CONTACT AREA,KEEL PIER-VESSEL CONTACT AREA, 

CAP BLOCK INCLINATION ANGLE. 


read (44,*) weight,h,Ik,dtau,kvs,kvk,khs,khk, gravity 
read(44,*) baseside, basekeel,htside,htkeel,ul,u2 

read (44,*) br,amp,plside,plkeel,sidearea,keelarea,zeta 
read (44,*) hull,nsys 
write (6,*) ‘do you want response plots? (y or n)' 
read(5,15) dec 

format (a) 

domilz 71-73 

de i3; j=1.. 3 

m(i, 4) )-0.0 

k (1,j)=0.0 

Ex 12 3) =0 0 

continue 

continue 


CALCULATE SYSTEM PARAMETERS 
mass=weight/gravity 


beta=asin (sqrt (br**2/(4*h**2) )) 


bee 





m(1,1)=mass 
m(1,3)=h*mass 
m(2,2)=mass 
m(3,1)=mass*h 
m(3,3)=Ik 


k (1,1) =(2*khs+khk) 
k (2,2) =(2*kvst+kvk) 
MS) = (eso *KVs™ (or **2.0) )—(weight*h) ) 
DETERMINE NATURAL FREQUENCIES OF SYSTEM 


Boome hoo) timo o) “KCl, 1) —m(i,3) *k(3,1)-m(3,1) *k (1,3) ) 
Vine) 3.3) —-m (1,3) *m(3.1)) 
Geena) ) (3,4) kh (les) *k (3,1))/ (m(1,1) *m(3, 3)-m(1,3) *m(3,1)) 


NATURAL FREQ. MODE #1 


Wi2—(>p>b—saqrt (bbb**2.0—-4* ccc) ) /2 
wl=sqrt(wil2) 


NATURAL FREQ. MODE #2 


wee—e(2,.2)/m (2,2) 
w2=sqrt (w22) 
NATURAL FREQ. MODE #3 


Wwi2—(—Dob+sart (bbb**2.0-4*ccc) ) /2 
w3=sqrt (w32) 


MCDE SHAPE #1 & #3 


moadel— (mits) *wie—-K (1,3) )/7 (—m (2,1) *wil2+k(1,1)) 
Mode s— (mili 5) Woeek(1, 3) ) 7 (—-m(1,1) *w32+k (1,1) ) 
DETERMINE Cli, Cis 631 ,C33 
mmxl=m(1,1)+m(1, 3) /model 
mmangi=model *m (3,1) +m(3, 3) 

mmx 3=m(1,1)+m(1,3) /mode3 
mmnang3=mode3*m(3,1)+m(3, 3) 
mmmmml=2*zeta*mmxi*wi 

mmmmnm2=2*zeta*mmx3*w3 

mmmmmn3=2*zeta*mmangl*wil 

mmmmm4=2 *zeta*mmang3*w3 


cx (1, 3) =(mmmmm1-mmmmm2) / (1/mode1-1/mode3) 
cx (1,1) =mmmmmi1- (cx (1, 3) /model) 
exs(2-2)—2*zeta*m(2,2) *w2 

cx (3,1) = (mmmmm3-mmmmm4) / (model-mode3) 

cx (3, 3) =mmmmm3- (cx (3,1) *model) 


READ IN ACCELERATION DATA 


dO 7s087 n=1,5001,5 


=140= 


300 


10000 


3000 


meaduit5,*) ac(m),acimt.)  ae(n+2Z) ,ac(n+3) ,ac(n+4) 
eeneanue 
ESTABLISH FAILURE CRITERIA AND FLAGS 


critl= min(ul, (u2*cos (beta) +sin (beta) )/ 
(cos (beta) —u2*sin (beta) ) ) 

crit2=min (ul,u2) 

crit3= (0.66*baseside-12.0) /htside 
crit4=basekeel/ (6*htkeel) 

ampacc=1.0 

counter=0.0 

ampacmax=0 .0 

continue 

flag1=0 

flag2=0 

flag3=0 

flag4=0 

flagS=0 

flag6é=0 

flag7=0 

flag8=0 

flag10=0 

maxx=0 .0 

maxt=0.0 

maxy=0.0 

mm=0 
x(1)=0. 
y (1)=0. 
t(1)=0. 
R=0 .0 
S=0.0 
TAU=0 .0 


ao © 


IMPLEMENTATION OF EQUATIONS OF MOTION INTO THE 
RUNGE-KUTTA FORMULUS 


@o 30277 1=1,6000 
DO 3000,11=0,4 
A(11)=0. 
Bill) =). 
C(11)=0. 
D(11)=0. 
E(11)=0. 
Bit) =O". 
G(11)=0. 
HH(11)=0.0 
CONTINUE 
mon=mmn+ 1 
DO 302, NN=1,4 
MEANNG EO. 1) THEN 
FE=0.0 
foe Pee nO. 2 .On. NN.EO.3) THEN 
FE=0.5 
PESE 1f (NN.EO.4) THEN 
FE=1,.0 
ENDIF 
A (NN) =dtau* (R+FF*D(NN-1) ) 


OO0O0000 


B (NN) =dtau* (S+FF*E (NN-1) ) 


= 4 


soZ 


one) 


C (NN) =dtau* (TAU+FF*F (NN-1) ) 


D (NN) =dtau* ( (-cx (2,2) /m(2, 2) ) * (R+FF*D(NN-1) ) - (kK (2,2) /m(2, 2) ) 
+* (y (mm) +FF*A (NN-1) ) -amp*ampacc*ac (1) /2.54) 


Ciinv=aGau~((-ex( 10) /mii, 1)) * (tbr *E (NN-1) )—(cx(1,3) /m(1,1) ) 
+* (TAU+FE*EF (NN-1) )-(k(1,1) /m(1,1) ) * (x (mm) +FF*B(NN-1) ) -ampacc*ac(1)/ 
+2.54) 

HH (NN) =dtau* ( (-cx (3,3) /m(3, 3) ) * (TAU+FE*F (NN-1) ) — (cx (3,1) /m(3, 3) ) 
+* (S+FF*E (NN-1) )—(k (3,3) /m(3, 3) ) * (t (mm) +FEF*C (NN-1) ) + (m(3,1) /m(3, 3)) 
Pe tecoa. 27m (22) ) “(R+re *D(NN—1) )—(k(2,2) /m(2,2)) * (y (mm) +FF*A (NN- 
+1))) * (t (mm) +FF*C(NN-1) ) —(m(3.1) /m(3, 3) ) *ampace*ac (1) /2.54) 


E (NN) =(m(1,1) *m(3, 3) *G(NN) -m(1, 3) *m(3,3) *HH(NN) ) / 
Hmis os) aml, 1) (1,3) %m(3,1)) 

F (NN) = (HH (NN) - (m(3,1) /m(3, 3) ) *E (NN) ) 

continue 


DETERMINING SYSTEM RESPONSE 

y (mm+1) =y (mm) + (A (1) +2*A(2) +2*A (3) +A(4)) /6 
x (mm+1) =x (mm) + (B (1) +2*B(2) +2*B (3) +B(4) ) /6 
t (mm+1) =t (mm) + (C (1) +2*C (2) +2*C(3) +C(4)) /6 
R=R+ (D(1) +2*D(2) +2*D (3) +D (4) ) /6 

S=S+ (E(1) +2*E (2) +2*E (3) +E (4) ) /6 


TAU=TAU+ (F (1) +2*F (2) +2*F (3) +F (4) ) /6 


MAXIMUM VALUES FOR TRANSLATIONS AND ROTATION 


if (abs (x(mm)).gt.abs (maxx)) then 
timex=dtau* (1-1) 
maxx=x (mm) , : 

endif . 

if (abs(t(mm)).gt.abs(maxt)) then 
timet=dtau* (1-1) 
maxt=t (mm) 

endif 

if (abs (y(mm)).gt.abs(maxy)) then 
timey=dtau* (1-1) 
maxy=y (mm) 

endif 


CALCULATE VERTICAL AND HORIZONTAL FORCES CAUSED BY VESSEL, 
TEST FOR FAILURE 


CALCULATE FORCES ON SIDE/KEEL BLOCKS 
rfl=kvs* ((weight/k (2,2) )-y (mm) — (br /2) *t (mn) ) 
rf2=kvs* ((weight/k (2,2) )-y (mm) + (br /2) *t (mm) ) 
rf3=kvk* ((weight/k (2,2) ) -y (mm) ) 
hf1l=khs* (x (mm) ) 

hf2=khs* (x (mm) ) 


Sie 





400 


410 


420 


430 


hf3=khk*x (mm) 
TEST FOR SIDE BLOCK SLIDING 


if (flagl.eq.1) then 
go to 400 
Clcemime(itl.lt-0.O.and.ril.gt.0.0.,and.abs(hfi/rfl).gt.critl) then 
timel= dtau* (1-1) 
flagl=1 
elcemurm(ni2.gc.©.O-.and.ri2.gt.0.0.and.abs (hf2/rf2) .gt.critl) then 
timel=dtau* (1-1) 
flagi=1 
endif 
x1=x (mm) 
y1=y (mm) 
tl=t (mm) 
continue 


TEST SPOR KEEL: BLOCK SLIDING 


if (flag2.eq.1) then 
go to 410 

eleerit irf3.gt-0.0.and.abs (hf3/rf3).gt.crit2) then 
time2=dtau* (1-1) 
flag2=1 

endif 

x 2=x (mm) 

y 2=y (mm) 

t2=t (mm) 

continue 

TEST FOR SIDE BLOCK OVERTURNING 


if (flag3.eq.1) then 
go to 420 

Srscemi te detec O.amd.1t1.gqt.0.0,and.abs (hfil/rf1).gt.crit3) then 
times= dtau* (1-1) 
flag3=1 

Glcemiucn( ni .ge.0.0c7anad.ri2.gt.0.0.and.abs(hf2/rf2).gt.crit3) then 
time3=dtau* (1-1) 
flag3=1 

endif 

x3=x (mm) 

y 3=y (mm) 

t3=t (mm) 

continue 


TEST FOR KEEL BLOCK OVERTURNING 


if (flag4.eq.1) then 
go to 430 

else if (rf3.gt.0.0.and.abs (hf3/rf3).gt.crit4) then 
time4=dtau* (1-1) 
flag4=1 

endif 

x4=x (mm) 

y4=y (mm) 

t4=t (mm) 

continue 


Geeieook. oLDE BLOCK LIFTOFF: 


=143- 





440 


450 


460 


470 


301 


60000 


ife(tiacgo.eq.1) then 
go to 440 
elsepuemintl.1teOoGr.on. fle-.1c.0.0). then 
timeS5=dtau* (1-1) 
flagS=1 
endif 
x5=x (mm) 
yS=y (mm) 
t5=t (mm) 
Eontinue 


DES WpOnecen lL: DEOCK LIP TOrE 


if (flag6.eq.1) then 
go to 450 

ClseslGrtGrs. 1t.0.0)) then 
time6=dtau* (1-1) 
flag6=1 

endif 

x6=x (mm) 

y6=y (mm) 

t6=t (mm) 

continue 


TESTerCk SIDE BLOCK CRUSHING 


if (flag?.eq.1) then 
go to 460 

else if (rfl.gt.0.0 .and. (rfl/sidearea) .gt.plside) then 
flag7=1 
time7=dtau* (1-1) 


else if (rf2.gt.0.0 .and. (rf2/sidearea) .gt.plside) then 
flag?7=1 
time7=dtau* (1-1) 

endif 

x7=x (mm) 

y7=y (mm) 

t7=t (mm) 

continue 


TEST BOR KEEL BLOCK CRUSHING 


if (flag8.eq.1) then 


go to 470 

else if (rf3.gt.0.0 .and. (rf3/keelarea) .gt.plkeel) then 
flag8=1 
time8=dtau* (1-1) 

endif 

x8=x (mm) 

y8=y (mm) 

t8=t (mm) 

continue 

continue 


PLOT RESULTS 
go to 999 
continue 


=A a= 





Shs hs. 
Boo 


4000 
4050 
4100 
4150 


4200 


4250 


4300 


4350 


4400 
4450 


4470 


4475 
4500 


4550 


4600 


4650 


+ 


+ 


+ 


i£(DEC.EQ.'N') THEN 
werce(o, )) 1 am finishing.’ 
GO TO 998 
endif _ 
write(6,*) ‘I am plotting. ' 
XSCL(1)=0.0 
XSCL (2) =30.0 
XLABEL="TIME IN SECONDS' 
YLABEL="ROTATION IN RADIANS‘ 
YYLABEL='RELATIVE HORIZONTIAL DISPLACEMENT IN INCHES' 
YYYLABEL='RELATIVE VERTICAL DISPLACEMENT IN INCHES’ 
CALL QPICTR (X,1,3000,QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) , 
QYLAB (YYLABEL) , QLABEL (4) ) 
CALL QPICTR (T,1,3000,QXSCL (XSCL) ,QISCL(-2) ,QXLAB(XLABEL) , 
QYLAB (YLABEL) , QLABEL (4) ) 
CALL QPICTR (Y,1,3000,QXSCL(XSCL) ,QISCL (-2) ,QXLAB(XLABEL) , 
QYLAB (YYYLABEL) , QLABEL (4)) 
go to 20000 
CONTINUE 


if(ampacc.eq.1.0) then 


write (46,4000) nsys 

femal, .2ox, **** Svstem ",12,1x,'*****) 

write (46,4050) hull 

Hommaciix.7 90x, ** Hull ©,13,1x,‘'**") 

write (46,4100) 

format (1x,//,28x,'* Ship Parameters *') 

write (46,4150) 

format (1x,/,5x, 'Weight',8x,'Moment of Inertia',9x,'K.G.') 
write (46,4200) weight,Ik.,h 

hOmuat (ioctl el- ix, Kips: ;-ix,f11,1,1x, 'kips-in-sec2', 


Tox tol, ix, ins ) 


write (46,4250) 

format (1x,//,26x,'* Drydock Parameter *') 

write (46,4300) 

format (1x,/,1x,'Side Block Height',3x,'Side Block Width’, 


+3x,'Keel Block Height',3x,'Keel Block Width’) 


write (46,4350) htside,baseside,htkeel,basekeel 
Commie ex ins ,lix,f6.1,1x, ins ,llx,f6.1,1x,'ins', 


+ Oceano. ioc, 1S ) 


write (46,4400) 

format (1x,/,1x,'Side-to-Side Pier Distance’) 
write (46,4450) br 
fowmac(ix,t7.f6,1,1x, ins’) 

write (46,4470) 

format(1x,/,' Total Side Pier Contact Area' 


+,3x,'Total Keel Pier Contact Area’) 


write (46,4475) sidearea,keelarea 

homme co, tid dx, in2’,14x,f11.1,1x,"in2") 

write (46,4500) 

formatiix,/,ix, Block=cn—Block Friction Coeff’ ,3x, ‘Hull-on-Block 


+ EGrction Coeff") 


write (46,4550) ul,u2 

hermuaw lx, 1Ox 67. 2,23x,f7.3) 

write (46,4600) 

format (ix,/,1ix, Side Pier Proportional Limit‘’,3x,'Keel Pier' 


7 eenopert ional’ Limit ') 


write (46,4650) plside,plside 
Poguarm(ixlOxete.o,1x, kips/inz 15x,f7.3,1x, kips/in2‘) 


ho 


4700 


4750 


47715 


4780 
4800 


4850 


4875 


4995 


4990 


4900 


4950 
5000 
5050 
5100 
5200 
5200 
x10) 6) 
5400 
5450 
5500 
5700 


6000 


6001 


6002 


6500 


+! 


write (46,4700) 

format (1x,/,1x, ‘Side Pier Vertical Stiffness',3x,'Side Pier', 
Horizontal Stiffness‘) 

write (46,4750) kvs,khs 

fommacie sx fla bce hips7 in. lix,fil.1,1x,.'kips/in‘) 

write (46,4775) 

format (1lx,/,1x,'Keel Pier Vertical Stiffness', 3x, 


+'Keel Pier Horizontal Stiffness') 


write (46,4780) kvk,khk 

penmamidewsx~, tli. 1 Kips/in’,l1x,f£11.1,1x, *kips/in') 

write (46,4800) 

format (1x,//,20x,'* System Parameters and Inputs *') 

write (46,4850) 

format (1lx,/,1x,'Input Forcing Function is Horizontal Component", 
of the 1946 El Centro') 

write (46,4875) 

format (1x,20x,' Earthquake Acceleration Time History.') 


write (46,4995) 
format (1x,/,1x,'Vertical/Horizontal Ground Acceleration Ratio' 


+,3x,'Data Time Increment’) 


write (46,4990) amp,dtau 

eeumat (ia 10x, f6.5,to5,.16.3,1X, sec’) 

write (46,4900) 

feqmaciix./,.1x, Gravitational Constant ’.3x, ‘Percent System 


+ Damping’) 


write (46,4950) gravity,zeta*100. 
Ponmawmere 7x 1o.2,1x- in/sec2 ;10x,f6.2,1«,°%") 
write (46,5000) 

format (1x,/,25x,'Mass Matrix',/) 
de5100° i=l 3s 

Wieetwe (20,5050) m(i,1).m(i1.,2),.m(1, 3) 
fomnar (ix, flo.4,5x%,f15.4.5x,f15.4) 
Gentil nue 

write (46,5200) 

format (1x,/.25x, ‘Damping Matrix',/) 

do 75300 i1=1,3 

write (46,5250) cx(1i,1),cex(1,2) .¢x (4,3) 
Femrmat(ix,f15.4,5x,£15.4,5x,f15.4) 
continue 

write (46,5400) 

format (1x,/,25x, ‘Stiffness Matrix", /) 
de S500 1-1, 3 

write (46,5450) k(i,1),k(1.2)-k(4, 3) 
RomMaciilx:, f15-4,5x f15.4,5x,f15.4)} 
continue 

write (46,5700) 

format (1x-,//) 

WRITE (46,6000) 

FORMAT (1X, 'Undamped Natural Frequencies',t35,'Mode #1',t5S0, 


+'Mode #2',t65, ‘Mode #3') 


+! 


write (46,6001) wl,w2,w3 
Pemmiaciix,tal,f7.,o,1x, rad/sec ,t46,f7.3,1x, rad/sec',t62,f7.3, 
rad/sec') 

WRITE (46,6002) 

FORMAT (1X, "Damped Natural Frequencies',t35,'Mode #1',t50O, 


+'Mode #2',t65,'Mode #3') 


+! 


WRITE (46,6500) wl*sqrt (1-zeta**2) ,w2*sqrt (1-zeta**2) ,w3*sqrt (1-zeta**2) 
Honmatdlx tol, t7eo,1x, rad/sec ,t46,f7.3,1x, ‘rad/sec',t6é2,f7.3, 
rad/sec') 


ai aG— 





10500 


+, 


25000 


endif 


write (46,10500) ampacc*100 
format (1x,///,1x,'For Earthquake Acceleration of °,f6.2,' % ' 


fon the Elecentro 7) 


write (46, 25000) 
format (1x, 'Maximums/Failures',t26,'X (ins)',t36,'Y (ins) ',t51i, 


+'Theta (rads)',t65,'Time (sec) ') 


S13 


write (46, 25001) 
fonmael 4 ———— = — TAPAS a coo = - =a oe 3 .cSOy 


foumat (ix, Maximum x ,t25,f£9.6,t65, £5.2) 

write (46,311) maxy,timey 

feqgiarm( ix.) Maximumey ,t35,19.6,t65,1£5.2Z) 

write (46,312) maxt,timet 

fermata(ix. Maximum Rotation’ ,t50,£9.6,t65,f5.2) 


Mem flagi.eq.i1) then 
flagi0=flagl10+1 
Wisltce (46,3135) xl,;yl,tl1,timel 
Pomnoeexeoide block Sliding ,t25,19.6,t3S,19.6,tS0,f9.6, 


+tGo fo. 2) 


314 


endif 


iio et lag2.eq.1) then 
flagl0=flagl0+1 
write (46,314) x2,y2,t2,time2 
fommacmixne\eccl block Sliding ,t25,19.6,t35,f9.6,t50,f9.6, 


+t65,f5.2) 


Sls 


endif 


if (flag3.eq.1) then 
flagl0=flagl0+1 
Wirtew (465915) x2,Y3,t3,times 
fomialml< side block overturning ,t25,£9.6,t35,f9.6,t50,f9.6, 


+t65,£5.2) 


S16 


endif 


if (flag4.eq.1) then 
flagl0=flagl10+1 
write (46,316) x4,y4,t4,time4 
feourarmixs eel Block overturning ,t25,19.6,t35,f9.6,t50,f9.6, 


+t65, £5.2) 


o17 


endif 


if (tlags.eq.i1) then 
flagl0=flag10+1 
write (46,317) xS,v5,tS,timeS 
fomMawil oldesplock Jiftoff {t25,19.6,t35,f9.6,t50,f9.6, 


+t65, £5.2) 


218 


endif 


if (flag6.eq.1) then 
flagl0=flag10+1 
write (46,318) x6,y6,t6,time6 
foemaue( ion kKeeleblock liftoff  (t25,19.6,t35,19.6.t50,f9.6, 


+t65,£5.2) 


co OT 





endif 


if (flag7.eq.1) then 
flagl0=flagl0+1 
Write (46,319) x7,y/,c7,time?7 
Sule, foumormt side bloeksemuch ing s.t25,f9.6,t35,f9.6,t50,f9.6, 
PLOD toc) 
endif 


if (flag8.eq.1) then 
flagl0=flagl10+1 
wrace (467570) xo- yS,t8,ctimes 
320 fepmoatelix, Keeleolcek crushing ~t25,f9.6,t35,f9.6,t50,f9.6, 
765, 15.12) 
endif 


if(flagi0.le.0O) then 
write (46,11000) 

11000 format(ix,/ ix. Ne tallures occurred. ‘) 
if(counter.eq.1.0 .and. flagl0O.le.0O) then 
go to 60000 
endif 
if(counter.eq.0.0) then 
ampacmax=ampacc 
ampacc=ampacc+.1i 
counter=1.0 
write(6,*) ‘I am in the secondary looping stage.’ 
endif 
endif 
if(ampacc.le.ampacmax) go to 20000 
if(counter.eq.1.0) then 
ampacc=ampacc-—.Q1 
else if(counter.eq.0.0) then 
ampacc=ampacc—.1 
endif 
geneo 10000 

20000 continue 
stop 
end 


= deer 





APPENDIX 3 

MODAL ANALYSIS OF THE TWO AND THREE DEGREE OF FREEDOM SYSTEMS 

Since the vertical equation of motion is uncoupled from 
the other equations, the three degree of freedom system need 
only to be analyzed to obtain damping coefficients for both 
degree of freedom mathematical models. Also, the maximum 
response of the two systems will be determined using response 
spectrum analysis with participation factors described in 
Section 6.4. 

The three degree of freedom equations of motions, 
undamped, as shown in Equations 4.4.lla, b and c are 


My] x 3 3 m,3 6 + OLEL x = -M)) Xg 


M2 ¥ + ko2 Y = -M22 Yg 
and 
m4 + m3, xX + k33% = -m3] <a 
where ™)] = M2 =M 
M}3 = M3, = MKG 
m33 = I, 
Ky1 = 2Ksph + Kxh 
Ko2 = 2Ksy + Kxy 
k33 = (B2/2) Key - W KG. 
To perform modal analysis, consider the free vibration system 
m7 X + m3 6 + ky] x = 0 
Moo Y + Koo y = 0 
and 
M33 8 + M3, X + k33 9 = 0. 


Assume the system response is in the form 


=149= 





= >) sin(t + P) 


va 
| 


y do Sin(wt + 0) 
and @ = $3 sin(wt to). 
Now, the equations of motion are 
Mj, (-dyw* sin(wt + o)) + m3 (> qwsinWwt + o)) 


+ LSI (Py sin@wt + 0 )) = 0 


Mo (-d5v2 sin(wt +0)) + ko5 (65 sinWwt +p)) = 0 
and 
M33(-$3 w% Sint + 0)) + m3, (42% sinWwt +0)) 
t+ kag (3 SinWtE +p ))= 0. 


The trivial solution to the above set of equations is 


Sin(jwt +0) = 0. Assuming sin(Wwt +o) = 0, then 
2 Z 
‘Si ae ee ee 4 : 
Z = 
2 yD oh 0 al a Oy}. (A.3.1) 
2 zZ 
—M 5 0 K337-M, 30 $3 0 


These equations are solved for, by setting the determinant of 
the first matrix equal to 0. However, the second row of the 


matrix is uncoupled. Thus one solution for is 


w* = ko2/m22- 


Now, the determinant equation reduces to 
eG ess 33 Ve) 1133 
F anes > ee 
iiiewsee Sl 31 Asha eee cauec 


Using the quadratic formula, the other values of y* are 


A.3.2 
ia Ga ( ) 


UG) = 
Z 


=o O= 





where b 


- (a 
= ease Shae ta ae 


i oe 


and 63 == SS 
hl 338 13, 31 


Using parameters found in Section 5 for the eleven representa- 

tive vessel-drydock systems, the natural frequencies wj, wo, 

and w3 for each system is calculated and listed in Table A3.1. 
The mode shapes are determined with the relationship 


(Ky17 - Mz w*) o2 - my? 3 = 0 
3 ASS a i 
With mode shapes determined, the equation or motion can be 
rewritten, including coupled damping 

M1 ¥ + m3 65 + cy, X¥ + c43 6 + kyy x = 0 

M22 Y + Co2 y + kg2 y = 0 
and 

M33 0+ M3] x + C33 8 + C3, X¥ + k33° = (O- 
At natural frequencies “,, “5 and”3 the coupled three degree 
of freedom system acts as a Single degree of freedom system 


with the following relationship: 


(04/03) x = 9 
(61/63) X = 6 
and (9,/%3) oS 68, 


Using these relationships, the three degree of freedom 
equations of motion at the three system natural frequencies 


are 


=—i>1= 





[m1 + (¢3/61)m 3) X + [Cz] + (¢3/61)¢31] * + ky, x = 0 
M22 ¥ + C22 y + koa y = 0 
and 
[($4/65)ma7 + m33]8 + [(67/63)057 + €33]8 + kg = 0. 
The damping coefficients for the problem are 5% of the 


critical damping for each mode, 


[C71 + (¢3/63)C13) = 2 EM) wy , w, = 1,3 (A3.3a) 

Coo = 2E M2 w2 (A3.3b) 
and 

(Con Ore C3) + C33) = 2 M3¥y, ¥ = 1,3 (Ase 3c) 
where ~ = percentage of critical damping, 0.05 

My = ™11 + (63/61)n ™13 


and M3 = m33 + (61/¢3)n ™31- 
Substituting known values and rearranging Equations A3.3a and 


A3.3c yields the following four equations: 


Mode l, Cyi t+ (93/61) 13 = 26M) 4] 
(61/03)1 C31 + €33 = 26M3 6] 
Mode 2, C11 + (63/41)3 C13 = 26M) %3 
ae) aecsietscss — 2° M3 > 
or C172 + (03/9 4)1 C13 = 26M) %] 
C11 + ($341)3 C13 = 2 EM, 3 


and (1/93)1 S31 + C33 = 26M3 0] 


(1/6 3)3 €31 + 33 = 2 E€M3 03 

Using the above set of equations and equation A3.3b, damping 
coefficients ¢c)], C13, C29, ©€31, and c33 are determined and 
listed in Table A3.2 for the eleven vessel-drydock configura- 


tions. Note that c)3 equals to c3),. In the two degree of 


freedom case, the damping coefficients Cj}, C19, C21, and C292 


Se 





are equal to those in the three degree problem, i.e. 


Two Degree Case Three Degree Case 
tal ea 
C12 SALE 
©21 S31 
C22 33 . 


Using the response spectrum method and modal 
participation factors, the maximum response of the two and 
three degrees of freedom systems can be calculated. Starting 
with the forced vibration equation 

[M,7 + (63/93) n ™3) x + Cc) xX + ky, xX = -mj} Xq 


or eos 3 ee -m 
1 11 -_ 11 P 


Pit (¢3/ 01) 3] Pit (63/41) ™ 3 


Equations 6.3.1, 6.3.2 and 6.3.3 predict a maximum response 


for a 5% critically damped system for a given mode as 


pein 


eC fOr 8), 0, < 2.24 
pax es na 3 
or Caen). - 26.03 , = s 7 = mOie 22 24) < Wy 1 W:; a 2 7, 
“ | 11 (3/1) 5 a 
Seb, 58: AGL 
ee aa) OD 1? 8a > 12-74 
max'n Wye pa tad 1 3 


Similarly for rotational response, the equation of motion is 


[m33 + @1/¢3)n M31) 9 + €39+ k339 = -m3] XQ 
or 


C20+k 36 “M34 


Q + = . 
ee oa arate 33°77 %3)5731) 9 


=—23> 





and the maximum response for given mode is 


m 
31 
(6 jae— 2262 nua fs for Ww diet 7 
max’n [m+ ( 1 3 n™1] i ees 





m 
Z6-.03 oyll 
foes (Oo  ) = ° pep eee Bee! <I E 
max’ n Ww Im. 3+(o, 3) ,M31] 1 3 
2M Ss zc! 
ns a Ef Pie aN 
else a ane 7 EO Was ”. pee 4 


[m53+(%)/%3) m4] 


Since being uncoupled, the vertical equation of motion is 


Mao ¥ + Cp ¥ + ko2 Y = -MQ2 Yq 
or es Cee y, os ‘ 
y + Sass = -y_ = -AMP x 
22 g g 
where AMP = ratio of vertical/horizontal El Centro earthquake 


accelerations. 


The maximum vertical response is 





(Ymax)2 = 11-62 * AMP for % < 2/24 
or 
(Ymax)2 = nO 20 * AMP for 2.24 < % < 12.74 
else 
Wane) = 228s * AMP for W5 > 12.74. 
W 


Now, the maximum response for each mode of the x and 0 
equations and for the uncoupled y equation are calculated. 

The configuration of the three degree of freedom system 
at any time is a superposition of the two coupled natural mode 
shapes along with the independent vertical maximum response. 
The absolute maximum response to a given earthquake is the 
numerical summation of the maximum response of each mode shape 


times its respective participation factor. The general 


=)54> 





formula for each modal participation factor is 
Ph =lMrorn/ Mr 6 “rn 
where [yn = modal particpation factor for the nt? mode 
r~ refers to which equation of motion is being 
considered. In the three degree of freeddom 
case, r= 1or 3. 
Using these particpation factors, the maximum response of the 


three degree of freedom system should be no greater than 
Xmax = (Xmax)a 2 + (%max)3! 3: 
Ymax = (Ymax) 2 

and 6 max = (8max)1 [1 t+ (8max)3! 3: 


These maximum response are determined and listed in Table A3.3 
for the eleven vessel-drydock configurations. Once again due 
to the uncoupled vertical response, the maximum response of x 
ands in the two degree of freedom is identical to that of the 


three degree of freedom case. 


=2005 





SYSTEM 


ui 


0 
ita 


NOTE: 


TABLE A 


Se... 


Three Degree of Freedom Vessel-Drydock System 


5 


3 


4 


2 


2 


2 


7 


Gre 


D 


Ga 


6. 


Natural Frequencies 


Wy 
.454 
.883 
.013 
eat 
.848 
690 
.953 

976 

863 

857 


177 


(rad/sec) 


W 2 
20.071 
36.823 
29.468 
Ng Pe ak 
Bon 737 
45.450 
ML Os 
39.247 
32.829 
34.832 


33.293 


W 3 
B72 a0 
34.367 
23.791 
2 Sal 
32.742 
40.863 
36.351 
34.045 
29.060 
29.921 


Zon OO 


Two Degree of Freedom System Natural Frequencies 
correspond to ©, and 3. 


=O = 





TABLE A3.2 


Three Degree of Freedom Vessel-Drydock System 
Damping Coefficients 
(kim-sec/inch) 


SYSTEM Cil C92 Can as eon 
at 65.398 166.826 1322882.125 3920.416 
2 58.612 157.226 938583.250 2874.701 
S 42.452 125.824 976657.438 2829.964 
“ 67.008 as SoS 682485.813 2055. 760 
5 54.303 152.588 687504.500 Zoos 17 
6 28.057 445.695 7879946.000 ens Om 227 
7 167.100 405.212 6430934.500 14269.135 
8 155.786 384.866 5622059.000 12660.506 
9 43.963 116.475 1196244.250 3340.509 

10 55,65. 86.437 836179.438 2408.758 
dd Sal 632 82.618 751062.813 2182. 306 


NOTE: Two Degree of Freedom System Damping Coefficients 
correspond to Cj)j, ©€13, ©€3 , and C33. 


SIS iS 





SYSTEM 
a 


Z 


10 
tak 


NOTE: 


Iwo Degree 
correspond 


me 
ae Seen 


Oo 3275 


0.27141 


0.65142 


0.53788 


0.21468 


0.45271 


0.46833 


0.46597 


0.64366 


0.65864 


0.67467 


TABLE A3.3 


0.10859 


0.12225 


0.19089 


0.21491 


0.12980 


0.08025 


0.09708 


0.10762 


0.15381 


0.13662 


0.14955 


= /30— 


Three Degree of Freedom Vessel-Drydock System 
Maximum Response Using Response Spectrum Analysis 


(its.) 


(rads. ) 
0.01646 
0.02348 
O02 550 
0.02350 
0.03202 
0.00683 
0.00866 
0.01009 
O01 25 
0.01342 


0.01522 


of Freedom System Maximum Response 
eonx and 0. 











DUDLEY KNOX LIBRARY ~ - 
NAVAL POSTER, NATE SCHOOL ~~. 
MONTEREY MYR ITA 959045-BOK 






































































































































8 oe. 
e ~ -, “ ¢ s@m- sa - Ps 
. s : ¢ s “s z Ld 
nd ¥ © of « + « fe ”_ 7 a rae 
rie “ 2 by be . °=“ ase rive ° 
“se * = ° vere - a eee 8s 8 « e lew os a - ° 
¢ o o ’ e ae e Sin 8 ® %e a . Py o - eee - a a Ca 
: o F " ve . e 9 le air rials a a on es ° ® ’ ete * 00g wee > *. © ° 
« « * a . ; bie ae Ld * . a s ° S 2 oe 8 ° To ve . . a ° oe 
= . : k s s 2 Z er yr . o * ~. ars a ot ss e x » a v er) 1 Py . wee mee a A od) @te 
: - Sie: Taine . 8 . - 2 2 » *neee - e. . e e Mg : x pi RN S = ss : cla oA sete AM : 
ES 52 oat ee es 7 Poe 7 z z i i . Es sa : - E : ’ ae ’ = = é ° es nw ew o Pa . 4 , ’ ’ . 7° ole ae © @¢o te « * 7 or 8 ° oe 
; Ps <a = . ~ eS cas ee ren = ¢ = = : ve + ¢ . hl - ( . ° . . Py Py . © . e é o @ete Pe e - Pet - Po fF was & oy bad = Cd 
~ ONO OS) Sige. gg , = on . a ee tow o e 7 ° Ai = ® 3 . 2 = : = Z i sty Se ese hae <2 Tama! eee) ‘a x ie ce. Tee ee eared a es ce ay Gales 
as ae be = Peet m o a ee. ee 2 « " 8 F o = = > fi ita = - ‘ cx ’ a * ° eof * @ vt we . ¢ ’ “ve a . * s "F Het e mms “2° o- 
so tm. on - eae” Jt Oe I ay - 5 ° re -.¢ < ° - ro ca oars = aa ae - - Nai = Pee bs » nf 5, * Cd ae fad s a . é Ce rere eats base se é¢e@ * am te ® * ¢ o Meets 2 reo ° ~- beset oof Ls = * 
eet ° . axe Pre same OP Os om. - : . * ’ 2 A be asa Pa bsg Sil Sted £ le SEE 5, be all ° Binds Ye eee ene enwiale 
Sete Sect eee Roa ail ee ERE pe ee OPIS eee oen ae r <6 , Cer ¢ * oe te os e , oF 6m. o “ear  F oz Kae opis * Ps muse we "me Ber ne . vo 6 
¥ aes eae: au ; , ee eae Saree =: ae ‘ .? ap a eNisms Ok P a - « P ™m « wa -_ rr ree s Ti Gentes Se OS Pats aa - ES 24 Cae wei a) eee.) s wet oe FEE 0 al 6 nem 
Sates icine ou 6 rs : a neta z of ’ ‘ ® - m@ ¢ ~ ° we ee - 5 : > , fe - we . ’ dn * - oe * . 
nies oa melee eee atari Sea Aes So , wail, ° de pes pr = Pins eo ae bal ad . cs een == fp 2 Wsy F Pa - Y « ‘ . ° © ane or s,s re os ’ rs oe ° ry ae 4°0 FUE or 050 £88 e Ooh. owe - 
; 2 a ait j fs ice en ¢. Pay oo oe oe ae i oe tae. pe. scarily *. = as e » an # < fue 7 ve ~ wie #8, . . « -* s oe « * ve . . a r oe s.¢ * *¢ vee 
pe ~ Be — A Pas = : * 2 = aire oo % 2&8 = 7 @=g « w ® ad =3 * « rd ow - 7 owe aoe - ee ae - - z ze ‘ as P e oer ° ; ‘ Fi dal : a ao she - # ‘* wee e rlee ews = eee 
Ke = cree Seles BAe a er i eal isn « - ai ae « ae -  « - - “- Pare e e re * a : - = iis an * f ps “4 i = hinge ed — gi pad eres © 
= = : oe pia: my wat“), (SSiere eA - eta We oy . wow - - as oo, , 7. . is vob pe i re” eielietaat oc es Mca, 
- - “- uo me -.° "ee - se - e MOO 8 wey on sere ea Sy eee eee Srna an : i ht . 7 ee - ee ag 7 Lf id d a =! - fe co o ‘ o wee were s * peels ee Roirsege ee 
2 « ee ee -- eta ein . © , ies ge: spears Bt oie ie A k i . 7 pel ze Che eat BEN aol slaeeie : "2 s . id = aC} ae es 
a e be " eae a = teties <j dc: ? ath ~ . % * ¢ @ a © s @# s@« - * * . ad . “gi Ss : i eee rere g e =a 7 4.8 
° ae e eae Sy ere > e nm ae ee = ee a me acute re -. > : é i, rd ~“*s ’ . . or . - BO or ae. & *« Pre ¥ £ @ “Hwee « OP te eran a Aa a 
- 7 oman o= ™e Ne eure eo .. +e * a ee mire te is o, Pe ae cere pateoes Se ae ay ee . prs - ‘i Set eee yeas, ae ee es PD er ete Ls : ot te Pee aw wy : © fe 4 we te ery S ¢ a i al ode 
a ee = “ ee aca ee am Sia) AMS lwig atte a ee: posers a is ae, Ae ate Bis ne, ep. . Ld pe . Te ses I erie we owe . rae Meh « me ts OU ee None wm | ae PY ts to o« 6 re "se STs i gL ah shea 
< . =) es ail 6s ae ece omen « vos 2 ie = : ¢ oy Fue 8 9M . . s# - a a ee . rng om x “ive é ° td 0 Le over wl fo 
epee . siete Rite ‘ ead pl aS = gece idl 7 i= wien Beaten nis Seer a. ee ‘: wis - a ar o fer = om acoeies *, © we a SOA > a , Satis <a. ere ciee eee mate eh aes avs he AeA epee « te rime e: : a e @: . aga _ Pee ie “ges is ae = = o - aif a 
ee ok “Nf Bw y oP Be te wee e+e OPS OA MMe op a im, @ F abt ue A <0. ae aie an Med Pong Oe ! dea OCR ams fe tls =e ie Cr aa ‘ eas 2 ONE Ot eng , pie oe IT PT a YO ~— oe Ome y oe 98 PRs ee Me TH 8 Ve ws is ia ends es 
- So Ses Me eS eae pee cain 17m me my dose, Biel pias mer ice "48" wwewn.. pi mptecelsacar, anita. Ne Rie, e A eAS eo . . Py . - e * +4 Pe @as Winn ¢ . Pres 6 ee CS sh g, ‘ ornlenn TN al p tp 5 ap o E- > Ru we Po e eee. fe 8 iy we elas een * 
ee se - 2 @ -~« FOR Cee aed ok Tt ee 7OF Am Lem. ere 7 ose ene er ee . - me ° . « 5s A < fed bes ~ as ° . wey va * mu " ae ae Fal) pane T° OER gg Fig fare ver: wy OF POs £54 +49 2 7! Sw of ae oo"? fee wo ow 10 o-at & ge v8 
ow, . New @ « ~ «® Neg ate ees "vee Oo Le wisem os. ee "e @: Mas boop es Seeuene z r *. re sah sk wwe ae ove es ole 0 sews, a, & . . “@ oP mse On ° a Evae oy Wis al ylegsy gas wee eae = 8 ry) we OOP ere ge oh Pe Mtee- 8 5 4, 
SF Nes ve?.4 = ~ ~~ os me --. & eon = = a . P oe eae ‘es cane * = one eistone Se zich = ae) r - «4 v,., 6 ore Huy el rte se Ios - a a © 7 ee eons « g Sae Te, oe * a oe ate o®. 1 2 ee Were fr ow 0 oe age lg ane LEP Magee 
=. . see, - = Ce stiee wes eM Fo) 8 wt rat o. So VOM eM Aimee wf os ine 26 we a ee ri wide = ba A bé ° el oy Toa - « « - ye om. ° Fee ¥~« o © owe ° He 6 gogswes YwKre we Hae So OSE PS) Sere ty aon gn el et nT Per tog Cogs 4 
ote oe + Ne em we te eels ° Me ¢ we Aer rs m\ ree oe j agar S 2 eS 2 TUTE see 8g the og a a or Yr 7 PP ¢ . oe wi = &— NOM Pe. . -  -##* @ oo” ae AF © twee yeee Oe Tr erne s Pe ge ge WON IT a8e® 08 ma 80, Sere 
a= i a Aig Nm Fo, owes ee Me ae — 1% om eas re il ' ee ne * Mle Spe ie 4 ae whee) eid e) Sten -@¢ Oh «6 oe alee ee > 7 ¢ 6 6€8T 8 6 7mm, *&re a a 0.0 gag rem 4 6 a ee ee a : s - “y Veale wwe ° oo 7 = os 
8 oe © © we ee e a ee em) eiacds A OWvesee we Cage tc Bias eee - i) oO we . ed et Pre Oe Nee aoe. . 6 ~ ow a ay i wae ot ony Se Aes soe e oes Spe ges <a % CROP RY FP e we, OFC + ow weErey 5 cg a oe. 
ce py °F 00 +. ower 2 Per “ent Ae iw, oF 4 we CAME RAN fw ow, i “ pate irs a a hes ais ~ wh 28 He ve al! mee = re me °° re FONG Ww tw ewe * LENE Se Fee prgey ng ome weet 8 em Hy gle te es 
po mee ween i Seats elt etapa as k oe. Skat eet eee Sree ke Pendle ape Y 46 tte. ° 4 oss --e WOrmrtrs 6 2s rPrem , Fe OT PITA wir of $0 Oretag FO. 8 Me 6 yp reeg Oi wore ase ep <0 Ue OP KE we A Oo dm -. 
Se ana iyi ana ee ee Figiaobah Bagel. $A meni nee og bad ¥ Y on, vow 9 Es form pe ee te Kee we 6, e Owe we comen - Me seve Cte We ee 8 eee wre « Fe pg, 
eke US oye dd Noreen aver o.0 Mee wg ® 6 we tmeg ee ; a gate. ~ Sareea a ie - * fiw w-ree ss wee PER e fF sews « OCR nae Poon ies . vs owe FE wee sag a g.d -y OPC OF EO FORE ge, F Sele wwe OP ow pean ces 
mae me ew om * "Fe Come o, . Pe moe s we iter remye 07ers ovneRempie tye. e 4 s td * we z - es * #om ove r VEE» eo wees v ~ Mee wuopae ye s opr ee eg te one ug tee ot A Ae ea PME” we el i) te 7 Oe nev ie 
© woe a. eo *y © Fete eye POP my, 3% 20 ws SOU tee re 8 re CPE. e Bey ot even Nee 's oe Mer: e . Pw cee ag FP OP ey Srna ee Oe! Ue alee enl ae ar OF seer te 8 oy Lt aes CAMO O wh we wiwy “ OF OF FE cere ton oe NEE ev seg "own re ee 
FTO DP Wosage Poi 2 e.g 7928.86 we * - «. on Mow — a ex Ne OMe — ome 2 , S*8 Pum se oe ee a ee? in -“ Ce od a * eet or, oef a o Pan wee Y ft ome eg Fee Metaleot wg AP eee ese 
7 AF.9 0.0. nnn orem we hy PP Fo Em Oe Smee, on Woe. 3 4g O. wevegue er a sy pid sea es ee : ele . gs . es © 4 SPR Ae wen, ". Ne Gad UNE UE ot er ee oo. PE ve We APOE otters OPMENT AF OER OT By Bau gr » 
Fue eye oe Flee -“ =m 7 ere § os * © wa Ww. FERED CR me nn Pi cra ; ; tae 5 ps eae ‘ ¥ : eerie tim Ft. ete rer ay Te 1 wy me oe ee fC 8 a wt te ? fen ey ¥ tf e «ye oer ae EWE} rags rEg ahy onraad = 
* wae OFF a4, Pathan OF = Wes FRY EP iecre oi 2, ee ve #y gam sr? wise *« wy » » aoe nf Capt , ad 2 PN PF et 2 ty Ln ae et = we « 0) 6 ee é Ries a a o eet rte x oe ew or gne oo ala til Adel ad” eet oO Tory awe Paes « TOTO HF Pe t-0 wget ? 
2 Gam F “Svcrs » 72 = _ a nd 5°? 4 ‘. Fe ‘ ve Hg ey re eo — <2 o oP %ene er, ‘ pied) se dpe 7 Pw ew Pie eer mp emo eC «Fetes os « ha tl Cel °F Pe re Woe « Faw VETO § eo hs ¢ Cee ee oe A oe ope OF! i eit wwe FOO TP we or oo eee Lag a al a ‘ae wigs 
teltiae ace og 5 con, ana ae rs ges wea we cate Md , owe ywurar “Fes em « ra rw i a a ary bog Fuss g¢ ayvew, “¢ Senses aPhcat! o a ee ot ee ve VU ow ateln VS Fr ee 8 Soars P°rEF Me te Oe wr ine £ee8 °F a ey te» 
© Pisa! MO pas nw. ge mat an Navmieccss ee eu meres ae one 7 paid hee ree pias Z ce fae, * wae == bf LOR Reem © 0 Fae « Cisvyue Ow od F ay gages °F trisa i -£e @e PL ERMAN ie 6 o wis + o- rhe @ , @ Re ewe FP OBIE eg ee fume - *¥e erg 
re we 7%, Dam 208, Few wen oe a. ity Mt Po ge 0 © Rm ge oneem es 4 re Oy AAAs oie ioe - Z 4 ane ae e e = € or ore SP P Free ww we wo, POON Foe Asa Po 6 Oy rr 4 ue op o FORE 6 ® os pes Fre ON. wo co weg 22 ao SPs a8 egng sige a mee Oe a oe ot gmge ge 
7 ie mw dt OWN mw raw, 5a. Heh on ected ° fu. * Pe ENA tnece ne 0 -IwEwER. ve tm ge fee if as Peta aa me CO ee me anes 0 oF 88 ans a Fhe 0 9 we +s ibid, a ae a Tree "A NOES Pavey sg S| weg Peeve oO Oe TNE es me ete @ op en te fee 
MN CP mama Rp Fy oe gee ae ae s oc Se OP FA! Oa Nees nh 4S OF te ih wo 2 ee Fe ert) ipa <_-«-._ ; aleheniaiaicge wig te thee og? e Wed eee Pel siniglalwiels wae 7 Wo AeNs mune wine o,f Pg OF 22 wee y WV gle 2092 zt yr 
Ss) meean sgn wien scl ota = CR a arene! ne Saat a elid pinion ies AX 2 ae ae eineeNesleMGlaren geass 59 ts we eo WAP ser s aw "» © MW sur ge ge See CO RE emir Vite whet ¢ Fe tagms ¢ Teer sea x oP 8 RO OE EY ay a OVO EPO & a i tea el tod - a me 
EA ee Oe os me i =e NORETNE Os = 8 oy 40%, Pee, aw, = 0 OP Oe ren 6 «0 ay ee = “in FW os TAL ae P negt “pas Si ww “oa 7 et 6 Ns At oy we we. ara ae ‘ a ¢ fo ey eel ek) Car aie aoe Pet 0 Hs ot CaS a8 ng OP 80-  wuiwe e oy Le ARE x ee 8 
+E nay at we are eon Ow “ROA row. oF ose eWvee > PR ee Le Wie yee x Pies or s # eee Hy ele & on ’ “Ie eg - nee Me e Tre Maa Dl a er a a “s- ft ate 9 ORPe=sWme Bg ult Gnd A et ae oe 
NF GDA OO. POM ow 78 eee egy TM, we WAL. és A he Oy Oem. ceorn Noe My farina pte potas rie Roni tcaeae Rstiie AB id " a SEM NC Neca ug St Tea ir sa ee 8 OF aie COP TS © One go we ote Oe we bo EF ae KL Neng Tes NI a Mc LAE 
NP NEGRO in wes < POP Mee en ne qi wee pe Suga? Socal eee Aan aeek me oe se “ahaa P-Mee eter 6 e ow SCONE OUP 0 te TI et fgg £ 5 8 oP sty 5 (2 F aguatete 4 w a 0°02 0 © eee PU ep py FO OO N69 8 a atne Hee IE Oe Fd ay 
ei 8/ ES’ Py Oe Nae enn =. hee ok Ca ee Ce ae 1h need of hs Tene e s a - E a ais Ld wg Vm say eI, vese es Fe Pm am ee 0 oe” ire © ggg UAE INS POH 0-0 great PF & gle : SNe oe oie og = me OE FOG eG LE gow roaatabetons 
Oe wees te &e ~ eve, Reale ee Lorca om wre ight AP ata eee ie ara MUSE -tinmwern ce ne ns : Reeds be ae pine it atl ot Mee re “ r'faa % BOWE Hs Crlg » we ’ Ee OP mates 5 faewerveg LEE 8 28a o Bees dpe CRORES DN Sa > wre Foor 
* Pe we ing “Us 2 WORM Ol DO see om PC @ ae rnin? wg me Bee ee ee Ue Ocean SA 0, : eho SOW * pad » . vere rw ~ teu es Bode. eue © atl Dee ee, ee rik \_MEL vo se wee P aww as ? a oD PHL Fe oF He a oe fe we 2° PM £00 2 ay aw Le 0 ee hogs # toe FUE OY & pn! Use yas FL © segs gy 
MP ON 2 Ob ee oe Far ek °» PAO ONG uses he tl cee, SIE aE ops TRicNacOl Gc, can , . ere. oR yee in FIRE lr ewe, on wee f duets o7 0 1s 0dN0 =n ory f Weeegtonuee-sa0y gag Pie-nee rae why She meas eters see PHONE 8 OP SPN PE Pig presen sees oF 
Sitio eeie ia F Nant © Repel eee . ww Po i Me 9 <n ete atts on phi cas pe “wm om MOE 2. OO Sig ae ty 4 bo Stearn i Sea et oR Bert SAT) Pe ete AMPs GOP, cele ap “ONE an Be “ «NOS dope og Reg | en e™ Fas Aarnsia > ine Pomel Midis 0s wd game Oe UE & O° pa 0g a0 intial belated ak dl ok del 
fF - mim ee og Se ee ees Fda eee Solo aen ve Maelo face et Mw ae ene eaeaae a maeeerene = s aie el . cere um bck a “ia =e ROP rage my eR a SKIT oe Pg Z NTO al aad de eee PD die ti ywe We try eos 8 whale we oye P MP--ihaid Megha ot Stat ee PE Sar ge 0 « ae ed ts bite 
Poet rou, eo mm earas a ? erouk ial dh i = : . - uv - bellacnad 9 ~ = Rey “APRN eee oN vais 5 “6 SR p UC ewee ge Z Pt “P06 8 pag Cy a Ne fem, Hy ~ « . Pte igs , o-Ps o - « SERVER eH EE ey Oe Peles re 
mu? pape oe ae we ee Eee ea ae Seeman a perpnk Lap afc * wat inna ae © FEF ES .c 0s Soi UO Haare ately De ~™ w- % o. eh le nd ee Oi rey tebe OOE™ 98 OX gabyeny : “We so Pure bine Pus eps a on ws ae see . ve esp e wwe 70 oe wr eee * tb spake agp Cs ae weatare: safe Paphos =p Pet gma Pee Or ys 1 ree nee omg | 
MeN te | ce fi Ame ca acdc W505 oe RN Re. of, pA eee Sr nt 4 9 aan So Oni ikl Ave eres ones SLE Se * ine “Ee ieee ee ay wes 1. ere A caer Wem tons ige ee ue OO. PSE a PL OLINE 6 4 ‘ MP TS suMisinens datece Praha ee Ah WP oF ah egre once ee Te TAINO 8 OPE SE bye 98 a8 Cay ye ps FW nee ee UTR talent 20 ate 
See Naeem aN my cena ers hi. he. won ge Oo> Pe mergeaery rn ee Wiwinceneeee ™ babe ed Talon beac NO uNe own TEC eA. Olivy Memwie ee. Pio lewis ewe y SUE ee SP ey ap + OEE arege gn SN Tw Be FIN ony WNP e vee ome he erie MM g= 5m 0 ni data cn I af deol hak ae apie dK shasta ied ated ae sphagt totale 
Be awe osu ayes - mA re hans 7 We mom Wi WO etl mee Ol mee. > © Femhiee ee By aye an NY wear on,” = r wl om = ms - sd ~- @vy a wa hh he er re FF idee wees a % . oe OT re iy ve ri) Pm. vem wy the Os Tn ome PRE oy © “et hove & sioaye OLE 0 t's Pw gt Oe LE o oP ye pares go g We Far Pee eet om “a SP weg 
“MA elw eee, cr 2900 See © sao FN MI, se oie Mette ean te Peieies § Oem Ka. a ae Sere pita PP dite Dea, SreGe © 8 Ges tueeihge a oa CASS NGA care WWTre Tut ° fe died tay Portes Vets Fe oy ne tne oe, gt NOOO A MMe 6. poe Fa Eo Ee ON ere pig » thekeliied che ee 
©. mee ee aN? Se ee SEAGRASS Pitan ads. te 72h Swag og Sag WA gwee Fe eee eee epee lady peated Ldn ee ay eee Seay COLOR) SP ating %, Tie oe ‘ NENT? OWE mein ee 8 5 gnoys Meese, ar oO) Lote te YE JP ae tras Be nes ¢ 8 Oe») Minny nde Od we energy PEF OND: FO oo Fea OP as we ee an 
2-2 Pee ele. ee Meee Achat sagt 9 Tet edhe su, saat bl ae tN ows ee Tg Oy XEON am WP ie SEPae ny i a ¥ pan rf sob es on oo ON Ow Re Re NET one OVIRE 6 Gf ON PONS 02.9 dN pap “VA o cw att od an a ee OF CAL OS 6 6° whee we ETEN ey ee MUO ges) LAF 
ao te, aa ee wats wvite an ta. oe biases, CM OTT a NOE oi dr gil - EOIN IS: oF ers wy ee Nee Nsee Spr a om ve ANNs. 6 ING te Weg Cw bie tly tops ad ON UND ro SO pe oy bitedaiubeatete’ Mall aoe Ae hited ait inital ot hed abn 
ena - ON OF Uh Hip yng an, CONN ewe, Wiss ain Mas te Ne Lot a FREs OS tte, gy hit) Abed et ART . - FP he ww Pre oeiye ol Soe $e 8g 18 Fe Crs ¢ rere OO 148 ot ice At eal a ee : s - WF yrot lig E> ge ge Oe Fray Medd hd dd A ee POO FOR FOE vy € 
‘ Sore se? Copa hla diedidhiiherkel st aliod at 2 he PUR eee 6 ny - sly 0 PNEE Sa. oye, - . I at Le ae SEO NOTE YE Ae goary oF oe =f Lee: cpr Wire ne pig heli ded nts mada oan Pao OSCE 8 Eee f 
itp lent Yee <2 hay . nie oe rs meh ac tah alts a fo weer sty a Saas ay wy -7os 5 PN Es “Up eerey OFT 6 Pham iP Woot 5s rn gos, P oy me» lh ddiokt dnd ton edhe ae Se Oye Se Peng ST gry FIOM POM FOE Gwe ge y 
—_ WPL ty sae ena Kg, EN oe Uriel mentation aie reciae See a eee pees lmhnie 7 INI ON EEE 8 tate 8 3 al . . oN of He 8 S Lene 80 mo Le pugeatig ¢ We og eg s ninth dipped aan Ne ee Frege EO Oe te Fy PAs 
= De mee PEN SR! tres -4Ue On ok EY OL TREN Sait e ayy Dy, ame Ep ae: weave eee rey un ae, pipidateten, Fe Otte 6 bet sity . Oe “ieee LOE FS Peg! 8 HONG: OF Ee Fi 4- ee Na et wb canes ge rn Nf oRe- he Pe Pas een ep of Gone ee Trine 
io) — = - pe NEAL! ey meen PEP wre Tea Pree MAS OS ee Lg wwrery Wyeen tee een be » - tas tt aise. ie btbeatch “sie * Letinedh od) SO a PF OER FP ee Cte © 6 preeshs t VOT pS Keng CMM Oar PPS POPE DP ENE ORY og LONGO EF a oe ys OF ee oe 
~_ MeN S068 Ouse hee em tesame, Sent tie on wi as Sees htt Be Cie co Te erable) amnion Dubin et eeael a mY EWN 8. POO INT 0 Cae he ay O00 Pens Sa el vip FPG wiF es pemah SA ee Fe #3 ilitntettndel ino-ahd didn oe Oe ra PO 9 eter Be 
= S= Petts Let ee bei beaded ee ee Tat pairs a ee ane re 7 ee tera “Eee ta sae pe one Oe gad eels Digna Or on Ns atten . eg ween OSS) WE got, Cais ew , 8 OEE EP T0808 0° ry UAE ome 8. SLE AF Pw ser yke 
———_—_—— tM, on PRAE hanson the Yet Pagans og ne aie. ash eee 38 papal ee Mtn OPEV en yy oe elegy wr “ye, PP TET ST Peer of ars ‘weue OWE SMR oy . RO! © SF. Seb te OTS FOr oy Prevh » : D Ver po) Ce ee SE OY weny ag iliacheliaMientet vad abd ont ant OF LT PEE oF Cree oe a 
e = = o SeAY OAR, 8 eu erarmre, beatbbatdes eat ee 8 am om oe #0019 AIN 2 oP PUO wv dae * Mean. apten get ane ae Saisie tia ual hoor beret Frew: 670 mime ae Se ep, ee pee WH! Pe e-4 yann-s Ream oe P FN Pe EE rea ay eigiin vy Fe eg VRP PAP e ely rag Se ee, Somewanes sary 
SS poke tend ddl 1, Bp lege pylon gE Bt RT ee ty tel A EP, # eee eee ete ed hoa We Or oh. oH TN ewe My 1g TURN v Bnet s PF PW ei pan pay SOLOS Wee pre. We" 10° d wenger ee EF ROR Myer am acne Ppp os W 0 P06 prety m gig 
OQ = 5% TE ne ae ase Mg SbF © esr ns” rey Pe ieee eee News. ike Fae ee See an orey oe net 8 ET Dt. oy ad re Boe og - TAILS) OF nny 5 Oe ae SE POPP OP ge tihe>- Gon Fe EE TPM PM EO Oat oon 7 SREB ee “50.9: wary 
=—=—=—=> Aww a eee * i as : ated ITs ladle al eat mew bi, . Linde eT) rire bide #4 ripe Te a ane é vee orie as bod eee Fee ga aay ee Fs wie Oe 6e* + ps " ¥ Pete. mes Pup y+ pw & We “ane ow me y re 
i) ————— 8 eee Ch ed Se? er Ore we et A ree *« [ os wt “PONE gag” ge fe ie OOF te- tm ree THe? Fe we perye. 9 Pg Pigs ae wig @ U4 yore. pe 
0 —_—_— oP env eywregn. tiie 0c pitied dada hat eee trey Te SU! © Sas rete Fae e NLT Oe w OE Hy y FOr Sy siete. 
————— ——— ce AN : vn de CP ote, Paw awe PUFA eS WP Oey Ow F OEE A PP eg Our eS Maen Oe Koy 
——= BRN Oconee oy ae a abe Sent iivhaiiead Sit oll Rae Lal eee mene Oe ROE Ss we Ageag oe: 
—_—— ee « s Cae ace ik ple 6 eres - e ENE Re FES ghey eee Fe 
———E 7 3 ". one ny 7¥e uo 
—_—_——— OEE Peas STURT ANS EManie ey, ~ Ur he ny Hy, apo Sane 
= es TR wk | ar ENON Ses. bee, ad ee -“~; 
oa — AINA W Ns MP OWE 6 Vere rary St et ed tn ee el ee 
——— Ss Ne AEN ETE OTL og ns oy os 
SS 
—————_— 
————— 





tye he 
“O8t7 PE Ue Sy Be “ne pers Fem ey 

“Tere ty Pedy ry ore Peres hse 
¥ SP PLP wy - 


0 FOF eer gow ge ypop 7 ater we 
'. so seep bhhded bd tale aad ol TTP he et roe O98 ge a 
LPRSEOE I OTD As Cet avg iiy'g yn 4 Ne Owe dG tmnt 8 ye'gn ig? CFP OMS ERY BR green pty gegen potgli ladicettl ate 1d, LY an 
SEP 18 oe AG a oo Yiye! 5 898 Fe fpr NOE ,. abn er anere 
OP PP Np Cre ng Ve 


PFS LP 8a try Oe ip wre th teeth od ad eT 
le gous Ot A ing OP O88 TRIAS Cy wy 


WN ew og a ie oer oe we 

Fits 

FO Rs grep ge a*2 

CFE ON fs Pe 

PAF ULE AY 8 TEE Por z Beste anf : sawnons Stee sa etal 

tse my cee 0 ONS B= 1hOr EN ae COUR HE otgiy fated cad Bat} : ~ ” ? _ 
OH mH eR. WA rams rpe Oe Nie Aerts ONES I-65 

= = EW Ug) ©) PO"§ quetee cS ww Boe oy Nitateag ed Lh 
eo owe smn 


FON 9 FEE ® bere 
Premteay % 


C8 ON POPS some 





mW 
| 
) 


i 


3 2768 000 67360 2 


Fee Pee gon, Cp 6 ON C8 PE rly og fog 
£8 Ps 0 we i ed oe oF lye Mg ye 
Oe Tong ~ wiagwe® 96 - 







ae =n WEIRTON Yo * 

2 EE WES 9 I eet Ore Oye 

edited Matin TL 

PEP HE ee se 

1 ree weg 
Re ese ew 

tl A 


— 

















| 
| 


ii 


OOF £8 FF oem Petree 
SOP D Fyne wee eo. 

= “ate POPP aire ¥.. Ga yo ge 

08 8 Pwr piey . “ Pas . ts 

YO POR ET he a psy 7 rong aa | 

"ees ey ux 5 908 PEE eu Pue ieee sa 


he ed ee ed 
VON rey 


ih a 
vere oy we 
is Fe > 


VAT ae Orv ee ey oe 








OOP FES ptm map 
POE FO Epa Ath shat Ad id dea =a 
TF 98 Ka wee 08 w DPE! we Oe 2 Ue PF of Fees wer ee OD FE FO ne oa? og 
EF FM oF Py yg POO Nt toy Mier air DOM eer ely Fw ee ae ge poge 
2 OD per Ow PP PP cy gt Wie wry 
puttlined it at LY 
* = tree FOr. ay eS ys, ee y~ 
Ohi <a "eg 4 Ne OE Pe re rey WO ee BIE a Fg 
sesh a ta ¥ re fo - "ON Pe Fee eo. ' 
eee ee Le eteres 5 . a 
: “ey ee 
ee =ty oe, 


Nbiidehdteton a ee ital 
OTS! Be Oe Pg © aw we Ba ae 


a ll a, 
PON “hae 


FEU POLY 0 UW Khe irs begets inact sd Pat as 4A re Te 
he of, VOL PNT Rete anes oy 89. OUP ote des a eee hdd Dele ae ie 
8. ebrhatets Mathilde tae ee Te bape Sled CP ee one Wish iret ss og 
a ek ok ee] es “Was waters wR Ties wet Ta » s C 

le ~ RO epee Oo 7 va eyes 
Om 28 OR EH 4 La onl ¥ 








a Pe bape yt yogis 
PP 20S Pa gt ge we 
8° * OO tap tow gin 
RIDGE FBG pm OE ine * 
ee ee we po F ne 4 OF 78 40 wweay=ge Ge 
AE CLP hg a eee ox 
8 FS we PB ome 
wie eo e a'ge Me B96, 


' 
| 


DUDLEY KNOX LIBRARY 


| 





¥ 


+ 0 PE saree a ree 
OS SEITE: Ga, 














OF EFT TOPO wy ee 
We et Peg dw we A 


OORT ee 
Star os VF eee 6° ot.g ripe p* 
the Mids et tn oo 
EO are Yl eng «gp OP AO Le Fon vipers 


OOP PMN EE ge sige +: re a 








l 





« 








Oe Ge 2 pS pawep OO pags Fit TOF 7 oF 5 hme 
Soe we rig igang FUP se Sy Ale ae Page @ ‘ety OY PEF oP orgs By ipa po am 
a? pe” orga Wepeete -gea er yp cow Pome ee te PP Or ows eS ae go Oe SPR © a8 omy we Fie gt ge, 
"0 SE ot ows Parle Pry ‘ss = or Fe & wit e Vir of ee st FRY Fri “gage 
CRONE ow pay AEE PE cel ad heee ay F000 Py oe 8 yee ge om 
¢ 00 AOS one gre Opp Cer carats s fF Mave 
el Per ay ay PIF eth, he, bP agiey 5 
to UN 0 he ite ay bee Pe oe tty, 
Wess 6 «wey WONG Nas yeaa, f 
———— 
—————— 


tat Ltr age 





fe wns 
ell hd did bd Bali eee a PITA © 





The coupled three degree of fre 


thesS493665 


"a tat Rg ote ny 





Fear er gg 
TF Py «ye © at pee p22 PP a FP. te 2 gre 
halk yhieils dda ad tak Ae ee dw oS nee a aa 
AIS Cd a” FO Gel eg at ate ee oe ye. 

Oo TS AP wey FOR © lh Aled del Le 
Oe UN, : Or ebe ry Wee ome gt ss 08 Fg 8 Migs 88.9 he erp 

le eg wy ty A OA NEP Py 6 e 

PUES egg A 


aw 
UPS We Pef Talc, 48 rin ws 

S OV" ee NE 00 Fees Bere ge 
es aol TT a | FEU rd ogni ves. ~ 7 
Fe sees el ee er ae 


itt ol all a) 
PF SPE TE Tee 8 go y- 





Oe ee TF EO eg a ogy ae 
PUPS gee owe 






















































“* ate Ware ore rere Ht go ld 
COE POP Rae eg ao wee oP ORP Fee eo Fs ee Ore es ae: wee Orgs! «ge 
ON FOe SM” 6-88 CoP O% Fe Prey we ye re OP? 2190 8 in We Hr Be pay FE 8 FP RNe 0 Pat atin Oe Fo wt, 
AF $e Pye YF WE Pip 8 Re ge oF My oe Or por edad teh de ot ht FORBES Bag gt ge og 6? PY He mn sidinhahee od aha rrr 
e fegry SEW 08 Or ote gf Ladi det et Be sone gm ‘~~ orn s PMA E EF £6 OF BO Peernony 7 
Vw ovr e U CS OOP Teme tive y pied od OF et won . PSP POS! at gw rg iy! 2-9°a Wipte 
Wad O56” augers UF os Feb core Fs - 5 CRT en Tee FOE EO UE Os SF ree Sees « ee Fat ie 
BA ete TG ee ee Tar Ode meee te: * a¥ve Were» ON OP wp geEr eT kewego co ee a oo IuOt BTU BD Panegue 9 emde toe get 

WONT. he ew Ow ee ee, «Pewmty oo VP ge tag Tay WE we ee pa a Ere PO fey Pew eye BOE LARP OF & ‘ge Ce ee Se 8 ws pee We wer 

thee eee piesa eared a Se erst PIE psa a 5 = NT OTP Cty 6 Ory FO 8-0 NO & AAs oe Pe wwe OO ae PAY Peres arene tg ore 6 eae 

wr Amengy Beg” Petite cel i on EO. POE Oy Ue PE TRCINY FUT tas heresy Fo OF BORE png: ie Oe ENS 8 PRD 0 pap gripe gs hc ee oat ee FERS? wm BE ph gt OTe Po SOK © 

smn ar . nae eee ety atiaiaed rw Ot Ne FES rie y 4 TE OF an eon are aT ay FE 8 Trang 14 SOOT FO & Wegee wege ecg LV AAP PWT O wey one ato or Re 
, ‘ We BY Pe Owe. wwe wr. oy @ EP set oo yy ge ee es NINE oe eres gig ES PW © keg ¥ BE” epCTe ages: co e's FP ee wrersverwsy SMT p= ee Hee ors wert tye iinet gy 
cuenta ’ Wer re nwrere € wer “eo. we "t= Te Pe | . hy Veoh WW5r pap .. wrete «e% WIE SP ym ay ge Perper wy ae te ol ad FO Wia-e 4 Pe ed Pe ey er Fg r Le deh Meena Catena Niel ae at Pee oe ee ade et ae BO FI wR a oe 
My Eo" Oe Tee Oe Ot i: ae boy NORTE 0 80 Uy: aA ait Nelch die abet ONE Play tC a eae 588 PUNE gg =n WORE opin dak dey Meath a af eee SRT Rie oy 6 gta 0-90 Or wy -99 pm Ome Pa 
papel oo STR atte bad gee e SENS NOW Ue ycnartgs Dadi tal tad pte nal lea oI etd SOR TECY GUL Ce vy! ote irre ing INO SOUT og) BH BEE PO te ae Wee wee Sw rr yang 
iia alte ads ih ek te) iil Laie a Ee Boe VM Hees, I Pr Oe Petree TNS OEM er arse uk Prag omg - iia denne a eh PP 084 0 rg: eT y sp Fb: Sg Ererenel grRe gigs. 
& thes okt a POO OE HE 89 IHS Gy- . Wa TEN ob whee OP Porm KAS OP FP Pig 8 OPP Err OA" Me Pr rey Tf ote POO 1? Fw gripe FOP OF F0- at. weeny. 
ue Veo We Ee, et ee HHO Tie te ie, Oey SONI BP Cre, of SOW TER vping’«, ty Fae Foe Us may om call Sel ae tn ne Se-thicttty Dell dt ak oa Pete dugg 
SON SEE EET Lane, OUENTE BRET OEP ri" oy FH oat ron Ne ret ci Ce Due be aw J " VOT 6 wn eran pug £0 ES 8 Ng een ne 16 ue 
tr ey, Perey an Ol ee Ree betel th hoe SL te “Wrert ig. PENwT POP We aly ibd fea 1 Nile hy, * EPP 9 P< Wao Way Bea, yen 
re SaaS PN LR y Ae Lk ee ee , : VEIT Fk ae, hd ta dal ten a ONE Ee. 9 aan oy 0h wewWer OL hey oot 
itbtbubste thane Ee eA ee RP eg eR ea TE re <w wy PRS Ele RIO EL Oy, VATE Wy epee 
eee ane Cedidned oh pet = Pee sT ee” “ob 8 
ae NONE ey 0 
TE Re ok, & 
FENTON: 02, ohewaeny, 
~ PUA ETS Yemen, yy 
Wie eA Uy Brae; yw, 


Fenn Ae EE CRIN BA oF weeny PT ew Pe Peele wy 

«46 2ag~ PPO Fe pt one we OES A 8 aes 5 ee ap PPT HN IM oy Hey 

EN NTE ONE Ne Wey" 09 MINOT. 6? mye. gree ue BPE PG EVER Grigg PEIN For BN O° ara 9? ge 
PD Ue oe SARL My ara he PEE god ae” YN (Mul ag 6 wee ¢ OCF Ww! viel, 

Py Poy pe VPP Ferre. 5 OF 8° Ve Pet werine og ower 

FAW Dg Farciwy SALTS bt ee OM 6? aT be 


POE FSO RR nemgre Ve ae Weimar ey geang's (54 Pee oe 
we CLD ENO PLO BIE KO 
Lid ded She ae DEF we irnly ri ray Pies gig 
Sdcdts Moot tae Te 


Ve som pores 
J HE 5 6 ape bv eray ce 

GPS TAG So eer re Fe Pn « 

Oy Gg ny aya 


wee Are ioe Fe a, 
Pe rh a ft 8 RUNS IW + ge... 
Wey eats 


oe aed EPG Rey, 
itt Den re ke ee 





#90 Sores 


bbabdbnithetiien et dehed ee biden te ee 


SL 98 es 0 ee 














SIRE Ot Fg" Pe Po ae Boe we VE 62 MEIN Pele 4. 
10 BS ES 88 Boa 9g gare aewreyr Eve bet oe Wet ot ge 
SPP Par wrote ONES rere 
EE ES em eghag? cy? rw re Wise y ov 
Zap WHF "an OE WE toy Pir Pot ww ery PEP am ye st Ml Akad bel 1 
ae le Om" ttl + oe Sep of OPV er 0 ore et an ee oe 
- "AEP ee Ge ape OW Ee Le ratty» rr ered § ¢- are PIPES Pweg Rode! scgee 
we oe Oey 3 4.62 ve Ca ae er) Meote® bay 2° ward Stee or Braengy 2 fe 
eNe 0 °C qe, Be De ae pyre wie NEUSE Ow TWH ayia. 9 PEW eyes. vo anette s marae Fo Pitas nit pol Mal ee 
\ : Sea 1 Nae ett Of -Fivervey 16 bbl LO gk va tes oP ee OW Of oor ween, Nee ti tamng ” 
; - art OTT UW RN wre 94 bran te PRM Aye Wee Pie a. 2 Si ae rman Sie dle ted ed aa ae 
= tts ; . ; OWES 5 are myer ls a raw @ me a ee. vy « . was 
botched eta Od . ON Sle Uae. 4, 35 on : PHO RA PORT y- titers stot ithe eta rey SA te Le : 
Lae ett tetehia hieci e Ltt ins mae ee UTEP w ir wie ripee OMAR Maen wi are ie Yiralrg Sener WNW neve Sia ete, ‘WEE we ee - 
co eee OA eeylieintohy tim ial se shill tested he ha) Peten, Piha tet ea hat ee) BNE Mer ey, a> re “Ms wefan sey, 
OE 8 Rs trey 6H, VWE Ra wy Terr rast o re ore . Par ps : Pe aa 
2 uns ie AFUE Oe on ent Duaietre j tia bet reyes SPs te aareue 
© 09 Nae neg RW oe Ore: : 
OTN BD Spy Fe rack 





OO Ode, 7 
oe AGG USI a sa tu 
Aces ere 

Witty AT eo 


SW OMe 8 Yttgn-ns TA aaa” FF aOR 948g” 
OF Tie ee sug ON Pager 

PO oF as 0 ov 
OPW Se 





O- Wetne Pitorae 


28 ie? SOP RI Rs gy wees 
8 He ye VE AY THU POT, ot 


PUPP aE Wialtg dee A ot 
Cet De at 
















ee gern bhi lialiaigh daha bia Leathe ik) #- 
NO Senge gh Fug tid det ie dated on oo 88 EM & pe Tbre Gym 
3 Er Par igo IMO Gtr 0g ots lb al det LP a od ek TAS SEY PED 8 we ep yew ba 
We 7 0° Aly wy jun dir dh illicdeed Le Lee Uy NEE 8D EB sop ang PF OOO ST 8 wore PMT PEG gs Bish en gp Gene Pern © 
a Cw Pe IT & ore Pye OP ER ig ry EEF OOS are gw gps peey- EPR” 8° 28 iS” Pee 190 po pee wiry Wes Parente Pr alive mi 
MY eee 69) aay Sitesi ds he = WEP OW arm: 5 ry PENT 06 oe gine Ure re ra eerie wrwtye 
C OVE Prem 9 5 oy hoe we. ee te ao ey Fore sry 8 ou trn: 
pert ‘ Ete ee at heidiel Bk me ba We Tyree ay ary, Feels OOo 6h We yeu (hy ue, a *. ff sae, Ky 
adeeb ‘7 dal a | . .- Wein i559 peer en ~ PA Lene "ie 
ee RANE Keer ee rm ORE a. ae < j : Pe Uh es tN ot oF LEAD 00a hy 
Se Wr wr wey. 46 PO OU . ra 
4 ‘ ah We Oe UNE WW emee en 
SEES BUCurt. wees Mee gin! ’ 2 
helmeted sth tbh kesertiee 





O Fie s8 so Pie u_y, 
VEAP a" 9 gras oywgtes 

FSF SRE. Shey spe igo 

Mil dele ak ad? he x 

te Res PE ary O WEA ee ow 5 









Fer peed 1° PF 6 pee a-gou Se gpg 
FOR" Pls El WH eso PEt SKY oy gory SUP re 0 Ue ee rE wee 
POL yma eo Ww fn py ipa We! £09k" of a-g oss eee 's~ OF POW Bet 0 80s ote po yell 
OSES “o-dbr ona ee 02 OF Bee ele el ka La yo) eng pee 
59 Sb OS AD gt 
it sleds i Pt EB 
wee 


oe oe" Narr Behe ~ 
Veh g ee 688 = 7) BAS alee ra Hin = has Bp tx gs Xe Fos " 
bites a ot Oe ae eect shes BY yaumea., oe 

Oe EL 8 OP U4 Aaya AS 0 VCU's WAN » 9) we 
8 LE IEY greg, an ore ree. Oe ETON ee tO ve erie aie o'ey. bhandliehchdndahall deta F PPB oe PEW pew © 
wt F en 8 "O99" a fe me age iy? 4 wer Rifai’ iad on oe te P88 Oo a 9 -< BAO Vise Oyen MaVyye -yryoynga De ee 
£ ONE OS re Irgen sbi ee MOLE EV Ne T-8* SA Bape KV ea oep ¢: 

ee Aye ig "FW Meds aig | EF Le o Fav yy sn wipe 
Ca a” 6 ate rine Cee 





WORE Oey 
he ate 





VOW oh gh gy VTP 
MENS Pe 9 ayn ee og esy 
eee ialideediee belt itt a Da 1 


oy. ‘Pe - WSO LP Parr on Oy By 
OS ee ry Wave, 


ee PY) 
OP ce ire a) ee Poe ot 
OO wr erin erm OF ele 
Lelie ad Sd 

















Wren ye: 
























OY Ope gg FTE 
OE TOS War ade We wna p SO IP apse rer yornts 
reer ee ee ON as BUTTE gay BOON O° F-98CUr 0: Be yeRryts Pee Te Te" 0 VU 9 3905013 O-#- 
Pe hdatth a ih ata ee ey AT Er ey ey BS 8 Fat 18 VON F 8 ar ve: pig: Oe Se eee 
"9" es Re Ree PRR Fat - ON 8 PUN @ Warege ga, . 
2 STNG LA SMe AT mpeg iy WO NSD OF ee We On eRe WES 4" KER" gs lh Se Seg? 
Oe ere Ty a Wud \awe oy POP IS oe 5 ge shew, Se Pvatcglnge'g ave Swe ENPMPOR- vem gece. 9 palette be 
ntti devonaetee ‘ eich eee . Re Gig é bas fi OO FOP 0°08? hig H.. PEP 2" OND ea Oe Pe =0 a” fe wry: 
= My « * wl ie “". or ns 7 ’ ro’ me — bak | » i" ‘ 
Pt lel hte TD) WR 6s Were ase, “ oot. 95 Da te Tl oe = FeV 0 EPR re phad TEE ee Py ie Si ne ache gal Pe phe eee 
NAV CLT ent e couarnes oP Wie 88 Wb ete Ou N Fone he THY A ee wey ie Pt dee Cal rt ae . 
FUT OAL oe, et tes a he ee 78 NYT AO are PO ei 9 Eps Oy CoN E Teo 1m 42, 
oS ost VES Ler A 5 e Sak Lt Tae ERM oR, at thr verso t P64 oA Ww 
bpblra ti, tat bade tall Molle PL PY vere, 
Pil he Sol OT | Leh ti tant 18 T 








Oe Vp gee 
Se orig a gas ie 
Oa by oT tees 9 fe | 





ITEP es Gley- Wegwae Li lee} or 
Pe OEE Sighs, garg 
©) OCT we ety 
oe avs @ eth oad nc 


7 Deewana) shear al 



































































































eA ic ae aT oe vi 
PEE 08 Te cerne EIU tay wae CHP oe eT em A FU es orp ee OIRO TF ie a eas WY wages 7 *OTES" 9 ar wa gy 
> PM alld ott oe oe teas | BPW FS aes Meg PMs ay BE ry bell eae dihatait i ox, tie hs CRrwte Fretacy # a-vey fi Sixveny ris TO Tre epe wea we OP Fee a! EER 
midds BE ERTS NI sercy OF PR oe oe tia oe FO OE EMS PRE vb et ig Fy -9: <8 Fi BEETS te TH gece SY 3° eens 5 8 horn: Bee Pree Peis, SPP EE ENE! 9 F Cte gee y 6y 
> a ee het lal es vite Was” = 0a Paring wang 4 we AS eee Oe eae y VP PY gy POG Fyne -gee idide te i eae 2 Va Sa- aor we, eee st. lp Tavs 2 rep ae ote ar eer e ee 
‘ea hes: > ahs Se. Sm wer oF Ow 2 oy HEX, rey, Tee Wh neh Re ye Peer ey Wee) eas Le - L" 8 @ive%ig” Penta. a: To pes # 2s & Hy yoy, ee ae Le 2 Bee gn Whe" Serarer te wt yi gee OF 0° Dew WY, ° i - re 
‘Ow eore OWT AT Hope my, ope eee iterate hd stadt, reeghe EY Me is prot een ye Maree te SMa RSET 5 Lite No FEE ep dE 8 otteys 5 10 mae Fee Vy he F Pye 114 g po fad (tintin a MER YG » 9° Forte! a rte 829 swear gig 58D PNA Prepegniegre L demdhtieded dd ade at aa 
Rd 0 Fadl gue rR at 45, : ee peste me ise ; Sabie llr ay ee by ele ate a ay Lat ren POE” Rae ge Avy EE TPT Me Eres: 6 Hom OO PO PEG 46 EF my, SO" FP Ey i g- ae gs Sie FV yP ot ps gee batndel 3d Se Bees he enews a oF STW pte 
ON eae sbede Scien pis pte Nin dls Pa Th Oe POT LS vee SOE PAW One prep TL rag SOE Ke OORT eT Eg | ign PEW IW F208 Sep gid POCO A EH 8s Pie ct ote es Woay: indent iiade aitlvdes om ans od at ne 
ested hic few Voted AB ie road POCO UNC y rare Ww we VO eg: “We ower re ED eee ga5 pegilis VET I Fiz OEE TEE ME F989 Ore ON yg wp co ota PE PDO UK, gee ay 98% - aH yey 
Adal ddedopin tt Abt, Hib eba Leienh ee ered Serna veo aah el ys bibs xed L sal JN Ga ate Ware = : Eres y= Ts Pe ee trey 0 wy Ldtilalihdahieds iomt tell oer el oa PAM wt hazy 696-8 wey Cet ie Lt hdlhdepel ad ae he ee 
arenes UTA SUNS wc WUN Cannone bide dhlin igs drth pete at beth Mant hada tet erg UE UM Weny remary bah ho ae ae io . ; EF 8 RMT Ber 8 0 8 -y ae | LPR A Fearn Wises Lenn Ue O08 Ege Oa w eg gry oe PETES 0 wer StS Ey ape 5 
et aanvedr ulate (tat eT Aealet ME Da on Tn Tes : ; ¥ é PF EONS PVE STE 186 ow yee cyry 2 Pema Fijeirsaoy? gq yige Rhee adsl on ale oe i pice a eee es 
WDlecdes ob Joh ten err Uk tiene aan) eho re ‘ atl “4 .” 2 . Lore ty 0° WP VO fe go ghey Kn, Ne Gh ashe My oy Oh Fe eee res Fen ae ages I~ 9 if fee pew 8 eet S Ny yee 
cee co neice ean rive e biktals v9 ee atei : WH So Gy SEE Fey & F- Pee we oS TU Gas eo 2° SPF Fao EG ares rey SOE. WPT Bs OE 90 Fw: 5s TOO EF ORE e Vigt 
i Pe 8 v8 Reon “ Aha) pts Di 1 Pegien, “OP wet hee - pial Mowe ile Th bee ie fo 216 SRT A gy ye "AE ys Pile PO yes ea nyey V5 84 8 Eee BLA Wa BH Gh gtijey 5° ener os ORI OTD 0 eB ew we Y WF are lewigiry 
Sad eee ETH, YS: pret A IB Tl italy gee ee ek a Dipl Ni nel Gnd Sa PO a ae Coa ee b eNO SRM Wy: sayy gree oe yy EW. 0-9-0" Paper e amen segm§-- = G 
- < MtSHe sana baer Sues ettgoes py f ‘ Ao etait ee eee plas He & tte Blade ety ~@ PUM TES Bees 5 OP lot oe EA rOert ep ay . TOLTs #00 Pa ge BY * 4° GS ve oF AOE 6 temy 
Pet iri ge Wragg 0 at td cal al @ aT Ay ala my ri : TIVE. Hts sy. “ayes y "vA St wg, Sy or aha sty OP 2 ft Peer 
avrutaceee wee sey Spa: ere se ETRE M wr aim ba ie Rela aa ae R heh eels lacs ats Pak oe LEN AY Serer cen twit 
SUNNY “Ou Diy. “aay eA oa eas © EVA EPPO Eur aiAD when het Se terPetattce-atny ne oR o “eben EES 9 6 Pye. aNeeg 
. : J % OW oot n, d * Rel ia Bl ae ie 2a aa “Ee wr 
blah Sabeeineatea SPA FOF Wey - He FOTN. 09 078 ace 
baal ree LEICA TE Newel ver En SUP, OF Lg WO perl. eae ree Y i 
yb dabeeloebteehante taster tobias "OTe, Po oer 
bleep ed eae Tee PRY: wryly PRE Mv ye Bq ad eg “hh 
Rie, Alpe, od ey te eG 








_ Pes Crlwerey reg 
wes Fe YS BVI OW 6s yore Cogs 
OF grew» e+ 08 Bh a4 gy ‘stir Wop .. 
OE WB By y 
Whlawd 









- O's wen we opt) 
o Feng ee FP BUR: POOVELY™ Ng aie O° 999 tai 
SS 'oe" ysis) gra where pthedietatie ch ie dak at an CPW eo serrate, FO erin > Sed hetiet to ka 
o oo RE VE mg op Wires bits Mebuiels del det oe ci ak cat na a O° BE VTD AO IS e* Tis! OT Bone a FINE OL ORe ee eee sh 
r EE OE 8 Se a Piers Tet ryey s+ yw, SW 2 a 80 Seep eyens bed Sr ede eae ee et PUNTER OE Yl Fir evs: yp are 
& PE awa” oman Vevey 0 FOF fh 0 eg ieee i ie Bd he Fee Fr blew FS wig Retell Aeinell athdek bl de ht ae bh dead O° F2 © 4 wet by 
| Pes y oF Bees ba ec oC Pr rsa-e" o 2 BO era F a taee OPE pe poy wren # 
‘h, deve: ae SU ate, Ratio na FOF 83 88 9 C8 HF gig ge Co ule ‘ 
+ ie iee ae eS peut re yee LEE Be WV a tea Sag 
: ~ fe Ne. ST At, Slr F 
‘ Lie deat 1 whys 
* 4 Fos FUP ery ay at ay Ve tf ee Ca 
ie eee hie bit Lar | 


TTS Tease 
‘eM. BTL Pee my, be 





aig oh? Cr oly 





0 TE ew FE Pye 





ON PFE ge PP ees 
3 tein? 
ss Pl “oD Oe OP "Eee 4 

















OP Wry peti Mndekon bene ae edith bt it ak a ee Th 
Lee Le Ste idk ae Bites ery re 9X 
1 SEV APT OO > EE SF PLT e es Vg Gy ii ee ee ay SO OF RAS 2-9 ee 
Oe YAS SEE BY eee eige Cte or iis SUOB 8" Figg 23728" Fhe 5 
= Se ee SPOR 9 ie Mies OF OG U8 tow Ba ele awe Gey 
Mary =a STAR bee ww gare fet aid ee Tee Poe wee E TS ven Try * a) ae ne ol 
P -¥ Gore. eles Hy RG ave 0) ed ee TPA 8 > oes ee 
MeO katt OS by aor ee wis ANE STW Oy vpn Be eR as to aWapeg oy, + Hated wip oe Facts THAI > he 
OI 6 oF CMY Oe 8b ¥ “al oF EEO 2% sang TE Fils 
rare Syn TOR. Mey, cd 
a) red a ee 


ee AG o ae" Feat pnetenee s nn le 
ee Rah A to 
8 OEE fe ee an ge 
» Yee <* “Sur ante 
MP WO ogres oy ch % — ian 
+ vow ey Byes - - 
Wy the ye Woe Uta oe 





